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ABSTRACT 

This study examines the possible use of hydrogen and 
oxygen isotope records preserved in tree rings as paleo- 
temperature indicators. The study consists of two parts. 
The first part deals with the methodology whereby the 
isotopic composition of tree rings can be determined. The 
second part investigates the significance of the isotope 
records as a paleo-thermometer. | 

A new method of measuring both D/H and /}°%0/!* ratios 
in organic compounds is discussed. This technique involves 
the pyrolysis of organic compounds in a nickel reaction 
vessel and utilizes diffusion of hydrogen gas through the 
thin walls of a nickel reaction vessel. This method 
enables quantitative extraction of hydrogen and oxygen 
from tree rings for isotopic analysis by mass spectrometry. 
The method gives reproducible and accurate results with a 
variety of organic compounds and standard waters. 

In order to eliminate any fractionation problems 
associated with chemical heterogeneity in tree rings, 
cellulose was chosen as the working material. For hydrogen 
isotopic analysis, the cellulose was nitrated to eliminate 
the exchangeable OH hydrogen in the cellulose. Two methods 
of nitration were tested. One method (W-C-CN) involves 
the extraction of a-cellulose from wood mill followed by 
nitration. The other method (W-CN) involves direct 


Nitration of the wood mill. “The first method (W-C-CN) 
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appears preferable to the latter (W-CN) in terms of repro- 
ducibility and correlation with climatic parameters. 

The D/H ratios of cellulose nitrate extracted from 
the wood of three spruce trees, which grew in Edmonton, 
were measured. All three trees showed a sudden shift in 
D/H ratios between juvenile and mature sections of the 
tree. In each case the juvenile sections are depleted in 
deuterium compared to the mature sections. The phenomenon 
is not present in the oxygen isotope records, nor is it 
due to climatic changes. Possible explanations for this 
Phenomenon are given. 

Both D/H and !°0/!°0O ratios from cellulose in tree 
rings show functional relationships with climatic para- 
meters. D/H ratios in cellulose nitrate show a strong 
linear correlation with temperature at the growth site. 
120/* °Ooratios intcellulose tare’ less welllcorrelated with 
temperature but are still significant. 

Strong positive linear correlations between isotopic 
compositions and ihecumcenserature? together with nega- 
tive correlations with the amount of precipitation, 
indicate that the isotopic composition of soil water is 
the primary factor controlling isotopic compositions of 
tree rings. 

Both hydrogen and oxygen isotope models are developed. 
The hydrogen isotope model incorporates the evapotranspira- 
tion model of Dongmann et al (20) in which no temperature- 


dependent fractionation is assumed. The oxygen model 
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assumes that oxygen in cellulose is derived from atmos- 


Pheric CO, which is in complete equilibrium with leaf 


2 
water and contains a temperature-dependent fractionation 
term between CO. and water. The models developed here 
are in excellent agreement with the experimental results 
obtained with tree ring samples from various localities 
of North America. 

The study shows that while the isotopic enrichment 
of cellulose with respect to soil water is a function 
of relative humidity, paleo-humidity information cannot 
be obtained from isotope data because paleo-isotopic 
compositions of waters are unknown. The best climatic 
information from isotopic records in tree rings aeeesns 
to be temperature variation, particularly that of mean 
annual temperature, which does not require a priori 
knowledge of isotopic composition of soil water. 

In summary, this study shows that isotopic composi- 
tions, particularly those of hydrogen, preserved in tree 


rings, provide a measurement of paleo-temperature changes 


occurring during the life span of the tree. 
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CHAPTER I 
INTRODUCTION 

It is generally accepted that global climate has 
undergone considerable variation in the past. It is to be 
expected that Such variation will continue in the future. 
eoreider ite ene vital impact of climate on human life and 
the uncertainties associated with future climatic variations, 
the development of an ability to forecast these future 
variations is an important and challenging task. In 
addition to natural variations which have been observed, 
there is a growing belief that human activities, especially 
recent industrial developments may have considerable influ- 
ence on future climate (7/5). 

Study of the past climate variation is essential for 
the development of climate models which may form the basis 
for prediction of future climate trends. Over the past few 
centuries, the development of meteorological instruments 
has enabled the collection of data which has allowed 
quantitative analysis of climate and climate variations. 
The study of climatic variations, however, cannot depend 
entirely on instrumental data which extend back only a 
relatively short time period. Additional climatic informa- 
tion must come from proxy records, such as fossil pollen, 
mountain glaciers, sediments from lakes and oceans, and 
isotope data. 

As early as 1948, Urey (99) proposed that paleo- 


temperatures of the ancient oceans could be estimated by 
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measuring oxygen isotope distribution between calcium car- 
bonate and water. The method is based on the observation 
that-thevoxygen isotopic composition of Caco. differs from 
that of water at equilibrium and that the difference is 
temperature-dependent. Since Urey's proposal many 
researchers throughout the world have used this principle 
to investigate paleoclimate variations. The study of 
paleotemperatures in the Pleistocene was pioneered by 
Emiliani (23,24,25). By measuring the oxygen isotopic 
composition of planktonic foraminifera from deep-sea 
sediments, Emiliani assembled a detailed record of tempera- 
ture variations ae the Pleistocene epoch. However, 
Shackleton et al (89) demonstrated that the observed 
isotope variations reflect primarily the changing volume 
of polar ice. The generalized paleotemperature curve 

thus revised by Emiliani and Shackleton (26) revealed 
recurring temperature fluctuations extending for about 
700,000 years into the past. 

Another successful application of isotope techniques 
to paleotemperature reconstruction has been reported by 
Dansgaard et al (17) and Lorius et al (71). They showed 
that the oxygen isotopic composition in each accumulating 
layer of polar ice caps is closely related to the tempera- 
ture at the site of precipitation. Isotope records from 
the ice caps are in excellent agreement with those of 
calcium carbonate from deep-sea sediments for the equiva- 


lent periods. Thus, they yield temperature information 
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tracing back a few hundred thousand years. 

Bio-organic materials that are deposited year after 
year have been reported to preserve climatic information 
in the form of stable isotope abundances (65). In general, 
D/H and '°0/'*O ratios in modern meteoric water depend upon 
the temperature at the site of precipitation (16). Thus 
the isotopic composition of biological systems utilizing 
mMeteoric water may preserve the record of past temperature 
variations. In addition, if the photosynthetic process 
leads to temperature-dependent isotopic fractionation 
between different photosynthates in the biological system, 
a detailed record of temperature variations covering 
several thousand years may be obtained. In particular, 
tree rings have great potential as climatic indicators 
for a number of reasons. Dendrochronology provides an 
absolute dating mechanism and relatively long continuous 
climatic records are possible. 

Like any man-made sensing device, natural paleo- 
climatic indicators must be calibrated and their distinctive 
performance characteristics must be understood before 
rigorous paleoclimatic reconstruction is possible. There- 
foreprthe®main*objectivecof this iwork istosinvestigate the 
existence of possible Sserape thermometers in tree ring 


sequences, their significance and calibration. 
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BACKGROUND 


a Isotopes and their Properties 


It is now known that almost all the elements have. 
more than one naturally occurring isotope. Isotopes of an 
element have the same number of protons but slightly diffe- 
rent atomic masses depending upon the varying number of 
neutrons in their nuclei. As a result, their chemical and 
physical properties are not identical. 

The chemical properties of an element are largely 
determined by the number and arrangement of their orbital 
electrons. Thus, the isotopes of an element, having the 
same number of electrons arranged in a similar manner, 
would have similar chemical properties. Research done 
Since 1930, however, has revealed that the isotopes of an 
element and their compounds differ slightly in their 
chemical properties and that the differences in chemical 
properties are due to the differences in their thermo- 
dynamic properties associated with their masses. It was 
also recognized that the isotopes differ in their physical 
properties that are directly related to mass, such as 
diffusion rate, rate of evaporation, vapor pressure, and 
density of gaseous and condensed phases. 

During the early stage of isotope studies it was 
believed that the stable isotopes of elements occur in 
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mass spectrometric methods for the determination of the 
abundance of the stable isotopes, it became evident that 
the isotopic composition of elements, especially those of 
low atomic numbers are variable because their isotopes are 
fractionated in the course of certain chemical and physical 
Processes OCCUrTIng In nature. This fractionation its pro- 
portional to the difference in masses between two isotopes. 
The higher the mass difference between two isotopes, the 
greater the fractionation effect. For this reason stable 
isotope studies have tended to concent rare on commonly 
occurring low atomic mass elements such as hydrogen (D/H), 
Garloont (Cy)? “Cy nitrogen (° N/* 'N)|, oxygen (7.0/0)! and 


Sel pnura tc. S/S). 


Ze Natural Isotope Variations 


In general, isotopic fractionation occurs during 


several kinds of chemical and physical processes (79,32). 


is Isotopic exchange reactions involving the redistri- 
bution of isotopes of an element among different 


molecules containing that element. For example, 


Sit *O, + H,'°O. 


NS) Toe 


1 asics Ieicyune 
5 oy GES 5 H Ox 


2 Unidirectional reactions in which the difference in 
reaction rate between two isotopic atoms or molecules 
leads to fractionation. The differences in reaction 
rates are due to differences in thermodynamic pro- 


perties, such as energy, entropy and the free 
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energy of the species. 
3% Physical processes such as evaporation and condensa- 
tion, melting and crystallization, adsorption and 


desorption in which mass differences come into play. 


se Equilibrium Constant in Isotope Exchange Reactions 


The thermodynamic properties of isotopic molecules 
and the resulting fractionation of isotopes have been 
studied by a number of researchers (5,99). Quantum 
mechanics shows that the energy levels of molecules are 
mass dependent. Thus isotopic species have different 
energy levels which result in a difference in the behaviour 
of the isotopes of an element in certain chemical reactions. 
The fractionation that occurs in such natural processes is 


indicated by a fractionation factor which is defined as 


Gag = Ry/Rp [1] 


where Ra isthe isotopic ratio such as, D/H or ‘20/a Oe 


molecule A or phase A and R, is the same in phase B. 


B 
A typical isotope exchange reaction between two 


molecules A and B may be written as 


pity ana 


1 ume Raed [2] 


2 i 


where subscripts 1 and 2 represent the light and the heavy 
isotopes respectively. Then the equilibrium constant for 


this reaction is given in terms of partition functions 
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AC equi Mibiaum,.the Iractionation factor’ @ is#equal* to the 


equilibrium constant k: 


Lf wche fractionatsaon.tactorm is: greater, thansunity,.tche 
products. of. reaction, 4 2lewild. best avoured ~whide,it dit. is 
less than unity, the reactants will dominate. Equation 
[3] indicates that the equilibrium constant or fractiona- 
tion factor is determined only by the ratios of partition 
functions... The, partition function is.defined as 


-E,/kT 


Q= > e fs)! 
i. 


where k is the Boltzman constant, T is the absolute tem- 
perature and . is the energy at quantum state i (3). 
EqUaCTONS (16) anc. loleinartcate that, fractionation 2s va 
function of the energy of the molecule and the temperature 
of the system. 

The energy of a molecule in a gas can be described 
in terms of the electronic interactions plus translational, 
rotational and ee eal components of the atoms in that 
molecule. Thus the overall energy of the system can be 


written as 
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Accordingly, the overall partition function of the molecule 
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becomes the product “of individual partition functions 


associated with the corresponding energies. 
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elect. trans. 


In general, the partition function due to electronic 
energy does not contribute to isotope fractionation except 
in the case of hydrogen. The excited electronic levels of 
the atoms or molecules of interest lie so far above the 
ground state that only at temperatures of several thousand 
degrees does the term kT become comparable with the energy 
term E. and the partition function assumes significant 
probability. Thus only the ground electronic state of 
degeneracy a contributes to the energy and the electronic 


partition function becomes simply the degeneracy 


When only the ratio of partition functions is considered, 
Lt specomes unity, sO enat no isotopic fractionation occurs: 


The. translational partition function isigiven by 
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where m = atomic mass, a = length of motion, Se = volume 
of gas, n = quantum number, and h = planck constant. 
Replacing the sum by the integral leads to virtually no 
loss of accuracy and the translational partition function 


becomes 
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where M = molecular mass, V = volume of gas. 


The=rotatzena. partition’ function* as" given by 
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where o = symmetry number is 1 for an asymmetric molecule 


and 2 for a symmetric linear molecule, I = moment of inertia 
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and j = rotational quantum number. For T >> 


equation reduces to a simpler form by making a change of 


Vecitaore 2 o=— y()-) poana- az. = (27+1) dy: 
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Phe form of ebhe vibrational wpartition, funckion 


neglecting anharmonicity is 
-U,/2 


Oui =. ) e7 (Wt) hy ./kT =e —_——— [12] 
; y | 


where U. = hv,/kT, Vig vibrational frequency, V = 


vibrational quantum number. 
Thus the overall partition function ratio between 


the light and heavy components of a molecule is given by 
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Tne tuanslacionaiepartition function ae, all temperatures 
is classical as is the rotational partition function at 
normal temperatures except for some molecules conta ning 


hydrogen. According to a theorem of Teller and Redlich (99) 


[14] 
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The ratio of the masses of the isotopes cancels out in any 
chemical reaction (5). The ratio of the symmetry numbers, 
regardless of its value, will not lead to isotopic frac- 
tionation since it merely represents the relative probabi- 
lities of forming symmetrical and unsymmetrical molecules. 
Thus the classical ratio COE is omitted in calcula- 
tion of the equilibrium constant. This omission leads to 
the equilibrium constant being the ratio of the equilibrium 


constants for the dissociation of the two isotopic molecules 


into atoms. A new partition function ratio is defined as 
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making whe natural Jogarithm of both sides, of equation 


[17] and utilizing series expansion yields: 


aoe haa 
! 1 = Fe hom Sos : = 
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thus) thefana iettorm of theifractionation® factor asioqiven 
by 


a ty gaye eee Dy sau, \ ~ ) GAU,',- [19] 
ea / " yf 


This equation indicates that the fractionation factor is a 
function only of the vibrational frequencies of the mole- 
cules involved and the temperature of the system. However, 
in an isotopic exchange reaction, the vibrational frequencies 
of the molecules involved are fixed. Thus the temperature 
of the system becomes the only variable affecting isotopic 
fractionation between the coexisting molecules in the 
reaction.» This is the basic principle underlying the iso- 


tope thermometry. 
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frequencies and low temperature (U,>20), the function G 
approaches the value 1/2 and ln Q5/2, iS DLOpOreToOnal.. co 
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frequencies and high temperature (U,<3), the function G has 


177. | Consequently, Ina is’ a function of 171s For Mow 


a slope 1/12 and ina, LS proportional, to Lye At even 


B 
higher temperature the fractionation factor becomes unity 
and no isotope fraction occurs. 

The theoretical considerations presented here were 
developed for ideal gases and extension of them to the 
condensed phase is complicated by additional factors such 
as lattice vibrational frequencies and interactions between 
the lattice and the internal vibrations. However, this 
theory serves as an excellent guide and similar temperature 


dependences of fractionation factors have been observed in 


both laboratory and natural systems. 


An Chemical Processes 


In the previous section, it was shown that vibrational 
motion of molecules is responsible for isotopic fractiona- 
tion in chemical reactions. However, the vibrational 
motion of polyatomic molecules is complex. Since the 
normal coordinates of a polyatomic molecule resemble inde- 
pendent vibrations, the study of simple harmonic motion of 
a diatomic molecule provides an adequate background for the 
understanding of more complex systems. The classical 


vibrational frequency for the simple harmonic oscillator is 
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[20] 


m.m 
reduced mass —— andk = force constant. 
m, +m. 


The force constant K is an approximate indication of the 


where u 


Strength -of* the “chemical ‘bond (85) Therefore, the iso- 
topic properties of a molecule depend largely upon the 
nature of the chemical bonds within the molecule and the 
reduced mass of the vibrating system in the molecule. Thus 
when isotopic fractionation between the molecules is con- 
Sidered, factors such as oxidation state, ionic charge, 
atomic mass and electronic configuration of the elements 
to which the isotope is bonded play an important role (79). 
In general, bonds to ions with high ionic potential 
and low atomic mass are associated with high vibrational 
frequencies and therefore have a tendency to incorporate 
the heavy isotope preferentially in order to lower the free 
energy of the system. For example, in natural equilibrium 
assemblages, quartz (Si0.) is always the most !°0O-rich 
mineral and magnetite (Fe,0,) is the most !%O-depleted one 
(11,4). The small highly charged Seer ion bonded to oxygen 
has a higher potential energy than does the relatively large 
moe ion and therefore incorporates !°O preferentially. 
Taylor (94) found a relatively large '°O difference 
(1.68°/,,) between grossularite (Ca,Al, (Si0,) 3) and andra- 
dite (Ca,Fe, (Si0,) 3) at 600°C. Since both Al and Fe are 
in the same oxidation state (+3) in grandite garnet, the 


isotopic difference comes primarily from the mass difference 
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between *’Al and °°Fe. 

The isotope fractionation will Usually. decrease with 
increasing temperature since the difference in free energy 
between two isotopic species becomes less significant as 
the total energy of the system increases. This is illus- 
trated in Table 1 where the fractionation factor for the 


CO,-H,O system is shown for various equilibrium tempera- 


tures (6). 
Tables la. Variation, of. o with temperature (6) 
H,O-CO, 
° 
a LEG 
1.04606 0 
1.04564 2 
IIA. 2 83 Ales, 
18. 040:75 25 
MLOBSSiy . 35 
am, 0.3696 % 44.7 
5 Physical Processes 


Physical processes are equally important in producing 
natural variations in isotope abundances. Among these, 
evaporation and condensation are the most important pro- 
cesses governing isotopic compositions of precipitation. 

To understand isotope fractionation occurring in such 
processes, the meteoric water cycle may be considered. The 
most abundant isotopic components of water are Hs 10% HD? &o 


and js ereeto Their average relative abundances are approxi- 
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mately 997680:320:2000 ppm respectively (16). The natural 
abundances of these molecules vary from site to ons Seernee 
evaporation and condensation processes result in isotope 
fractionation. Fractionation factors for such processes 
depend upon the temperature and rate of reaction. If the 
process proceeds slowly, then equilibrium conditions are 
practically realized at the boundary between the liquid and 
vapor phases. The fractionation factor between the liquid 
and vapor phases then becomes simply the ratio of the vapor 


pressures of the light and heavy components: 


a = P1/Po (P,>P,) 


where subscripts 1 and 2 represent the light and heavy 
components respectively. At normal temperature, the 
fractionation factors for HDO and Fa are approximately 
1.08 and 1.009 respectively. Thus when water evaporates, 
the vapor i1s;depleted in Dp and "0. The remaining water 
which is subject to further evaporation becomes steadily 
enriched in heavy isotopes D and '°O. As the vapor condenses, 
the reverse process occurs and HDO and Hetene will tend to 
concentrate in the condensed phase. 

The Clausius-Clapeyron equation predicts that “the 
vapor pressure of isotopic species of water is a function 


of the molecular mass and the temperature of the system: 
InP = -H/RT + constant = -mLv/RT + constant [202i 


where m is the molecular mass, Lv is the latent heat of 
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Vapoetzarron ana Hets the enthalpy (85) 0) "The™errect of 
temperature on the fractionation factors for water is 
illustrated in Table 2. 

Table 2. -Fqualibrium fractionationutactors for 


hydrogen and oxygen isotopes at 
various temperatures (16) 


Tse aD BS 
100 1.029 100s 
80 ta037 1.0045 
60 1.046 1.0058 
40 1.060 1.0074 
20 iss0:79 1.0091 
Oras 1.106 ib oye 
-10 igelyes TRAE 
~20 1.146 1.0135 


The oD and a/°O values for temperatures t>20°C are those 
measured by Merlivat et al (74) and Zhavoronkev et al (110). 
The fractionation factors below O°C are calculated by 
extrapolation. 

Equilibrium conditions in this system are attained 
only inégpractically infinitely slow processes. In reality) 
non-equilibrium conditions are prevalent and they have 
important consequences for several steps in evaporation 
processes in nature. Non-equilibrium processes are often 
observed in the situation where extremely fast evaporation 


from a limited water surface occurs as in the cases of 
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closed basins and lakes. The isotopic fractionation in 
such cases becomes more complicated due to additional 
kinetic effects introduced during the phase change. 

The evaporation from open surfaces has been consi- 
dered as a molecular diffusion process into the atmosphere 
(9,73,56). The equation describing the diffusion process 


is given by 


See. q x 
ee Dem Ceeee) [23] 


where j = water flux, k = wind velocity coefficient, q is 


a function of the aerodynamical behaviour of the surface, 
Ce and Co are saturated and unsaturated molar concentra- 
tions of water in the air respectively, and D = molecular 
G22ffusivity,: eBhenuthe -kinetic»firactionation factonyis 
given by 


Cpa 
ok So fe : (2n} 


‘Gos 


Cratgq etal (14) have estimated the ratio lof diffusion 
rates from the kinetic theory of gases: 
D M M. +M 7 2 : 


1 MZ. Se Mra Me 


where Di, Dao, M, and M. represent the molecular diffusi- 
vities and masses light and heavy water and Ma is the 
molecular mass of the gas into which the water vapor 


diffuses. Their model values of isotopic fractionations 
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do not account for simultaneously observed D and O 


fractionation. Thus, Merlivat (73) proposed an alternative 


model in which he assumed a rigid elastic spherical mole- 
cule model for water. The ratio of diffusivities from 


the kinetic theory of gases is given by 


\ / 
Hee are CCE 1/2 /T, +14 2 i 
1 ae. a Gi 
where TI is the diameter of spherical molecule model. The 


equation is in relatively good agreement with observation. 
The interpretation of the equation indicates that not only 
the relative mass difference between isotopic water mole- 
cules but also the displacement of the center of gravity 
of the water molecules on isotopic substitution plays an 
important role in determining the kinetic fractionation 
effect. 

This kinetic effect, however, appears to be a unigue 
characteristic of evaporation. No kinetic effect associated 


with condensation processes has been reported (16). Expe- 


riments show that deuterium content is much less sensitive 


to kinetic effect than the }80 component since aD deviates 


. . ‘ 18 
no more: than 15% from equilibrium value, whereas a° ‘°O 


deviates 200% (16). 
6. Natural Isotopic Variation in Meteoric Water 


The study of natural isotope abundances in meteoric 


water is very important because meteoric water is involved 
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Gusatmost yall btochemical reactions occurring, on vland. The 
isotopic composition of substances formed in such reactions 
is sae tineoen by the isotopic composition of the source 
Water. Theegisotopic composition of meteoric water at 4a 
given site depends not only on temperature but also on the 
thermodynamic conditions prevailing during evaporation and 
condensation and on the initial isotopic composition of the 
source water. 

Condensation is the most important process that 
determines the isotopic composition of precipitation. The 
first small amount of water condensed from the vapor in 
equilibrium with the source water will have essentially the 
Same isotopic composition as the source water from which 
the vapor was derived. By further condensation the vapor 
becomes progressively more depleted in heavy components and 
so does the newly formed condensate from the vapor. Dans- 
gaard (16) has examined various condensation processes and 
the corresponding isotopic fractionation. If condensation 
occurs in a closed two-phase system with equilibrium 
between the liquid and vapor phases at any stage, the iso- 


topic composition of liquid and vapor can be expressed as 
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and Ee =the remaining fraction ofthe vapor phase. (15). 


In the case of isothermal condensation, 


Gees Aaa 

Teal Ue a (2.9) 
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a values of condensate and vapor formed in isothermal 
processes are shown as two practically -straight lines in 
Big. 2. | if the condensation is” caused by cooling, %iwill 
increase as the process proceeds. The isotopic composi- 
tions for such cases are shown by the dashed curves in 
PLC) dis 

Rayleigh condensation is a slow process with imme- 
diate removal of the condensate from the vapor after 
formation. Under Rayleigh conditions, condensates and 


vapor will change as 


p mt Fat i eoney 
p M1} 22 


respectively, where a, 6 and a, are fractionation factors 
at momentary temperature t, initial temperature to and 
average temperature (ttt) /2, respectively. These are 
shown as 28 and oF inwhig. 2." oIn contrast” to condensation 


under closed two-phase system conditions, fractionation 
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Figure 2. 


Isotopic fractionation of the remaining vapour 
and of the condensate as a function of the 


remaining fraction, iter 


Of Vapour loo 


ogee) = equilibrium between the total liquid 
and vapour phases during isothermal 


condensation 
6° (62). = 
VeuC 

cooling 
65(6,) = sublimation by 


composition of 
6 (S) = sublimation or 
under Rayleigh 
composition of 


same as above during condensation by 


cooling. 6s is the mean 
the solid phase 
condensation by cooling 
conditions. dc is the 
newly formed condensate. 
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under Rayleigh conditions increases drastically as conden- 
Sation proceeds. In nature, however, exchange will more 
or less smooth out the phenomenon and the actual isotopic 
compositions of water will follow curves in between the two 
extreme éf and oF curves... .iIn fact, isotopic) vapor pressure 


data collected by Craig et al (14) show that 


6D/6 (O = seat —10°C 


SD/6- 220 5 at 100°c 


corresponding to Rayleigh processes at liquid-vapor equili- 
brium. Craig (12) has obtained essentially the same 
relationships between 6D and 6'°O with water samples from 
various localities. 

Pesporacian at the site of precipitation is also 
important in determining the final composition of precipi- 
GacLOn.. .During thelr tall. trom cloud ‘to ground the rain 
drops are subjected to evaporation and exchange with the 
environmental vapor. The evaporation is relatively high 
in a dry region and probably proceeds under non-equilibrium 
conditions (21). When the humidity is high exchange 
between the environmental vapor and the newly formed con- 
densate will possibly be the dominating factor (42). 

A number of researchers have studied the fractiona- 
tion occurring during evaporation and condensation of 
water vapor as <itsis «carried in an alr mass. 


i Epstein et al (27) found that oceans, being large 
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reservoirs of water, have a faivly tunLiomm Ls0bop ie 
composition disregarding those parts which are | 
directly mixed with fresh water. 

2 The deuterium and oxygen 18 concentrations of 
méteorie water vusually*vary#in parallel (4027 22016) - 
In general, deuterium and oxygen 18 variations in 
water from various rivers, lakes, rain and snow 


follow a relationship 
6D = 86180 + 10 (Craig a2. [33] 


except for water from closed basins in which evapora- 
tion is the dominant factor governing the isotopic 
relationship. 

3% The heavy isotope content in precipitation decreases 
as the condensation temperature drops. This is 
related to several factors such as latitude and 
aleLeudes (16). 

4. There is also a Paneonat effect on the isotopic 
composition of rain water; isotopic composition is 
lighter in winter precipitation than in summer 
precipitation because of cooler condensation tem- 


peratures in winter (54). 


Tae Isotope Thermometry 


The isotopic composition of molecules formed in 
natural processes depends not only upon the temperature 


of the system but also on the overall isotopic composition 
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of fluids and other solid phases in equilibrium. When two 
solid phases formed in equilibrium are considered the iso- 
topic fractionation between the two will be determined only 
by temperature and it will be independent of the presence 
or quantity of any other substances. This suggests that 
for a pair of co-existing molecules, if the relationship 
between the fractionation factor and the temperature is 
known (either by experiment or theoretical calculation), 
the measured fractionation between the pair will yield the 
temperature of last equilibrium. This is the basic 
principle behind the so-called "thermodynamic thermometry" 
of stable isotopes (65). Most of the geothermometers 
dealing with the formation temperature of minerals belong 
im chis category.* for example, Bottinga-et al) (7), found 
the fractionation factors between quartz and magnetite in 


the temperature range 600°C to 850°C to be related by 


6 


5.57 x 10°r? 


[34] 
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where Q = quartz, M = magnetite and T = absolute temperature. 
In practice, instead of measuring absolute isotope 

abundance ratios, it is customary to measure the isotope 

ratio relative to a standard and to quote the result as a 


6 value defined by 


1} x 1000 [35] 
{ 


where R_ = isotope ratio of the unknown sample and Re 
>, 
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isotope ratio of the standard sample. The most commonly 
used standard is SMOW (STandard Mean Ocean Water) whose 
@-"O "and 6D values are ee yeciaa to zero. A positive value 
for 6 means the sample is enriched in the heavy isotope 
relative to sea water and a negative § value means the 
sample is depleted in the heavy isotope relative to sea 
water. 

The fractionation factor between two substarices A 
and B can be expressed in terms of 6 values: 


R {=6 sg, nals \ 
a epee sm ie 1/380 a bl 


AB R , 1000 1000 


Taking logarithms on both sides, 


If 6, and 6,<<1000, 


* - : econ 
1000 Indap Sy 6. 
Since most biological reactions occur at near normal 
environmental temperatures, the fractionation factors 
involved in such processes are expected to be an inverse 
function of the absolute temperature: 


a 
~ = me ES CE 5} (372) 
1000 Inada, on é, - 


where a and b are constants. In addition, thes cotcaws 
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variation of temperature in such reactions is very small 


(a few degrees). Thus the equation can be simplified to 


1000 inoa, a on = 5. Sio-Cle cs sb [38] 


. where C is a positive constant. Thus lna turns out to be 
a negative linear function of temperature. The equation 
indicates that by measuring the difference in isotopic 
composition of the two solid phases formed in equilibrium 
One can deduce the temperature of the equilibrium. 
PLrequehely, 1c. 15 Gitticult to find, cwo. co-ex1sting 
solid phases in equilibrium which show a significant 
temperature-dependent fractionation between them. For 
example, equilibrium assemblages secreted by organisms, 
such as calcite, silica, and phosphorite show distinct 
temperature dependencies of oxygen-isotope fractionation 
between water and the minerals (13,61,70). Unfortunately, | 
their temperature coefficients are practically identical 
resulting in no oxygen eereene fractionation between the 
co-existing minerals. In such cases, absolute temperature 
information cannot be obtained and a variety of environ- 
mental factors controlling the isotopic composition of 
the individual substances within the system must be taken 
into consideration to obtain temperature information. 
This is the so-called “environmental thermometry", intro- 
duced by Lerman (65). Most of the isotope thermometers 
obtained from biological systems belong in this category. 


To determine past climate it is necessary to find 
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natural systems which have recorded the isotope variations 

in sequence. Such systems must satisfy a number of 

eniternia, (48)... 

k. The isotope variations found within the system must 
be -determined solely by climatic .factons onlin «such 
a Way 7that,.non-climatic,influences,canwbe .filtered 
Out. 

As The isotope record must be permanent. Thus it must 
be proven that no further fractionation processes 
Such as re-equilibrium or exchange have occurred 
after the strata have been deposited. 

cy The record preserved by the strata should be 
reasonably complete and continuous. 

4, The strata must be datable so that a time scale can 
be attached to the record. 

awe The time resolution of the stratified system must 


be compatible with the climate information required. 


Natural systems When, to a greater or lesser extent, 
satisfy cas above criteria include: deep ocean and lake 
sediments, polar ice sheets, speleothem deposits, peat 
deposits, annual coral rings and annual growth rings of 


trees. 


S. Tree Rings 


Living organisms consist mainly of the elements 
oxygen, carbon and hydrogen which are incorporated into 


their tissues from the environment by metabolic processes. 
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In consequence, the abundance of stable isotopes of these 
elements in plant and animal tissues depends both on the 
isotopic composition of the environment and on the isotope 
Eract Ponationpoccurring! during photosynthesis» and other 
metabolic processes. The use of biological systems as 
climate indicators presently presents complex isotopic- 
geochemical problems. However, high rates of deposition 
and good time resolution suggest great potential for these 
methods. Furthermore, biological systems present intrigu- 
ing possibilities in that they contain three potential 
isotope indicators which, in any given system, may be 
independent of one another. These are the isotope ratios 
189/169, D/H and !3c/!2¢c,all of which are measurable in 
a Single sample. Thus multiple over-determination of 
temperature using three isotope indicators may provide an 
excellent means of cross-checking the result with some 
level of confidence. It is also possible that non- 
temperature effects will be demonstrated by the use of 
multiple thermometers. 

One type of sample which appears to be suitable as 
a recorder of past climate and environmental conditions is 
matter of vegetable origin, e.g. tree rings, peat, humus, 
etc. Other samples of organic origin which might also 
preserve isotope records are those of animals e.g. teeth, 
antler, bone, skin, eggshell, etc. Animals reflect basi- 
cally the isotopic compositions of vegetables on which 


they feed (18,36). 
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The*systems studied to-dateareVtree rings “and to a 
lesser extent peat deposits. Tree rings have great poten- 
tial as climate indicators for a number of reasons: 

Ee Annual growth rings of trees provide a time resolu- 
tion of at least one year. 

at Absolute dating by means of dendrochronology is 
accurate to one year. 

Sis A long continuous time period can be investigated. 
The bristlecone pine chronology from the White 
Mountains in California, for example, has a length 
of 7000 years (33). 

4, Tree samples are widely distributed throughout 


almost all continental area. 


SL D/H Ratios in Tree Rings 


The use of the isotopic composition of organically 
bound hydrogen as a paleothermometer has been recently 
suggested by. Schiegl ._(90,91)..and Libby (66,67). .Schiegl 
et al (92) measured D/H ratios of both marine and land 
plants and found systematic variations of deuterium content 
in plants. For example, the hydrogen of both marine and 
land plants is significantly depleted in deuterium compared 
to, the water in .whichsthe ,phants ,qrew..They,.also, found, a 
systematic relationship between the natural deuterium 
content of organically bound hydrogen and that of precipi- 
tation. Since the deuterium content of precipitation 


depends mainly on the climate, particularly on mean annual 
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mpemperature, (16) ,Vochnegl! (92) concluded that plant 
Material should Garry information on past climate. In’ an 
attempt to check the conclusion, Schiegl (90) measured 

the deuterium content of peat from The Netherlands and 
found a relatively good relationship with mean July tem- 
peratures. The method was subsequently applied to growth 
rings of picea from southern Germany (see Fig. 3). Schiegl 
(91) found that deuterium variations in tree rings of picea 
are essentially a function of annual air temperature. 

Short isotopic variations are partly influenced by changes 
of relative humidity. On a long term basis, however, it 

is the mean annual temperature variations which are res- 
ponsible for deuterium variations. 

The isotope thermometer introduced by Schiegl is an 
example of an environmental thermometer which cannot 
provide absolute temperatures but yields empirical tempe- 
ratures which may have been altered by various environmental 
factors, sucheas 6D values of water, and relative humidity. 
Nevertheless, this was the first isotope study that 
examined the possible use of stable isotope data stored 
in tree rings as a paleothermometer. 

Libby (66) ).0n the other hand, utrlized a thermo- 
dynamic approach. Libby calculated temperature coeffi- 
cients of several isotope thermometers in bio-organic 
material using statistical mechanics. The major assump- 
tions made in this calculation are as follows: 


ibe Cellulose is formed in equilibrium with water. 
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Figure 3. Comparison between summer temperature (May- 
September), annual mean temperature, deuterium 
content of Picea, relative humidity (May- 
September) and the calculated deuterium 
content of the leaf sap all in the 5-year 
means. Dashed lines indicate the trends in 
the histograms (91). 
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Pe Only the internal vibrational frequencies and the 
bending frequencies associated with rotational 
energies are included in the calculation of 
Darnkut tomslunctions. 

ce The internal vibrational frequencies used were 
assumed to exist in cellulose molecules. The 
actual presence of such frequencies and any other 
possible frequencies were not verified experimentally. 

The fractionation factors and their temperature coeffi- 

cients thus calculated are shown in Table 3. 

Tables. ~ Fractionation factors computed tor 


isotopes of hydrogen and oxygen in 
cellulose (66) 


Fractionation Temperature Temperature 
a aK coefficient 
273 298 fool °C 
a}®0,-cellulose on aleyse 0.96 
oko -cellulose 185200 Irs 4 Teas 
CO, 
8g) Veeinitose Tap ey ds 1056 ogo 
H0 
*PHO-cel lulose (OH) 1.5881 7.50163 0.4 
Py O-cellulose (CH) 1.8095 18606 2.0 


For example, the temperature coefficient of the 
fractionation factor between water and carbon-bound 
hydrogen in cellulose was found to be 2°)/ ae GEDDY, 


has compared this value with the measured temperature 
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eoetiicient or -3°/..,/°C obtained by Schieg] (85) mtonra 
peat deposit in The Netherlands. However, this comparison 
is not valid because the value 2°/,,/°C represents the 
Eemperature coefficient of hydrogen isotopic fractionation 
between water and CH groups in cellulose whereas the value 
3°/,./°C is the temperature coefficient of deuterium 
content in whole wood. 

For the purpose of comparison, the properly calcu- 
lated temperature coefficient of 6D values of cellulose 
would sbe 1. 1° 7.5/7 °C provided that. the 6D of water is 
constant. This temperature coefficient appears to be too 
small to be used as a thermometer because the average 
precision of §D measurement is about +2°/,,. Furthermore, 
the magnitude of the variation in annual mean air tempera- 
ture is generally small. For example, the total variation 
in mean annual temperature in southern Germany for the 
period Schiegl examined is less than 2°C. The corres- 
bonding, 6D. variation would, then be 2.2.97, which is 
within the random experimental error. Therefore, the 
small temperature coefficient calculated by Libby may mean 
that cellulose—-water in tree rings cannot be used as a 
thermodynamic thermometer. However, D/H ratios preserved 
in tree rings may provide an environmental thermometer. 

It is apparent that more work should be done to examine 
this temperature-dependent fractionation of hydrogen 


isotopes in tree rings. 


To compensate for the limitation of tree rings as 
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a thermodynamic thermometer, Libby (66) has suggested the 
use of multiple thermometry in tree rings. In photosyn- 
thetic production of bio-organic materials, all three 
elements C, H and O may yield independent thermometers. 

In contrast to the relationship between D and !°0O in 
meteoric water, the three elements are not functionally 
related in photosynthesis. The multiple overdetermination 
of temperature using a set of thermometers may be useful 
for cross checking the isotope result. In particular, it 
may be possible to detect some effects that are not caused 
by temperature and to filter out non-climatic effect. 

Libby et al (67) measured the isotope ratios of C, 

H and O in tree ring sequence of a German Oak and corre- 
lated them with the existing weather records from England, 
Basel and Geneva to evaluate the empirical temperature 
coefficients. The temperature coefficients are listed in 
Table 4. 

The isotope data 2cbeltts with winter temperatures 
recorded in England are plotted in Fig. 4. Unfortunately, 
it is impossible to compare these results with others 
because the isotope ratios were measured relative to an 
internal working standard instead of an internationally 
accepted standard. The magnitude of the temperature co- 
efficients obtained is large. This was attributed to the 
fact that the temperature coefficient is a combination of 
equilibrium fractionation during photosynthesis, and both 


equilibrium and non-equilibrium fractionation during 
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precipitation. Epstein et al (30), however, demonstrated 
that the combined effect of thermodynamic and environmental 
temperatures could not account for such large temperature 
coefficient. 

Table 4. Temperature coefficients computed from 


Measured isotope ratios and weather 
records (67) 


Isotope ratio Temperature coefficient No.2 OL 
Samples 
EGY ase) (5.2920.68)T yang he GN: 68 
(2.9140.41)T, 2 561 9 emus sak 56 
(2.8640.52)T LO ova ily ay ON 56 
D/H (89.5416.-0)Tegiand Pe eC 68 
(SC UE Aa ay Ay AL 56 
(Ml Ae Bc DL es es nig W/otci/ a 56 


Nevertheless, subsequent work done by Libby et al 
(69) demonstrated qualitative correlation between air 
temperature and stable isotope ratios in tree rings... Their 
subjective choice of temperature data for correlation pur- 
poses was again criticized by Epstein (30). 

The isotope results of Schiegl and Libby were obtained 
using whole plant material and suffered from several comp- 
lications. Plants are extremely complex chemical systems. 
Wood for instance consists of approximately 40-50% cellulose 


and 20% hemicellulose, the rest being mainly lignin (30%) 
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Figure 4. 


Winter temperature 
in England (-C) 
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Deuterium and oxygen ratios in German Oak 
compared with English winter temperature 
(30-year running averages) (67). 
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with small amounts of other material (see Table 5s 


Table 5. Average % chemical composition of 
softwoods and hardwoods (95) 


Component Softwoods Hardwoods 
Cellulose 42 45 
Hemicellulose 27 30 
Lignin 28 20 
Extractives 3 5 


Isotopic differences between chemical substances in 
plants may affect the D/H ratio determined for the total 
System. Smith and) Epstein (93) found that the: lipid 
fraction of plant material is approximately 100% depleted 
in deuterium with respect to the plant total hydrogen. The 
lignin fraction of whole wood is some 10°/,, lighter in 
¥®O than the cellulose fraction (51). Furthermore the 
percentage composition of chemical substances is not 
climate dependent (51). Tor reso ive the problem associated 
with chemical heterogeneity of plants, Epstein et al (31) 
and Wilson et al (107) used a single substance, namely 
cellulose, instead of whole wood. 

Cellulose is the only substance in plants whose 
chemical structure is relatively well defined. The defini- 
tion of cellulose has been a point of discussion among 
workers in the field. There is a general agreement that 
cellulose is composed predominantly of §-linked D-gluco- 


pyranose units (53). Thus, extraction of cellulose in the 
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pure form from woody tissue is physically impossible. The 
term cellulose in this study, therefore, is referred to as 
a=cellulose which is alkali-resistant cellulose retaining 
O7orisnosiOr non=—gilucan’ carbohydrates. ‘The molecular 


Structure of cellulose can be represented by 


CH.OH 


In a cellulose molecule, the hydrogen atoms fall into 
2 groups. One is carbon-bound hydrogen (CH) and the other 
is oxygen-bound hyceeeen) (OH) comprising 70% and 30% res- 
pectively. The bond strength of a CH bond is greater than 
that of an OH bond. Thus their hydrogen isotopic composi- 
tions may reflect two different fractionations. An 
additional complication is the problem of isotopic exchange 
of hydrogen. Carbon-bound hydrogen in cellulose is non- 
exchangeable under normal conditions whereas the hydroxyl 
hydrogens exchange readily (43,72). In addition; tine 
structures of cellulose also influence the hydrogen iso- 
topic composition of cellulose. Cellulose consists of 


highly ordered crystalline as well as disordered amorphous 
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tegions (59). "With different celluloses, the relative 
amount of material in the crystalline and in the amorphous 
regions, as well as the degree of order within the crystal- 
line region, are not the same. Infrared spectral analyses 
By tannin et val 1(72) show that about 25% of the hydroxyl 
groups in cellulose are in the crystalline region and that 
about 30% of this hydroxyl hydrogen together with all the 
hydroxyl hydrogen in amorphous regions exchange with water. 
In other words, 25% of the total plant hydrogen in cellu- 
lose comprising 22% from amorphous OH group and 3% from 
crystalline OH group exchanges with the hydrogen in water. 
To resolve the problem of exchangeable hydrogen in 
cellulose, Wilson et al (10) equilibrated the hydroxyl 
hydrogen with water of known isotopic composition. Using 
this technique, they analyzed D/H ratios in aye cellulose 
from a Monterey Pine (Pinus Radiata) growing in Hamilton, 
New Zealand (see Fig. 5). They found that tree rings 
formed in summer were more depleted in deuterium than 
winter rings. Since the isotopic composition of the cellu- 
lose varies in the opposite direction to the isotopic 
composition of precipitation, they concluded that the D/H 
ratio of cellulose in the wood of Monterey Pine changes 
with temperature due to the temperature effects on bio- 
chemical reactions leading to cellulose. Taking into 
account the isotopic variations of the atmospheric pre- 


cipitation, they obtained a temperature coefficient os ig GA 


per °C. 
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Figure 5. 
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This work has been criticized in that the negative 
temperature coefficient is opposite to that expected from 
thermodynamic arguments. According to thermodynamics, the 
isotope fractionation in a chemical reaction normally 
decreases as the temperature of the system increases. To 
Satisfy the thermodynamic condition, the D/H ratio of 
Belliiose should increase as the temperature of the system 
increases because the hydrogen isotopic composition of 
meteoric water is a positive linear function of air tem- 
perature. There are several possible explanations for this 
contradiction. The phenomenon may be explained in terms 
of a time lag between the actual photosynthesis and cellu- 
lose deposition. A conifer manufactures photosynthates at 
all times of the year except when climatic conditions are 
unsuitable. These photosynthates are stored for short 
Or extended periods before they are laid down as wood. 
Thus the cellulose in late wood of Monterey Pine might have 
been produced in summer instead of winter and thereby show 
aonighs D/H ratio retlecting the D/H. ratio, of summer pre— 
eipitation,....Wwhen, this time. lag is considered, Che rsysten 
may show a positive temperature coefficient. Another 
possible explanation was given by one of the authors. 
Wilson et al (107) assumed that the fractionation between 
the hydrogen of the water and the OH hydrogen of the 
cellulose was small. Wilson (106), however, pointed out 
that the fractionation may be very large, perhaps due to 


the extensive hydrogen bonding in the cellulose. The 
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fractionation must be of the order of 200°/,, such that 
the OH hydrogen is enriched in deuterium. The correspond- 
poo Cemperacure=coctriebent 4s, in the Srder@or™ 10eveeysc. 
Thus if the crystalline OH hydrogen atoms have a very 
large negative temperature coefficient, outweighing the 
positive temperature coefficient of the CH hydrogen, the 
net result would.be a negative temperature coefficient. 
However, the temperature coefficient found by Libby for 
the fractionation between OH hydrogen in cellulose and 
hydrogen in water had not only a very small magnitude but 
also opposite sign. 

On the other hand, Wilson (106) proposed a possible 
thermodynamic thermometer in tree rings, which utilizes 
both CH and OH hydrogen of cellulose. He hypothesizes 
that the OH hydrogen in the crystalline region of cellulose 
may provide a sensitive paleothermometer which would be 
independent of the cell water composition. The isotope 
ratio of OH hydrogen could bé determined by measuring the 
6D of non-exchangeable CH hydrogen by replacing OH hydrogen 
with “NO, groups and non-exchangeable hydrogen (both CH and 
OH) determined by the equilibration method. Using this 
technique, Friedman et al (41) calculated 6D values of OH 
hydrogen in cellulose and correlated it with the altitude 
of growth sites to find the temperature relationship. The 
non-exchangeable OH hydrogens thus obtained were more 
enriched in deuterium than sap water by as much as 650°/,, 


indicating that the fractionation between OH hydrogen in 
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cellulose and hydrogen in water is*Targeme in -contrasti"'to 
Wilson's prediction however, temperature dependence of the 
fractionation was not detected, which may be in keeping 
with Libby's prediction. It is of interest to note from 
the result of Friedman et al (41) that non-exchangeable CH 
hydrogen in cellulose is more enriched in deuterium than 
that Ofistrean water byabout .30°/,.>).. .n.wone,. case 6D 
values of CH hydrogen were found to be even higher than 
that of sap water. These are opposite to the results 
obtained by Epstein et al (31) and White et al (103) in 
which CH hydrogen, regardless of temperature and relative 
humidity, is more depleted in deuterium than environmental 
water and sap water by about 22°/,,. However, it is 
apparent that further work with a well defined technique 
is required to resolve this contradiction. 

A second method of dealing with the exchangeable 
hydrogen in cellulose has been reported by Epstein (31). 
Plant samples were nitrated and pure cellulose nitrate 
was extracted by acetone dissolution. Subsequent tests of 
the nitrated product have demonstrated that the nitration 
extraction procedure eliminates the OH hydrogen in the 
cellulose but does not alter the D/H ratio of carbon-bound 
hydrogen. 

Using this technique, they measured the D/H ratio of 
non-exchangeable hydrogen in cellulose extracted from 
aquatic and land plants. They found that plant-extracted 


cellulose nitrate D/H ratios are systematically related to 
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the D/H ratios of the associated environmental water (see 


Fig. 6). The equation thus obtained is 


J 7 SDE 


igi 


where cell and w represent cellulose and environmental 
water respectively. 

From the fact that the overall relationship is 
linear with slope close to unity, they concluded that D/H 
ratios of cellulose in tree rings reflect mainly the D/H 
ratios of environmental water. In support of this, White 
6€ alg(i03) found a similar linear relationship between 
the D/H ratios of the non-exchangeable hydrogen in tree 
cellulose and the model-calculated D/H ratio of tree sap 


for eastern white pines: 


=" OD - 21 [39] 


soot sap 


where sap represents sap water. 


The slope and the intercept of this equation are the 


same as those of the equation obtained by Epstein et al 
(31) despite the fact that White used sap water whereas 
Epstein used environmental water. This may indicate that 
absorption of water by plant roots involves no fractiona- 
tion. The temperature distribution over much of the 
earth's surface has been correlated with the pattern of 
D/H and !80/!*o0 variations in meteoric water (16). Thus 


natural systems which record the variations with time of 
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Figure 6. 
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The relationship between the 6D of the 
cellulose nitrate extracted from plant 
material and the 6D values of meteoric 
waters sassoczated ‘withthe ;plants 7631). 
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the D/H or '°0/'°oO ratios of meteoric waters at a given 
site should preserve a record of past temperature variations 
at the site. 

Epstein et al (31) also found non-systematic but 
Significant differences between the values of plant total 
hydrogen and non-exchangeable carbon-bound hydrogen (see 
Fig. 7). Unlike non-exchangeable carbon-bound hydrogen, 
plant total hydrogen has a poor functional relationship 
with environmental water. They attributed this poor corre- 
lation to the effect of chemical heterogeneity on 6D value 
of plant material. 

Epstein et al (29) subsequently applied this technique 
to paleoclimate reconstruction. Measuring D/H ratios of 
non-exchangeable hydrogen in cellulose, they obtained avdD 
record spanning 1000 years from two bristlecone pines, 
whose growth periods overlapped by about 100 years, from 
the White Mountains in California. This record is compared 
with Lamb's winter temperature curve for central England 
and ring width variations “ina “bristlecone™pine’ from the 
White Mountains in California (see Fig: 8). There is 
qualitative agreement between records. For example, all 
four records in.Fig. 8-show a cooling period between ,1450 
and 1800 AD which. is often. called "The Little.ice Age”. 
Epstein et al (29) also obtained 6D records from a Scots 
pine located near Loch Affric, Scotland (see Fig. 9). The 
average 6D value of the scots pine is about -60°/,. as 


compared to -95°/,, in the bristlecone pine for the corres- 
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B0cuow OF RAW WOOD VS. CELLULOSE MITRATE 
OF BRISTLECONE PINE, 


WHITE MINS., CALIF. 


TiGOe Pasoses1GO™™ =13092 1209 UNO EAENDG OE=g0 
SDemow (Yeo) OF CELLULOSE NITRATE 


The relationship between the 6D values of raw 
wood from a single bristlecone pine and the 
corresponding cellulose nitrate extracted from 
the raw wood. Each point represents a diffe- 
rent 10-year interval of growth (31). 
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ponding time interval. The difference is in keeping with 
the difference of the SD values of meteoric waters in their 
respective localities, -42°/,, for Loch Affric and about 
1007/4. ©Or ehefwaters in the bristlecone pine area. & 
comparison of the 40 year running average of the 6D values 
of the two pines produces a relatively good linear relation- 
Ship which leads Epstein to suggest that long term isotopic 
climatic trends, as recorded by the two species of pine are 
Similar, while the short term fluctuations are not related. 
Epstein et al (29) conclude that widely different types of 
climatic environments produce 6D trends which may reflect 
large scale long-term climatic changes, superposed upon 
which there may exist short term 6D variations which could 
reflect temperature changes over a restricted local area. 
The ultimate purpose of paleoclimate reconstruction 
is to identify possible periodicities in climate variations 
and the causes of such variations. Libby et al (68) cal- 
culated climate periods from Stable isotope ratios of 
oxygen and hydrogen measured in a Japanese cedar spanning 
1800 years. These climate periods were compared with those 
from 800 year Greenland ice (17) and with periods computed 
from the tidal stresses of the Sun-Moon-Earth system. 
According to Libby a good correlation was evident. Also, 
Fourier analysis, done by Epstein et al (29) on the isotope 
data of the bristlecone pines from the White Mountains, 
showed a 22 year periodicity in the power spectrum. They 


Telated it to the 22 year drought cycle in the midwestern 
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Bristlecone Pine, White Mts, Colifornia- Cellulose Nitrate 
10 year intervais 


(a) 


-90 Bristiecone Pine, White Mts., Coliforma - Cellulose Nurcte, 40 veor 
running average of decadal means 


45 Winter Temperatures (D-J-F) of Central Engiand 


Ring Widths of Bristlecone Pine, White Mts., California 


Ring Width (mm) 
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Figure 8. 


a) The 6D record of isotopically non-exchangeable 
hydrogen of cellulose (cellulose nitrate) in two 
bristlecone pines which grew in the White Mountains 
Of «Ca luifornia-. 


b) The 40 year running average of the 6D record. 


c) Winter temperatures of Central England (Dec-Feb) 
from a variety of indicators. 


d) Absolute ring width variations in a bristlecone 
pine from White Mountains, CalLiorniay (63529). 
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Figure 9. 


a The “oD record of cellulose nitrate from 5-year 
segments of a Scots pine growing in Scotland. 


b) The 40-year running average of decadal means for 
winter temperatures (October through March) of 
Edinburgh, Scotland. 


c) The 40-year running average of decadal means of 
6D Gt the Scots pine (29).. 
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U.S. (84). The 22 year cycle of drought in turn has been 
attributed to the Hale sunspot cycle which seems to occur 
at about 22 year intervals (76). The validity of the 22 
year cycle in the power spectrum of the bristlecone pine 
isotope data, however, is not firmly established due to 
possible aliasing problems. In any time series analysis 
of discrete experimental data, the maximum statistically 
meaningful frequency is the Nyquist frequency defined by 
Ve a [40] 
where T is the sampling interval (60). The sampling 
interval chosen by Epstein et al (29) was 10 years. The 
corresponding Nyquist frequency becomes 1/20 which is very 
Close to the frequency of the 22 year cycle. Thus, the 
22 year cycle may be an artifact of the Fourier analysis 
of data contaminated with colored noise. Sampling interval 
of 5 years or less may resolve the presence of the 22 year 


eyele periodicity if the cycle is indeed real. 


See2 18071°O Ratios in. Tree Rings 


Oxygen isotope ratios have been used extensively in 
paleoclimate studies for a number of reasons. Firstly, 
there exists a large reservoir of oxygen in the oceans. 
Secondly, oxygen is very reactive and forms compounds with 
most other elements and lastly, the natural abundance 
variations are relatively easy to measure. For example, 


ives natural abundances of the three atmospheric oxygen 
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meotopes;| = ©, O, and *°O are 99..759%., 020372, -and 
0.2039% respectively (78). 

Oxygen isotope ratio measurements on cellulose 
extracted from trees have been reported by Gray et al (49, 
50), and-Epstein et al (28). Gray et al (49) have reported 
‘89 values of cellulose from successive 5 year groups of 
rings in a spruce tree (Piecea glauca) which grew in 
Edmonton, Alberta (see Fig. 10). A striking linear rela- 
tionship between 6!°O and mean annual temperature was 
obtained (see Fig. 11). The temperature coefficient thus 
obtained is 1.3 0.1°/,,/°C with a correlation coefficient 
of 0.97. The correlation between §'°O and mean winter 
temperature iS good whereas it is poor for mean summer 
temperatures. The best correlation occurs with September 
to August mean annual temperatures. Since spring and 
summer are the time most conifers lay down growth rings, 
the poor correlation between 6!°O and summer temperature 
requires some explanation. ane et al (49) proposed two 
possible explanations for this observation. First, in 
addition to any temperature effect on the fractionation 
factor during photosynthesis, there is also variation in 
the isotopic composition of the source water being used to 
produce cellulose. They assumed that the cellulose oxygen 


is derived from CO. which undergoes isotopic exchange with 


2 
HO before photosynthesis occurs, and thereby reflects 
189/155 ratios of the source water. Thus it is possible. 


that winter precipitation stored as soil water after melting 
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may play a significant role in the production of cellulose 
during the growing season. Second, it has been shown that 
net photosynthesis and growth in coniferous trees depend 
markedly on the production of stored foods. Thus winter 
temperatures may have a direct bearing on the isotopic 
composition of the cellulose produced the following year. 
This linear relationship between §1°O of cellulose 
and mean annual temperature has been criticized by Epstein 
etal (29) on the basis that 6/°O values of cellulose are 
affected by evapotranspiration and humidity as well as the 
© °O of the meteoric water. In support of this), Ferhi @t 
al (35), in experiments on bean plants grown under con- 
trolled conditions’, found that the variation of 6/°O value 
of cellulose with temperature was very small if not negli- 
gible. In a subsequent series of experiments, Ferhi et al 
(34) concluded that the most important factors in deter- 
maining. O77 °O.ratzo-o4 plant cellulose were the isotopic 
composition of water absorbed in the roots and relative 
humidrty. To clear this apparent. contradiction; Gray eteal 
(52) measured 6*°O of cellulose for\a number of [trees fvom 
various sites in North America. The overall mean §'°O and 
annual temperature data for each site are DLOLeEea in fig. 
12. The least squares estimate of the equation defining 


these data is: 


S28 Oo c= OL o47 + 22.9. [41] 


These data are compared with 6'°O data obtained from 
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sphagnum moss. The slopes of the two lines are the same 
within experimental error while the intercepts are diffe- 
rent. The fact that 6'°0O values in cellulose of tree 
rings, and sphagnum moss have a linear relationship with 
air temperature suggests that in many of the regions 
temperature information is stored in the 6!°0O of cellulose. 
Furthermore, the virtually identical temperature coeffi- 
Clients shown by two different biological indicators suggest 
that enrichment by evapotranspiration must be relatively 
constant. Evapotranspiration is known to be controlled 
mainly by water availability and relative humidity, and is 
only indirectly related to temperature. The sphagnum moss 
samples were taken from peat bogs where humidity is 
generally high and evapotranspiration will be minimized. 
Thus the displacement of the tree ring line may be explained 
in terms of enrichment by evapotranspiration. This is 
possible only if evapotranspiration is assumed to be 
constant over a wide range of mean annual temperatures. 

ini Support OL this, a similar functional relatzonsnip 
between mean annual air temperature and §'°O values of 
cellulose in trees ge aT coastal stations in North 
America has been reported by Burk et al (10). They found 
a linear relationship between 6*°O of cellulose from tree 
rings and the latitude of the collection site. When con- 
version from latitude to mean annual temperature was made, 
a linear relationship between mean annual temperature and 


189 of cellulose was established. 
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EXPERIMENTAL EQUIPMENT AND PROCEDURES 


ie Gas Extraction Methods 


To determine D/H and '°%0/!°O ratios of organic 
material, it is necessary to extract the desired elements 
in a suitable form for mass spectrometric analysis. 

Oxygen is normally converted to CO, while for hydrogen 
isotopic analysis, hydrogen gas is the most suitable form. 
The extraction of these gases from organic materials 
involves various chemical and physical processes which may 
introduce undesirable secondary fractionation and alter 
their-original isotopic records. 

A number of methods for the separate extraction of 
hydrogen and oxygen from biological material have been 
described. For the extraction of oxygen, there are a 
variety of methods (83,96,55). For hydrogen, only one 
method is generally used. This involves the quantitative 
conversion of hydrogen in a sample to water followed by 
reduction of the water over hot uranium metal to hydrogen. 
These methods are time consuming and often impossible for 
small samples. In some cases it is desirable to measure 
bothuD/H.and..'°07)2O ratios,on axygiven sample. Bor this 


reason, an attempt has been made to analyze simultaneously 


hydrogen and oxygen from a single sample. 
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1.1. Pyrolysis of Cellulose in a Quartz Tube 
Se ee ees in a Quartz tube 


fils. Procedure 


The first method used involves the pyrolysis of 
cellulose in a quartz tube followed by separation of the 
desired gases from the pyrolysis products. About 30 mag 
@f cellulose isplaced in the bottom of a-quartz tube 
ie emiin diameter and ll cm*long. © The quartz tube 
containing the cellulose is then evacuated and sealed with 
a heating torch. The cellulose in the quartz tube is then 
py Lolyzed in’a furnace,at 1200°C for.20 minutes, @arter 
which the quartz tube is mounted on the separation line 
(see Fig. 13). The gases from the pyrolysis are released 
by breaking the quartz tube mechanically using a drop 
weight and the individual gases which include Ho, HD, CO. 
and CO are separated by means of chromatographic adsorp- 
feu ON. 

The gas separation Wine, consiststof a spiral ‘trap, 

a silica gel trap and a toepler pump connected in series. 
As the gas mixture from the pyrolysis passes the 
chromatographic columns, CO. is held in the spiral trap 
kept at liquid nitrogen temperature and CO is held in the 
Silica gel trap kept at liquid nitrogen. Hydrogen gas is 
collected with the toepler pump. The yields of gases thus 
separated are measured manometrically. 

An alternative method for the separation of 


hydrogen gas from the mixture of gases was used. This 
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involves diffusion of hydrogen through a 75% palladium - 
25% silver alloy wall maintained at 500°C. Hydrogen 
diffuses through the thin wall of the palladium-silver 
mube and 16 collected, on charcoal at liquid nitrogen 
temperature. Isotopic ratios of hydrogen and oxygen gases 


were then determined by mass spectrometry. 


ote. . ReESULt Vand. Discussion 


To determine the optimum experimental conditions for 
Dine complete extraction of hydrogen and oxygen from 
cellulose, samples of Whatman filter paper cellulose were 
pyrolyzed in quartz tubes under various conditions. The 
results are shown in Table 6. 

The yield of H, gas was about 90% of the theoretical 
value indicating incomplete decomposition of cellulose. 
To enhance the yield various combinations of reaction time 
and temperature were employed. At a temperature range of 
1280-1320°C, increasing the pyrolysis time produced 
essentially no increase in the yield of H, gas. The vari- 
ation of temperature in the range of 1280-1320°C appears 
to have no effect on the yield of hydrogen. A temperature 
higher than 1300°C is physically difficult to obtain due 


to fast oxidation of heating elements in the furnace. 
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Table Gye s°Yields*of Hydrogen from Pyrolysis of "Cel lulose 


fiteaket- Y1e1d "lemp, “uration D/H ratio Cold Trap 
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The low yield of hydrogen may be due to the produc- 
tion of methane gas during pyrolysis. Mass spectrometric 
analyses of the pyrolysis products showed the presence of 
methane including deuterated compounds, such as CH,D., 
CHD, CHD, and CD,. Approximately 10-20% of the total 


3 3 4 
hydrogen in cellulose was in the form of methane at the 
end of pyrolysis. A similar result has been reported by 
Holliday et al (58). For thermal decomposition of methane 
at temperatures up to 1300°C, they found that about 90% of 
the methane converted to Ho and the rest remained unaltered. 
Thermal decomposition of methane at this temperature is 
believed to take place mainly on surfaces. Thus the 
Bpeaction 1S;rapid initially and the reaction is virtually 
complete after 10 minutes at which time the reaction is 
strongly retarded in the presence of hydrogen gas. The 
retardation is due to the preferential absorption of 
hydrogen which protects the surface so efficiently that 
eventually further decomposition of methane is prevented. 
Incomplete conversion of cellulose hydrogen to free 
hydrogen gas inevitably leads to isotopic fractionation 
between the hydrogen gas molecules and the compounds 
containing hydrogen. For example, methane exchanges 
hydrogen isotopes with hydrogen gas molecules at virtually 
all temperatures. Therefore the D/H ratio of the hydrogen 
from the pyrolysis is different from that of the original 
hydrogen in the cellulose. 


The apparently high yields of hydrogen fore tital 2o- 
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41 in Table 6 are likely caused by the presence of 
impurities, mainly methane gas, in the hydrogen gas phase. 
Depending on the temperature of the silica gel trap, some 
methane gas may escape the cold trap and mix with 
hydrogen gas increasing the apparent hydrogen yield. The 
methane gas, when ionized inside the mass spectrometer 
geses part of hydrogen gas to form radicals such as CH., 
CH, etc. The isotopic composition of the hydrogen thus 
released from methane is different from that of hydrogen 
Mase yelhis 15 allustrated in Table 7; 

Table 7. The effect of methane content in hydrogen 


phase on the overall isotopic composition 
of hydrogen 


fora. H. Yield CH,/H, D/H RacLo 

% % a /aeks 
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38 ro? Si il 134 


The apparent D/H ratios of cellulose hydrogen increase 
with higher methane content. This indicates that methane 
is more enriched in deuterium than hydrogen. Unfortunately, 
the relative amount of methane in the final products of 


cellulose pyrolysis could not be controlled. 


The hydrogen gas from trials 42 and 43 were collected 
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by diffusion through a Pd-Ag tube which results in pure 
hydrogen. However, they show a rather significant 
difference (10 °/..) in 6D values between the two whereas 
the difference in hydrogen yield is only 1%. The relative 
amount of methane in the final products of pyrolysis and 
the isotopic fractionation associated with it appears to 
be sensitive not only to the average temperature of the 
pyrolysis but also to the cooling rate of quartz after the 
pyrolysis which is almost impossible to control. This may 
partly account for the poor reproducibility of 6D measure- 
ment of cellulose. It was also found that oxygen appears 
to exchange isotopically with quartz (Si0,). Therefore, 
the method of pyrolysis of cellulose in quartz tubes 


cannot be used either for oxygen or hydrogen isotope 


analyses. 
ite2s Pyrolysis of Celiulose in a Nickel Reaction) Vessel 


The second method of obtaining hydrogen and oxygen 
COryrsOtOpiCc analysis) srom/ organic Matter made use otva 
nickel reaction vessel. The principle of the extraction 


method is similar to that previously described by Thompson 


and Gray (91). 


1.2.1. The Reaction Vessel 


The reaction vessel in Fig. 14 is designed to 
separate hydrogen gas from carbon dioxide and carbon 


monoxide produced by pyrolysis of a sample in the lower 
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combustion chamber. Gases produced in this chamber are 
contained by a stainless steel plunger which forms a 

vacuum tight seal with a gold gasket placed on the rim of 
eestainless "Steel insert. At the pyrolysis temperature of 
moc ©, Hydrogen ditftuses through the 1). 5, mm thick: wets eou 
the combustion eineA Tiare and is adsorbed on a coconut 
charcoal trap immersed in liquid nitrogen. 

The upper part of the combustion chamber is water- 
cooled to avoid damage to the teflon vacuum seals. To 
contain the hydrogen a quartz sleeve surrounds the combus- 
tion chamber. A viton O-ring placed between the cooling 
jacket and the quartz produces an air tight seal when the 
brass collar is tightened against the brass ring. The 
addition of a quartz sleeve to the reaction vessel has two 
distinct advantages. 
dings Hydrogen gas diffusing through the vessel wall can 

‘be collected quantitatively. 
rs Since the nickel reaction vessel is never exposed to 
air at nigh temperature, Oxidationvot the vessel 

Surface does not take place and the vessel has an 

indefinite life. 

However, prolonged heating of the quartz sleeve at high 
temperature causes gradual devitrification of the quartz 
and the lower part of the quartz sleeve becomes opaque and 


starts to deform. This does not appear to affect the 


performance of the apparatus but the sleeves are replaced 


after about 30 analyses. 
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ike2yore i Procedure 


ieee c. b. “Organic: Material 


20 to 30 mg of dry organic sample are weighed into a 
nickel boat and transferred to the combustion chamber. The 
gold gasket is positioned and the vessel attached to the 
vacuum line. It is evacuated and then the quartz sleeve is 
opened to vacuum. The lower chamber is then sealed off and 
the quartz sleeve is connected to the charcoal tube. The 
sample is then pyrolysed at 1150°C for 2 hours. Hydrogen 
gas is continuously collected on the coconut charcoal trap 
immersed in liquid nitrogen. During the pyrolysis, two 
pirani vacuum gauges are used to monitor the hydrogen pres- 
sure in the sleeve and any possible gas leak through the 
gold seal. When combustion is complete the plunger is 
raised and the oxygen-containing gases are collected ona 
silica gel trap immersed in liquid nitrogen. Between 
pyrolyses, the combustion chamber is cleaned with a rotary 
brush and then with a nylon rotary brush using acetone as 
solvent and lubricant. The vessel is then flushed with 
acetone to remove carbon and ashes. 

The CO fraction. is converted to CO. and. Anewa nron 
voltage discharge maintained between platinum electrodes. 
The CO is introduced slowly, in aliquots, to the discharge 


tube and a 4000V potential difference applied using a neon 


luminous tube transformer. A characteristic blue discharge 


Lseproduced and CO is converted to CO, gas. When the 
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conversion is complete, the CO. gas formed is combined 
with the CO, fraction from pyrolysis and the yield measured 
manometrically. Hydrogen yields are measured using a 


Toepler pump. 
2.2.2. Water Analysis 


1 eiestleg OQuartbzCapillary, Method 


For water analyses, a quartz capillary’ tube was used 
to transfer water samples into the combustion chamber. 15 
to 20 mg of water are drawn into a quartz capillary which 
is then sealed and placed in the reaction vessel along 
with 20 mg of carbon powder which has been degassed by 
heeting at 1150°-C for 2 hours. ‘The Carbon” and” capi Vary 
are heated to about 500°C for 1 hour to degas the carbon. 
At this temperature the quartz capillary does not break. 
The combustion chamber is then sealed and the temperature 
Batsed to 1150°C. The capillary. explodes and reaction 


between the H,0 and G“sStarts: 


4 
viene. 3H,0 Zin COs ty:GOs eat 3H. - aa 


Care should be taken to use excess carbon to promote 
complete reaction. Insufficient amounts of carbon results 
in oxidation of the nickel wall, which is difficult to 
remove. The remainder of the procedure is exactly the 


same as in pyrolysis of organic material. 
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An alternative method of transferring water into 
the combustion chamber was tried. After degassing the carbon, 
dry nitrogen at about atmospheric pressure is introduced 
into the quartz sleeve. The lower part of the combustion 
chamber is then immersed in liquid nitrogen. After waiting 
10 minutes for the nickel vessel to cool, about 15 ul of 
water are injected from a micro-syringe through a viton 
diaphragm and frozen into the combustion chamber. During 
this operation the remainder of the apparatus is warmed to 
avoid adsorptive losses of the water. It is essential that 
the transfer of the water be quantitative to avoid 
isotopic fractionation. After the water is frozen into the 
combustion chamber and the chamber sealed, the nitrogen is 


pumped out of the quartz sleeve and the pyrolysis started. 


ie oe Result and Discussion 


Nickel has been used ae purify hydrogen by diffusion 
(37,64) .. In-general,. «the .rate.of diffusion .o£, hydrogen 
through palladium is greater than through nickel under 
identical conditions. However, palladium at high tempera- 
ture (>300°C) quickly develops a coarse crystalline 
structure and breaks down after some hours of operation. 
On the other hand, nickel is stable even at 1200°C. 
Furthermore, continuous use of the nickel reaction vessel 


has produced no signs of physical damage or any less 


effective performance. 
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Ene rate Or -diriusion of hydrogen through the nickel 
wall is a function of temperature, pressure difference, 
and thickness of wall. For a constant temperature, the 


rate of diffusion through the wall per unit length is 


given by 
q = const... x Sie rd (r>>d) [43] 
where r = radius of tube, d = thickness of nickel wall 


and p = differential pressure. 

Whatman cellulose was pyrolysed in the nickel 
reaction vessel and the hydrogen which diffused through 
the nickel wall was collected with a toepler pump. The 
yield of hydrogen was measured manometrically. The 
pyrolysis was stopped after 70 minutes when hydrogen ceased 
to be evolved. Fig. 15 shows yields of hydrogen as a 
function of time at various temperatures. It is apparent 
thaty the rate, of diffusten through the wall offthe’ nicket 
reaction ‘vessel is. a functionm.of temperature. |The rate iof 
diffusion increases as the temperature increases. Thus, 
for rapid and complete pyrolysis, a high temperature is 
desirable. However, the melting point of quartz sets an 
upper limit to the reaction temperature. The optimum 
pyrolysis temperature for isotopic analyses was found to be 
1150°c. At this temperature, no significant physical 
deformation of the quartz sleeve was observed. A higher 


temperature rapidly deteriorates the furnace by oxidizing 


the heating element. 
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AGalio0 CC, pyrolysis and’ Separation: of hydrogen by 
qefeuston- are complete,in 1.5 hours. To prevent possibile 
incomplete hydrogen extraction and isotope fractionation 
associated with it, most pyrolyses were allowed to proceed 
for 2 hours at 1150°C. Progress of the reaction was moni- 
tored by measuring the residual pressure of hydrogen in 
the quartz sleeve... Mass spectrometric analysis of the 
pyrolysis products of cellulose and cellulose nitrate show 
no measurable amounts of hydrogen or any other hydrogen- 
containing gases in the vessel, indicating complete 
pyrolysis. Unlike static pyrolysis in the quartz tube, 
diffusion of hydrogen through nickel walls forces the 
reaction to completion. 

To test the nickel reaction vessel method, a number 
of water samples were analyzed and the results were compared 
WiltChecnose: Obtained by conventional methods.» The results 
are shown in Table 8. 

Within experimental error the 6D and 5)" Or valuesmnor 
the standard waters agree with the accepted values. The 
reproducibility of the method is about +1°/.. for hydrogen 


and £0.2°/.4. for oxygen, which is comparable towother 


methods. 
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Table 8 Hydrogen and oxygen isotopic compositions of 
Standard Waters 


(A) 6D values of Standard Waters (°/,.) (SMOW) 


Water Sample 
Method NBS-1 GISP SLAP suw4t 


Carbon Reduction 


(1) syrinde - -189 (2) -418 (1) - 
i yequartz capillary -48(8) =13°7°(2) -424(2) -117(4) 
Uranium Reduction -48 (2) -187(2) -423(4) -117(2) 
(B) 6189 values of Standard Water (°/,.) (SMOW) 


Water Sample 
Method NBS-1 GISP SLAP 


Carbon Reduction 


syringe =i Je Bal2:) = 245 (29) M547. .8 (2) 
CO,-H,O equilibrium Lee -24.7(2) -55.4(2) 
i 


SHW4 is one of the internal working standards 
“value taken from Thompson and Gray (96) 


Number in parenthesis is the number of analyses. 
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Therevare;two distinct applications of this carbon 
reduction method for isotopic analysis of water. One is 
the simultaneous measurement of 6D and 6!80 values of water 
Samples. The other is the measurement of 6§!80 values in 
milligram quantities of water. This method uses 0.015- 
0.02 ml of water whereas the CO,-H.,0 equilibration method 
requires relatively large amounts of water (0.5-30 my 
(2731): 

The results of replicate analyses of various organic 
compounds are shown in Table 9. 

While the method is quite suitable for compounds 
with a low carbon to oxygen ratio in the molecule (C/O), 
it is less suitable for compounds with a relatively high 
C/O ratio. For example, the hydrogen yields of citric 
acid, cellulose nitrate and glycine are 100% whereas those 
of cellulose and dextrose are lower. The reproducibility 
of SD measurements on citric acid and cellulose nitrate are 
excellent (£0.5°/.o.). On ‘the other hand,“reproducibility 
of 6D measurements for cellulose and dextrose are relatively 
DPOOr. 

The low hydrogen yield and poor reproducibility can 
be related to tar formation during pyrolyses. Compounds 
with a high C/O ratio tend to form tars during the early 
stages of pyrolysis which condense onto the relatively cool 
upper parts of the combustion chamber and are never com- 
pletely combusted. The problem was partially solved by 


adding about 50 mg of ancl. to the reaction vessel along 
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Table 9. 6D and §!°0O values for some selected 
Organic compounds 


Sample H,O Yield 6D Gin2O: wuG/O Seatac 
(3) (°/ooNBS-1) (°/o.SMOW) 

Celi (5)2 94+] =—3 725 450 at 0e2 6/5 

+ZnC1, (8) 96+1 = 2559) sO ht0 3. 6/5 

dext (17) 96+1 AT At 35 7 oe fF 0s a 

+ZnCl, (14) 99+] 59.3120) ale 0 1 

citric(4) 100+1 2952085 - 67a, 

CN (8) 10247 =5 7 10e5 ~ 67 ial 

glycin (8) LOO2 a 49+3 - i 


*Number in parantheses is number of analysis 
Cell = cellulose, dext = dextrose, and citric = 
emer vo acid 


CN = "cellulose nrtcrare 


Lewis acids such as ZnCl. and HgCl, have been used 
extensively in the nydrocracking of ‘coal (ili) HCl. trom 
the thermal reaction between cellulose and ancl, catalyzes 
the nonglycosidic condensation and charring reaction of the 
Very qeanly stades of neating, (300° C)) (62) sehurener shearing 
of nonglycosidic condensation products leads to carboniza- 
tion. Thus. the role of znCl. is to start the decomposition 
of cellulose at low temperature and to prevent condensation 


of volatile products on the upper parts of the reaction 


vessel. 
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Addition of ZnCl. increases both hydrogen yields 
and 6D values, indicating considerable isotopic fractiona- 
tion during the early stages of pyrolysis. For example, 


the D/H ratios of cellulose with ZnCl. present were 26°/,4, 


higher than without ZnCl. when the hydrogen yield increases 


by only 2%. Since virtually no oxygen is included in the 


tars, 6'°O values are unchanged when ZnCl. is added. This 


2 
also demonstrates that neither extraneous oxygen nor 


moisture from the air is introduced with the ZnCl. even 


2 


though ZnCl, absorbs moisture readily. 


2 


The use of ZnCl. however was not entirely satisfac- 


2 
tory, particularly with cellulose, as yields never reached 
100%. The relatively high standard deviation of replicate 
analyses reflects the uncertainty of the method. 

If hydrogen isotopic analysis alone is desired, tar 
formation in conjunction with the pyrolysis of compounds 
of high C/O ratios may be reduced by introducing either a 
small amount of oxygen gas or other oxidizing agent with 
the sample into the reaction vessel. 

Nitrogen does not appear to interfere with the 
quantitative extraction of hydrogen or oxygen from organic 
samples. No attempt was made to separate the CO and CO. 
fractions from nitrogen or nitrogen oxides when the nitro- 
gen containing compounds were analyzed. If nitrogen oxides 
were formed during the pyrolysis, oxygen isotopic frac-— 
tionation between CO, CO, and nitrogen oxides would be 


inevitable. However, the mass spectrometric scans on the 
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pyrolysis products of glycine and intronaphthalene showed 
no traces of nitrogen oxides such as NO or NO... This 
suggests that pyrolysis results in nitrogen being in the 
molecular form, N,- This makes it possible to measure the 
oxygen isotopic composition of the nitrogen-containing 
compounds. 

Table 10 shows the results of test for memory effect. 


Table 10. Memory effect of pyrolysis of cellulose 
nitrate in nickel reaction vessel 


Date Sample 6D (°/. .NBS-1) 
uvuly: 26 Whatman =O 
Aug. 1 Whatman = 5 
Aug. 4 Spruce (1) Salles! 
Aug. 10 Spruce (ITI) -128 
Aug. 14 Whatman =) 
Aug. 20 Spruce (III) -90 
AUG ao spruce (IIT) -91 


Alternating pyrolyses of cellulosesanitratés-of ;differ- 
ent isotopic composition show no memory effect in the 
respective 6D values. Similar tests were carried out using 
two cellulose samples differing in D/H ratios by 170°/oo. 
When the vessels were cleaned between analyses with a 
rotary wire brush, no significant memory effects were 


observed in spite of the tar formation. 
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2. Chemical Treatments of Wood Sample 


Zou. —Extraction of Cellulose 


For isotopic analysis, cellulose was isolated from 
woody tissue according to the standard sodium chlorite 
method described by Green (53). The procedure begins with 
the cutting of wood samples to pass 40 mesh in a Wiley mill 
followed by removal of the extractives with 1:1 benzene and 
methanol for 24 hours. Treatment with the organic solvents 
is necessary because the benzene-methanol soluble fractions 
are so depleted in deuterium that even relatively small 
amounts of extractives in the sample might alter the final 
measured D/H ratio (31). As an additional precaution, 
samples are further washed for 24 hours with acetone. The 
extractive-free wood samples thus treated are vacuum dried 
ace 40°C: 

About 1 gram of the dried sample is suspended in 300 
ml of hot water in a 1 litre Erlenmeyer flask containing 
Ssriecor glacwal acetire ‘acid. "Then -H0'-grams of technical 
grade sodium chlorite is added. The flask is stoppered 
with an inverted beaker and heated for one hour at 90°C. 
The whole chloriting operation, including washing, should 
be done in a well ventilated hood because chlorine dioxide, 
generated during the reaction, is very toxic. “After the 


chlorite treatment, fresh portions of acetic acid (3 ml) 


and sodium chlorite (10 g) are added. This step is repeated 


twice resulting in a total of 3 steps of chloriting treat- 
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ment. The final residue should be nearly white and retain 
the woody structure of the original sample. When a satis- 
factory degree of whiteness has been attained, the solid 
residue is filtered on a Buchner funnel and thoroughly 
washed, first with hot water to prevent redeposition of 
oxidized lignin, and then with large amounts of cold water. 


This solid residue is often referred to as holocellulose. 


Zee.) Alkaline Extraction of Holocellulose 


Holocellulose thus prepared contains an appreciable 
amount of non-glucan polysaccharides referred to as hemi- 
cellulose. Hemicellulose consists mainly of xylan and 
manan with small amounts of other polysaccharides, such as 
galactan, araban and glycurons. The isotopic compositions 
of hemicellulose may not be the same as those of cellulose 
due to the difference in chemical structure, molecular mass 
and chemical bonds os the molecules. To eliminate the 
possible fractionation associated with inhomogeneity in 
chemical constitutents of the sample, holocellulose is 
further treated with sodium hydroxide to give a-cellulose 
whichwis the purest torm of cellulose one .can obtain 
experimentally. 

Holocellulose is transferred to a beaker with 50 ml 
of 17% NaOH solution. The mixture is stirred at room 
temperature for 40 minutes. During this treatment, hemi- 


cellulose dissolves in the water phase leaving only 


a-cellulose behind. The mixture is often slightly brown 
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in color because of the presence of some hignaney later in 
the alkaline treatment, extracts become colorless. 

Holocellulose in general contains up to 2% of klason 
lignin. The subsequent alkaline extraction removes most 
ef the Lignin Weecing only traces. This alkaline extrac- 
tion serves as a means of checking whether delignification 
with the sodium chlorite technique was successful or not. 
When the holocellulose, on contact with sodium hydroxide, 
euris fed in color, delignification is not complete. in 
this case, one more step of chlorite treatment was neces- 
sary for the complete removal of lignin. 

The mixture is then filtered by suction on a fritted 
glass filter. a-cellulose thus extracted is washed first 
with 17% sodium hydroxide to prevent redeposition of lignin 
and hemicellulose and then with 10% acetic acid to neutral- 
ize the sodium hydroxide. The cellulose is then rinsed 
with a large amount of distilled water. The cellulose is 
then dried by solvent exchange with acetone on a buchner 
funnel and further dried in vacuo at 40°C for 2 days. The 
dried cellulose was then used for oxygen isotopic analysis. 


For hydrogen isotope measurement, the cellulose LS nitrated. 


3. Nitration Method 
For hydrogen isotope analysis, cellulose and other 
chemicals are nitrated in order to eliminate the problems 


associated with exchangeable hydroxyl hydrogen in cellulose. 
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3.1. Preparation of Nitrating Mixture 


The nitrating acid mixture was prepared according. to 
the method of Alexander and Mitchell (2). This involves 
the slow addition of 404 grams of phosphorus pentoxide to 
1000 grams of cold 90% nitric acid to produce an acid 
mixture of 64% HNO,-26% H,P0,-10% P50... 

PUO0sg camsgor 902) fuming nitricpacid, isutransternced 
Soa Zeliter 3 neck-flask. The flaskgis kept ice cold by 
immersion in an ice-water bath. 404 grams of phosphorus 
pentoxides are added to the nitric acid. During the addi- 
tion.of Pon, the acid is swirled continuously with a 
mechanical stirrer made of a glass rod and teflon wings. 
Care should be taken not to allow too much atmospheric 
moisture to get into the acid mixture. 6 or 7 wide mouthed 
powder funnels are recommended to prevent the pentoxide 
from sticking to the funnel, absorbing moisture and getting 
into,the acid /mixture....With. occasional swirling, the 
solution is complete in a few hours and the acid mixture 
is then filtered through glass wool into a glass-stoppered 
bottlesand._sconeduingawcoot. dark uploce. 

Alexander et al (2) also have shown that the age of 
nitrating mixture has a direct bearing on its nitrating 
performance. The nitrating performance of the acid mixture 
was poor before 3 days and after 2 months of age. Thus, the 


acid mixtures used in this study were aged between 3 days 


and two months after make-up. Deterioration of the acid 
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was indicated by a marked increase in orangecolor due to 


oxides of nitrogen. 


Bec. Nitration of Wood 


Milled wood is nitrated according to the procedure 
described by Goring and Timell (47). Approximately 500 mg 
of extractive-free wood are rapidly immersed in 40 ml of. 
the nitrating acid mixture which has been cooled to -16°C 
in a ethanol-dry ice slush. Due to the violent reaction 
particularly between lignin and the nitrating acid, 
incomplete immersion of wood samples in the acid often 
results in burning of the sample. The reaction is allowed 
to proceed in a refrigerator kept at 0 - 5°C for about 
24 hours. 

The mixture is then cooled to -16°C and filtered 
through a sintered glass filter. Suction is applied and at 
the same time the sample is gently pressed with a stainless 
steel tamp which is roughly the area of the filter. This 
gentle pressing of the sample with the tamp not only 
removes the maximum amount of the acid but also minimizes 
the moist air drawn through the sample. The solid residue 
is immediately washed with a 1:1 mixture of glacial acetic 
acid and water cooled to -16°C. Then the sample is immersed 
in cold water (0°C) and neutralized by adding a small 
amount of powdered sodium carbonate. The nitrate is washed 
three times with water. The nitrate is stabilized by three 


5-minute treatments with 60 ml each of boiling water. The 
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filtrate is then drained, soaked for 5 minutes in 50 Tissor 
methanol and drained again with suction. The resulting 
nitrated wood is dried in vacuo at room temperature and 
is then divided into two aliquots for further analysis. 
HOrTeNe extract ton jor ‘cellulose nitrate one aliquot 
of the nitrated wood is dissolved in 80 ml of acetone with 
vigorous stirring for about one hour. During this process, 
acetone dissolves mostly cellulose nitrate with a trace of 
hemicellulose leaving others precipitated. The viscous 
Mixture is then centrifuged and the clear solution decanted. 
A large amount of cold water is quickly poured into the 
solution to give a fibrous precipitate of cellulose nitrate. 
Cellulose nitrate thus obtained is washed on the glass 
filter with a large amount of water followed by washing 
with methanol. The white fibrous cellulose nitrate is 
further dried in vacuo at room temperature for 2 Ee, 
Unlike nitration of cellulose, cellulose nitrate 
obtained from the direct nitration of wood is often Ligne 
brown in color probably due to impurities, such as coloring 


matter and a trace of nitrated hemicellulose and lignin. 


No attempt has been made to remove this problem. 


Peo wnt racion of Cellulose 


The cellulose samples prepared by the sodium chlorite 


method were nitrated in a manner similar to the arrect 


nitration of wood. 100-200 mg of cellulose is nitrated 


with O20 mil of the nitracing acid mixture cooled to 0; Cuin 
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an ice water bath. Unlike wood, addition of cellulose to 
the acid mixture does not require further cooling with the 
ethanol-dry ice slush. A reaction time of 3 hours is 
Suitttcient” for the nitration of cellulose. The rest of the 
procedure is exactly the same as the procedure employed 


for the direct nitration of wood: 


See Discussion Of Nitration 


The purpose of nitrating cellulose is to eliminate 
the exchangeable hydroxyl hydrogen in cellulose enabling 
the 6D value of the non-exchangeable CH hydrogen in cellu- 
lose to be measured. The nitration process, however, 
involves a series of chemical treatments which may intro- 
@uce secondary isotope fractionation into the system. 
Therefore it is essential to ensure that the nitration 
process does not alter the original isotopic composition 
of CH hydrogen in cellulose. 

Epstein et al (31) have shown that the benzene- 
methanol mixture has, within experimental error, no 
Sionigrcantwerrect ron tne, b/ iM ratio, OL. tNceittrated 
cellulose. 

Two cellulose nitrations with acid mixtures of dif- 
ferent isotopic compositions showed no difference in the 
6D values of the resulting cellulose nitrates, that is, no 
isotope exchange between the hydrogens In CH “group: Or 
cellulose and the nitrating acid (31). The effect of ace- 


tone on the isotope ratios during cellulose nitrate 
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dissolution was examined by nitrating pure cellulose before 
and after acetone treatment (see Table 11). The acetone 
treatment has essentially no effect on the 6D values of 

the samples. 


TableelL lL, 6D values of cellulose nitrate before 
and after undergoing acetone dissolution. 


Sample Nitration +Acetone 
éD 6D 
Cellulose = At 25 lee.) =r hee) 


It is not possible to obtain complete nitration of 
OH group in cellulose without affecting to some extent the 
Original isotopic composition of the cellulose. A longer 
nitration time may increase the degree of substitution of 
OH group but will inevitably degrade the cellulose structure 
and thereby alter the original isotopic composition of the 
cellulose. Thus it is essential to find the optimum 
nitration conditions that give the maximum degree of substi- 
tution but minimum degradation of samples. 

To find the optimum reaction time, cellulose and 
samples of wood were nitrated for various lengths of time. 
The results are shown in Table 12.. The duration vot 
nitration depends upon the type of material to be nitrated. 
The reaction time of one hour appears sufficient to complete 
the nitration of the exchangeable hydrogen of cellulose. 

Similar results have been reported by Timell (98) in 


which one hour reaction times were found to be sufficient. 
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Table 12. The effect of duration of nitration 
on 6D values of the resulting cellulose 
nitrates. 
Sample Duration SDG 7s) 
Whatman -Cellulose CE NOuUr aoe (3) 
3 2 hours -54+1(3) 
vh 3 hours = aga.) 
- >) nours apoae CS) 
spruce-1 24 hours -146 
spruce-1 40 hours -145 
spruce-2 24 hours -140 
spruce-2 40 hours mlSo 
spruce-3 24 hours -143 
spruce-3 40 hours -140 


The direct nitration of woodmills requires a longer reaction 
tame “(2,53 ,97 ,31)47- andra time of 24 hours appears sufficient, 
The same nitration time has been used by Epstein et al (31). 
They also have shown that repeated nitration after the 
first 24 hour-long nitration has no measurable effect on 
6D values of the resulting cellulose nitrate. 

Temperature also has an effect on the optimum dura- 
tion of nitration and degradation of the resulting cellulose 


Nitrate: eethetis (illustrated ingTable 13. DP represents 


degree of polymerization. Higher values of DP correspond 


to increased degradation, therefore a temperature of 0°C 
gives slightly less degradation but requires longer reac- 


tion time because at low temperature the acid becomes 
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Syrupy and is slower to penetrate the sample. 


Table 13. Effect of nitration temperature on 
nitrate degree of polymerization (2). 


Temperature °C DP a(k=75) 
0 L210 
20 1140 
40 1045 


To ensure the completeness of nitration and minimum 
degradation of sample, a reaction time of 3 hours and a 
reaction temperature of 0° to 5°C were adopted for the 
nitration of cellulose samples. For the direct nitration 
of wood, a reaction time of 24 hours and a temperature of 
0° to 5° were chosen. 

The reproducibility of nitration of cellulose was 
tested. Table 14 shows that nitrating acid mixtures 
prepared separately readies essentially the same 6D values 


Of .cellulose. 


fable 14. Reproducibility of cellulose nitration 
for Whatman cellulose 


Nitrating Acid Yield (3%) 6D(°/o. NBS=1) 
Maoxture © (3).* PO =o OD 
Miscture ©. (4) O02 -6.4+0.5 
Misccire f£l1 (2) 102 -4.8+0.1 


*Number in parenthesis is the number of analyses. 


A hydrogen yield of more than 100% may result from 
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nitrations which have not completely replaced all the 
hydroxyl hydrogens in cellulose. To check the degree of 
substitution, the % yield of nitration was calculated 
according to the equation 
We.” Of} Product 
% Yield = ———__________ x 54.6 [44] 
Wt. of Cellulose 
The % yields of nitration of Whatman cellulose were in the 
range of 96-100%, which is comparable with the result 
(96-98%) obtained by Green (52). 

There are two different ways of obtaining cellulose 
nitrate; one is the direct nitration of wood mills followed 
by acetone dissolution (W-CN) and the other is the nitration 
of cellulose which has been prepared by the standard sodium 
chlorite method (W-C-CN). It is important to check whether 
these two different methods produce cellulose nitrate of 
the same isotopic compositions. For this, 5 year groups of 
tree rings of white spruce (WS1) were divided into two 
aliquots. One aliquot of fe group was directly nitrated 
followed by acetone dissolution (W-CN). The second aliquot 
was delignified by the sodium chlorite technique including 
alkaline extraction to produce a-cellulose. The resulting 
cellulose was then nitrated (W-C-CN). One soft wood and 
a hard wood were also nitrated in the same ways. The 
hydrogen isotope results are shown in Table 15. 

It is interesting to note that there is a systematic 
difference in 6D values between the two different nitrating 


methods. Cellulose nitrate prepared by W-C-CN are enriched 
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Table 15. Isotopic compositions of cellulose 
nitrates prepared from the direct nitra- 
tion (W-CN) and the nitration of 
cellulose prepared by sodium chlorite 
method (W-C-CN). 


Sample W-C-CN W-CN AéD 
SDAP /ss) Yielake PODlS7 ey) Yield eyes) 
WS1l-1l ALS ee =~146. (2) i) 
WS1-2 a2 9: (2) =~ o oe 2) LO 
WS1-3 sie (2) -142 (2) 8 
wS1-4 -141 (2) -150 (2) 9 
WS1-5 See Gar) =] 30 ec) 6 
WS1-6 =—P32-(2) =137 2 (2) 5 
WS1-7 =L3o 5a) -148 (2) 9 
WS1-8 “131 (2) =139 (2) 8 
Soft Wood -78+1 (6) 1002 meals a tocal DO Oa) Ma Wie ta i! 
Hard Wood -68+1 (6) 1002 -82+3 (6) 104% 14 


Number in parenthesis represents number of run. 


in deuterium relative to that prepared by W-CN. The 


average difference is about7°/,... The phenomenon is even 


more obvious for the hard wood where the isotopic difference 


is, AST aie: 

The lower 6D value of cellulose nitrate prepared by 
the direct nitration of wood may be attributed to the 
presence of impurities, such as hemicellulose, coloring 
matter and lignin whose 6D values are lower than that of 


cellulose. As mentioned above, the final products of the 


direct nitration of wood are slightly brown in color 
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indicating the presence of impurities. 

In particular, the presence of hemicellulose 
especially manan in the cellulose nitrate from the direct 
nitration may account for this phenomenon. Timell (97) 
has shown that cellulose nitrate produced by the W-CN 
method contains small amounts of hemicellulose such as 


xylan and mannan, and that W-CN contains slightly more 


mannan than W-C-CN. Xylan, when nitrated, becomes completely 


insoluble in all solvents and partial solubility in acetone 
is realized only after the molecular weight has been con- 
Siderably reduced. On the other hand, nitrated mannan and 
presumably gluco mannan are easily soluble in acetone. 

Thus mannan in W-CN, having a low deuterium composition 
would lower the apparent 6D value of cellulose. 

Due to previous treatments such as sodium chlorite 
and alkaline extractions, impurities are largely removed 
from the cellulose which is further purified during the 
subsequent nitration followed by acetone dissolution. Thus 
nitration of a-cellulose shows better reproducibility of 
6D measurements than the direct nitration of wood followed 
by acetone dissolution. 

Incomplete nitration of wood is also a possible 
problem. For example, the hard wood shows a rather large 
difference in 6D values between the two different methods. 
The high yield (104%) of hydrogen for the direct nitration 


of hard wood indicates incomplete replacement of OH hydro- 


gens in cellulose. This may result from an insurerc vent 
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amount Of the nitrating reagent or insufficient reaction 
time. Wood requires a longer nitration time because the 
nitrating reagent must break the woody structure in which 
lignin and cellulose are interbedded (46). Timell (97) 
used 120 hours of reaction time for the nitration of the 
hemlock wood. The optimum conditions for the direct 


nitration of wood depend upon the type of wood samples. 
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CHAPTER IV 


EXPERIMENTAL RESULTS AND DISCUSSION 


of. Identification and Description of Tree Samples 


Three of the tree samples chosen (denoted WS1, WS2 
and WS3) are white spruce (Picea Glauca) from Edmonton, 
Alberta. WS1 and WS2 were collected from White Mud Creek, 
a Site having a small Be eeauicatiya WS3 was collected 
from a flat site with no immediate source of water in the 
vicinity. White spruce is a natural species in the 
Edmonton region. WS1 and WS2 are believed to be natural 
stance while WS3 may have been transplanted at an early 
age. These tree samples are listed in Table 16. There 
are also 5 conifers from various growth sites in North 


America. 


2 Hydrogen Isotopic Analysis 


To investigate the functional relationship between 
6D values of tree ring sequences and temperature, 6D 
values of cellulose nitrate in tree rings were determined 
for the three white spruce samples, WS1, WS2 and WS3. 

Each tree was divided into five-year groups of growth 
rings. Thus, the D/H ratio of each group represents the 
mean D/H ratio of the five consecutive tree rings. The 
isotopic composition of the five year average is therefore 


biased toward that of the wider rings. No attempt has 
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been made to normalize the amount of wood contributed by 


each of the rings to the materials Subsequently analyzed. 


261: Juvenile Rings and Mature Rings 


inkag. L67bhe oD values of five year groups of the 
spruce samples WS1, WS2 and WS3 are plotted against the 
corresponding growth periods, and compared with mean annual 
temperatures of Edmonton. The general patterns of SD 
variations for the three spruce trees are in relatively 
good agreement with the variation of the mean annual 
temperature. 

One striking feature is that the early growth rings, 
referred to as juvenile rings, of all three tree samples 
have lower 6D values than the mature growth rings indicat- 
ing the possible existence of two different isotope scales 
within a single tree. For example, the juvenile rings of 
WS1 and WS2 from the periods 1883 to 1903 and 1888 to 1903 
respectively have lower 6D values than the mature rings by 
an average of 15°/,4 and<the juvenilewrings) of WS3 are 
more depleted in deuterium than the mature rings by about 
30°/,.- This phenomenon can be seen more clearly by 
plotting 6D values against temperature. 46D values of 
cellulose nitrates for the three spruce samples are plotted 
against mean annual temperatures in Figs. 17, 18 and 19. 
From these figures, it is clear that there exist two 


different ring groups and that each of them has different 


isotope scales. 
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Figure lo. 
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D) Mean Annual Temperature 


~170 


2 
TEMPERATURE (°C) 


1883 1903 9273 1943 1963 198 


YEAR 


Comparison between 6D values of cellulose 
nitrate and mean annual temperature obtained 
from the meteorological data A) WS1, B) WS2 
and C) WS3. * is the boundary between 
juvenile and mature rings. 
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(A) Mature Rings 


(B) Juvenile Rings 


.60 2.90 8.00 3.00 


2.00 2.30 2 
TEMPERATURE (°c) 


The relationship between 6D values of non- 
exchangeable hydrogen of cellulose (cellulose 
nitrate) in tree rings of WSl and mean annual 
temperature (Sept of previous year to August of 
growth year). A) Mature rings, B) Juvenile 


TINGS. 


ec cnins 


Pees a 
rie Pas 
- ‘ 
, - 
aa 


ig 
i 20 #eauUl av 470. fee 
ee nnd aL dees De 
Beant les) A a og ( are 
: _ Seuans nese: ine peo 


(& ,Reaia "aut o 


=t10 


mre 


=FU 


-140 


(A) Mature Rings 


(B) Juvenile Rings 


=190 


-160 


—170 


Figure 18. 


2.00 2.30 2.60 2.90 3.20 
TEMPERATURE (°c) 


The relationship between 6D values of cellulose 
nitrate in tree rings of WS2 and mean annual 
temperature (September of previous year to 
August of growth year). A) Mature rings, 

B) Juvenile rings. 
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) Mature Rings 


(B) Juvenile Rings 


2.00 2.30 2.60 2.90 3.20 3.50 
MEAN ANNUAL TEMPERATURE (°c) 


The relationship between 6D values in cellulose 
nitrate of WS3 and mean annual temperature 
(September of previous year to August of growth 
Year.) vd) Matbure, rings ~.B)suventle wings. 
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Wood Samples 


Sample 


WS1l (white spruce) 
WS2 (white spruce) 
WS3 (white spruce) 
RW1 (redwood) 

YSP1 (white spruce) 
YSP2 (white spruce) 
FVS (white spruce) 


Vie (cedar) 


Growth Period 


LES 5-196 


TSO 0= [9 G5 


16950 Sie ous 


1890-1974 


Del OOS 


E90 = LOT 


10a S975 


1S se Ce oes: 


100 


identificationWand’ Deseripteion of 


Location 


White Mud Creek, 
Edmonton 


White Mud Creek, 
Edmonton 


Mill Woods, 
Edmonton 


Grant Grove, 
Calviorni a 


Decoeli Mt., 
Yukon 


New Kluane Lake, 
Yukon 


RorteVermails! jon, 
Alberta 


Vancouver 
Island, B.G. 


The two groups of isotope data cannot be regarded 


as random scattering of isotope data for a number of 


reasons. 


i 


yA 


Each group consists of the consecutive growth rings. 


The difference in 6D values between the juvenile and 


mature rings are in the same order of magnitude as 


the overall 6D variations shown by mature rings and 


therefore cannot be taken as scattering of data. 


All three spruce samples from Edmonton show the 


same phenomenon regardless of their age. 


This large difference in 6D values between two ring groups 
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does not coincide with any climatic event happening at a 
particular time. For example, the low 6D values of WS3 
for the period 1903-1913 are not reflected in the 4D values 
of WS1 and WS2 for the equivalent time period. There was 
no record of unusual or sudden change in climate between 
these periods of time. The answer may be found from the 
physiological factors controlling growth of trees. 

During the life of the tree there are normal varia- 
tions in the structure and the chemical composition of a 
cell wall, which are associated with the development from 
the juvenile stage to that of maturity (15).. For example, 
the average cell length increases through the successive 
growth rings from the center of the tree outward until a 
more or less constant value is reached. There are also 
changes in chemical composition as the tree grows. It has 
been shown that the amount of cellulose increases through 
thier juvenile“ period (101). "Fall these* changes* are” such 
that as the age of the tree increases the wood formed 
becomes more uniform and shows a greater degree of 
stability throughout the period of maturity than that 
observed in the juvenile stage (15). Due to the isotopic 
competition between chemical constituents in wood, changes 
in structure and relative amounts of chemicals will inevi- 
tably alter their isotopic compositions. Thus, such 
differences between juvenile and mature sections of tree 
rings may give rise to the large difference in 6D values 


observed. Dendroclimatology encounters a similar problem 
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associated with juvenile rings which are often eliminated 
from the analysis because they provide the least reliable 
elimatic information -(44)-: 

To check whether this large difference in 6D values 
is truly related to changes in physiological properties of 
plants, growth ring widths of the three spruce samples are 
plotted against the growth periods in Fig. 20. The ring 
Wrath In’Frig. 20 isan “arbitrary unit which shows. the 
relative size of the ring width. The data were smoothed 
with an 11 year running average which was chosen arbitrarily 
to illustrate general growth trends. Juvenile rings in 
dendrochronology are the first 10 to 30 year portion of 
growth wings, which often show increase in ring’ widths or 
the three spruce samples, the juvenile sections of tree 
rings are well matched with the time periods which show the 
lower §D values of cellulose nitrate. Thus the sudden 
drop in 6D values of cellulose nitrate in early growth 
rings appears to be a unique characteristic of juvenile 
Kings Of white spruce. ‘Thus, ring width, datarcanbe ca 
useful tool to differentiate the juvenile rings from the 
Mature. No attempt has been made to find whether other 
species have the same phenomenon associated with juvenile 
Linds. 

In the functional analysis with climatic data, there 
is no doubt that ee two. rang groups. representing dit te- 
rent isotopic scales should either be treated separately 


or the juvenile ring data be discarded completely. In 
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this study the isotopic analysis of juvenile and mature 
rings were done separately. 

Examples of this sudden drop in 6D values associated 
with juvenile rings may be found in the literature. For 
example, fYapp,eb.al.(108).measured,.éDianded!°Cc values_of 
a branch of a tree from Two Creek Forest in Wisconsin. 
Most of the 6D values varied around a value of -125°/,,. 
However, the oldest 5-year group sample had a 6D value 
L8°%/e5 lower than this. This was not substantiated by the 
673C data. They attributed the low 6D value of the oldest 
sample to utilization by the tree of ambient precipitation, 
rather than the water from Lake Chicago, due to the low 
level of glacial Lake Chicago during the Valders advance. 
However, the low 6D value may just be the unique property 
of juvenile rings mentioned above. 


ive Relationship between Temperature and 6D Values 
of Cellulose Nitrate 


The choice of growing period over which: the tempera- 
tures are averaged is rather ee Most trees have 
a growth period limited to several months rather than the 
whole year. Thus, the functional analysis of isotope 
data with temperature involves the selection of the proper 
temperature group in a statistically meaningful manner. 
For this purpose, linear regression analysis with 6D 
values of cellulose nitrates and various growing season 


temperatures was carried out. 
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Table 17 shows the correlation coefficients between 
6D values of cellulose nitrate and various mean monthly 
temperatures. Proper interpretation of the correlation 
coefficient should include a confidence level which pro- 
vides information about the probability of the relationship 
found being pure coincidence. The confidence level of the 
eorrelatronm coefficient is a function of the number of 
datum points; the more the number of datum points, the 
higher the level of confidence for a given correlation 
coefficient. The 95% confidence level is widely used to 
specify the significance of the functional relationship 
sought (109). The correlation coefficient equivalent to 
953 confidence level in this analysis is 0.58. In other 
words, there is a 5% chance that the functional relation- 
Ship found is purely coincidence. Any correlation 
coefficient greater than 0.58 may be regarded as showing 
a significant correlation. Despite the fact that trees 
Abebamonton, which is 4 ery dry region, grow in 
spring and early summer, it is apparent that the monthly 
temperature of December and January of both the previous 
and the two previous years show the best correlation with 
SD values of the tree rings. The monthly temperatures of 
September, October and November are poorly correlated with 
6D values of cellulose nitrates for all three spruce 
samples. This suggests that winter temperature plays an 
important role in determining D/H ratios of cellulose 


produced in the subsequent growing season. 
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Table I). ~2 Correlation” Coefficients of the iLainear 
Regression Function between §D values of 
Cellulose Nitrate and Mean Monthly 


Temperature 

Sample Year J F M A M J J A S O N D 
WS1 n Some coe ol. 29IR2 OF 38 lS 19 

aa s Git LOVE 24 18" F498" 20.7.9 Ze (oe 13 48 

New 24° 410 jo 19 
WS2 n Oe = Heeee0o 584 "26" 52 745 58 

nad Some 22) 19°. 565-30 Soeee 2) gl Cees 

n=2 694 43 Die oo 
WS3 n > i ima 0 Symes’ See eee 

ge SaaS LON oS”) Soe wee jm eerie! 3:6 

n-2 ToS 4a Oe OS 
Note: n = Contemporary year 

n-l = Previous year 

n-2 = Second previous year 


The importance of winter temperature is also seen in 
the correlation with seasonal temperature. Table 18 shows 
the correlation between 6D values of cellulose nitrate and 


various seasonal temperatures. Winter temperature is 


defined as the average temperature calculated from December 


of the previous year to March of the contemporary year. 


Summer temperature is defined as the average temperature 


of April, May, June, July, and August, which represents the 


growing season temperature. The samples WS1 and WS2 show 
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Table 18. Correlation between 6D Values of Cellulose 
Nitrate and Seasonal Temperatures (éD=at+tb) 


Sample 5-year tempera- Slope Intercept Correlation 
turexqroup Gy he eek vee Ne coefficient 
oe EE A ae Od ee See emer ease ace a 18 eT 
WS1 Winter (Dec-Mar) t Dice 4 Ope Oro. 
| Summer (Apr-Aug) 7 WA CPs 6) =Teo@ 0.47 
WS2 Winter (Dec-Mar) 4.6+0.6 -74+46 OF..92 
Summer (Apr-Aug) OF #3 50) a3 ise 42 0.64 
WS3 Winter (Dec-Mar) Sg cpead ie i) =8/7i14 OS65 
Summer (Apr-Aug) TOL3 =2 036 0.74 
i winter (Dec-Mar) = Average winter temperature calculated 


from December of previous year to 
March of contemporary year 


Average Summer temperature calculated 
from April to August of contemporary 
year 


pormer (Apr-Aug) 


excellent correlation with winter temperature in contrast 
to the rather poor correlation with summer temperature. 

On the other hand, WS3 shows a slightly better correlation 
with summer temperature than with winter temperature. 
Nevertheless, correlation of 6D values of WS3 with winter 
temperature is significant and is comparable with the 
correlation between 6D values and summer temperature shown 
by WS1 and WS2. Most conifers in the Edmonton region grow 
in spring and early summer, so that it might be expected 


that SD values of tree rings would be better correlated 
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with this growing season temperature. However, excellent 
correlation with winter temperature and the poor correla- 
tion with seasonal temperatures shown by the three spruce 
trees indicates that the direct effect of the growing 
season temperature on the hydrogen aR: fractionation 
during photosynthesis is not the major factor in determining 
theiD/H ratios*of.<tree rings. It appears: that D/H ratios 
in cellulose of tree rings reflect that of ground water 
taken by the plants... The isotopic composition of soil 
water taken by plants in turn may depend upon the mixing 
ratio of winter and summer precipitation. For example, 

the high correlation between SD values of cellulose nitrate 
and winter temperature shown by WS1 and WS2 may be attri- 
buted to the importance of winter precipitation in forming 
the soil water taken by the trees. 

It is unlikely that photosynthesis leading to pro- 
duction of cellulose is vigorous during the winter season 
especially December and ey However, it is possible 
that winter precipitation may play an important role in the 
production of cellulose during the growing season. This 
may be the case for WS1 and WS2 which came from Whitemud 
Creek where a natural drainage system is well developed. 
Thus the summer rain at the time of precipitation drains 
out to a nearby stream whereas winter snow, due to poor 
mobility, stays and thereby has a better chance of pene- 
truatingethesground,forsaglong peniodsofetimeshethusythe 


soil water taken by these plants may consist mainly of 
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winter precipitation. Consequently, cellulose produced 

by trees utilmizingithis soil water will reflect mainly sche 
isotopic composition of winter precipitation. On the 
other hand, for WS3 whose growth site was flat, the 
contribution of summer precipitation to the 6D values of 
cellulose might be compatible with that of winter precipi- 
tation. Therefore it is the isotopic composition of soil 
water which plays the most important role in determining 
Dyieratto of cellulose in ‘tree rings. 

However, summer temperatures prevalent during the 
growth year do show reasonable correlation with 6D values 
of tree rings. This may be partly due to the effect of 
summer precipitation, although small, on the D/H ratio of 
soil water. The existence of factors other than ground 
water affecting D/H ratio of cellulose is not completely 
ruled out. For example, temperature-dependent relative 
humidity in the growing season may affect D/H ratios of 
both ground and leaf waters through evapotranspiration, 
and thereby D/H ratios of tree rings. Unfortunately, lack 
of relative humidity records excludes rigorous functional 
analysis between relative humidity and 6D values of tree 
rings. 

Various mean annual temperatures were correlated 
with 6D values of cellulose nitrate in tree rings (see 
Table 19). The mean annual temperature averaged from 
September of the previous year to August of the growing 


year shows a better linear correlation with 6D values of 
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Table 19. Correlation between 6D values of cellulose 
hnhitrate and various mean annual temperatures 
(éD=at+b) 
Sample 5-year tempera- _ Slope Intercept Correlation 
ture group Vie Pe ele yas coefficient 


$$ SSS) 


WS1 Previous year Tree ~155+5 0.80 
(Sep-Aug) 
Contemporary year 5.71.3) tba GO. 79 
(Sep-Aug) ” 
Contemporary year Or O22 le ete) Onde 


(Jan-Dec) ? 


WS2 Previous year deg te) = OLE 4 OER 
(Sep-Aug) 
Contgenpobany  Vears wLO.3+1).:3 -149+4 O93 
(Sep-Aug) 
Contemporary year lee ta. = 19513 £6 O86 
(Jan-Dec) 

WS3 Previous year =O ora 0;..9.3 

E522 

(Sep-Aug) 
Contemporary year 11+4 —15AS5 0.88 
(Sep-Aug) 
Contemporary year 14+3 Sine bse: Uao3 


(Jan-Dec) 


*Previous year (Sep-Aug) represents mean annual temperature 
measured from September of the n-2 year to August of the 
n-l year. 


“Contemporary year (Sep-Aug) represents mean annual tem- 
perature measured from September of the n-l year to 
August of the contemporary year (n). 

3Contemporary year (Jan-Dec) represents the mean annual 


temperature measured from January to December of the 
contemporary year (n). 


cellulose nitrate than the January-December mean annual 
temperature for all three spruce samples. Taking into 


consideration that conifers in the Edmonton region grow 
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from April to August and that December-March winter tem- 
peratures showed excellent correlation with 6D values of 
cellulose, the mean temperature of September-August is 
thought to be a reasonable choice for mean annual tempera- 
ture. 

It is of interest to note that mean annual temperature 
averaged from September of the n-2 year to August of the 
n-l year has the best correlation with 6D values of cellu- 
lose. In other words, the mean annual temperature of the 
previous year plays a more important role in determining 
D/H ratios of cellulose in tree rings than the mean 
annual temperature of the contemporary year. This suggests 
that the plants analyzed here utilized meteoric water 
stored in the ground and that the time lag between the 
aepualliprcecipitation, and fixation by the: plants is about 
1 year. Thus the proper grouping of the mean annual 
temperature that shows best correlation with 6D values of 
tree rings may depend upon rene at factors such as 
continental or marine environment and other factors such 
aselatitude,, almitude .arywom wet, clinateserc ~The tack 
that the 6D correlation with mean annual temperature of 
the previous year is better than that with the contempo- 
rary year may be an artifact caused by the use of 5-year 
groups of rings. However, the fact that all three spruce 
samples show the same trend makes this possibility less 


likely. This will be further discussed with oxygen 


isotope data. 


cif 


SE uA Fas ta As ve a inom eft 
moo Lavine ibang, eokoite: éidanoanss 6) od ps 


a 
— 


he - r) 


> 


Exagmed Lanne msei ae ston oF sestedal 26 ‘eh. 
wpeA, Oa io G+ suis io vadetedan® moxt & 

~iiieo: Bo sen tev\Gy) ee adele tas | , atta ais 
ij io axusssb6ame2 ei | Saori eas , sow XeRtO cas = 


aa 


(Sia ayalg say | 


ne ea li ea" ry »eiodiaag 03 ois Se. ninmaaieln 
054 Sut began ‘eee host iene etasig 4 

at apes arts send bad BSbOD ads, ak 
bg sri add i taxi+ Sas noid 
Was con era to Dapaeebly. imam wie aeet 


es tun aaoga6 “eins tons Yanetant 4m emai, 
ioue s*orostessig6! bate Stipnnorsvee anitse« to (adasnidaos 
jos? SAD) Ado tt Php dey $0 nb yotwteets oath ant 4 
14 Saabyaie anna? Sepitie mae Sa meLagiaraae ge ‘ada 3 
—oanehion stig ity +sile nisi ‘yavdod al sEey ws. 4 
xesy-? 10 1a od. Ye BOM at age am et yaa 
Spite. sents Tis suet ‘Sea3 Ee i Sees -mpcis 20 — fi 
Sei viilldbesog- eth Bacien bret settee ks : Woda: antigen 
nsuyen Adie ileeaaae sors uh, a Ew must alent 


_ 


The linear regression equations describing the 


relationship between 6D values of cellulose nitrate and 


mean annual temperature are 


6D = (7.342)t - (155+5)°/,, for WSl [45] 
CUm healt etl 5624))° 2. forenc? [46] 
6D = (1522)t>=- (16528)°/,, for WS3. [47] 


Both the slopes and the intercepts of these equations are 


adtviferent. In addition, the 
three spruce samples for the 
are not identical. This may 
paleoclimatic reconstruction 


different growth periods are 


absolute 6D values of the 
time equivalent ring groups 
introduce serious problems to 
especially when two trees of 


used. Despite the fact that 


WS1 and WS2 are from essentially the same site, the 


average 6D values of WS1 are 


lower than those of WS2 by 


14°/,, whereas WS2 and WS3 have essentially the same SD 


values. These may partly be 


attributed to the difference 


in the’ D/H’ ratio of the” soil water taken by the tree 


samples. . The mixing ratio of summer to winter precipita- 


tion in the formation of soil water may depend upon factors 


such as lithology, permeability of the soil and topography 


of the®qrowth Site." Due™to differences in“ such*tactors, 


the relative amount of winter precipitation in the soil 


water taken by WS1 might have been higher than that taken 


by WS2. Thus the average SD 


values of WS1 would be lower 
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than that of WS2 and the temperature coefficients diffe- 
rent. In Support of this, it is interesting to recall the 
6D values of WS1 were better correlated (0.89) with winter 
temperature than with mean annual temperature (0.80) 
whereas WS2 shows better correlation (0.95) with mean 
annual temperature than with winter temperature (0.92). 
Since they reflect two different temperatures, both the 
slope and the intercept are also different. 

The number of 5-year ring groups in juvenile 
Sections Of WS1, WS2 and WS3 are 4, 3 and’2' respectively. 
Thus, the number of ring groups may be too small for any 
meaningful functional analysis between 6D values of juve- 
nile rings and climatic parameters. However, they appear 
to exhibit some important trends that are useful in inter- 
preting isotopicratios in tree rings. Table 20 shows the 
correlations between 6D values of juvenile rings and 
various temperature groups. Seasonally, 6D values of 
juvenile rings have better Ce eon with winter tem- 
peratures than with summer temperatures. This is in 
keeping with the result obtained from mature rings, 
namely, that trees utilize mainly winter precipitation 
stored in the ground. 

WS1 and WS2 show that 6D values of juvenile rings 
are better correlated with the mean annual temperature of 
the contemporary year than with the mean annual tempera- 
ture of the previous year. This is in contrast to the 


result obtained from mature rings. The result may indi- 
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Table 20. Linear correlations between 6D values in 
cellulose nitrate of juvenile rings and 
various temperatures 


aay -ecmperature group Slope Intercept iG 
WS1 

Annual (previous) 6. 621...9 =a 5 0.912 
Annual (contemporary) Boal 6 oS ie taal Ss) 0.97 
Winter (Dec-Mar) Oe Wis aes -9454+17 0.93 
Summer (Apr-Aug) ZU 2 ~402+143 Ona 
WS2 

Annual (previous) 290212 ~£502434 0.16 
Annual (contemporary) k5#3 -183+8 0.98 
Winter (Dec-Mar) 54422247 =91424 Cheeks 
Summer (Apr-Aug) =~LBe29 +730¢2355 =O, 52 


Annual (previous) represents mean annual temperature 
measured from September of the n-2 year to August of 
the n-1l year. 


Annual (contemporary) represents mean annual temperature 
measured from September of n-l year to August of the 
contemporary year (n). 


Winter (Dec-Mar) represents winter temperature measured 
from December of the n-l year to March of the contemporary 
year (n). 


Summer (Apr-Aug) represents summer temperature measured 
from April to August of the contemporary year (n). 


cate that there are compositional and chronological diffe- 
rences between soil moistures taken by juvenile rings: and 


Mature rings. For example, they utilize two different 
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soil moistures whose isotopic composition and time of 
peccipitation,are: not) identical. »iWhen trees. are) young, 

the roots are shallow and use soil moisture from relatively 
Shallow strata. The shallow soil moisture in turn may be 
composed mainly of contemporary winter precipitation. On 
the other hand, the fully developed roots of mature trees 
absorb the soil moisture stored in deeper strata in addi- 
tion to shallow soil moisture. The soil moisture in 

deeper strata may be composed of precipitation from pre- 
vious years. Then, it is expected that juvenile and 
Mature rings within a single tree may contain two different 
isotope scales that represent two different mean annual 
temperatures. However, the difference in the source of 
soil moisture cannot account for the large difference in 

6D values between juvenile and mature rings because the 
annual yariation in 6D value of water is relatively small. 
Since the number of juvenile rings is few, compared to 
mature rings, the era ee of 6D values of juvenile 
rings may be fortuitous. 

As a result of this experiment, it issapparent «that 
there exists an empirical relationship between 6D values 
of non-exchangeable CH hydrogen in cellulose and mean 
annual temperature and that it is a positive linear rela- 
tionship with reasonable correlation coefficient. 


25 Relationship between D/H Ratios of Tree Rings and 


the Amount of Precipitation 


The possibility of a functional relationship between 
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0D values in cellulose nitrate of tree rings and the 
amount of precipitation was examined. Correlation coeffi- 
cients of linear regression analysis between them were 
calculated. The grouping of average amounts of precipita- 
tion was done in a similar manner to that of temperature. 
Table 21 shows the correlation between the amount of 
monthly precipitation and 6D values of cellulose nitrate. 
Phe waMount Of Precipitation in the month of January of 
both previous and contemporary years shows the best corre- 
lation with 6D values of cellulose nitrate for all three 
spruce samples from Edmonton. Unlike temperature, the 
precipitation of December is poorly correlated with 46D 
values of tree rings. 

Table 22 shows the linear correlation between 6D 
values of tree rings and the seasonal and annual precipi- 
tation. Compared to the correlation with temperature, the 
D/H ratios of tree rings are poorly correlated with the 
AaMOURnGYOL Precipitation an all cases. | forgexamp le, ene 
best correlation coefficient obtained between winter pre- 
cipitation and 6D values of WS2 is only 0.68 whereas its 
correlation coefficient with winter temperature is 0.92. 
Nevertheless, the result of linear correlation between 
6D values of tree rings and the amount of precipitation 
is qualitatively in agreement with the result from the 
temperature-é§D relationship. 

For all three samples, WS1, WS2 and WS3, winter 


precipitation averaged from December to March shows a 
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raple Zl. * Correlation Coefficients of the Linear 
Relationship between 6D Values of 
Cellulose Nitrate and HoptaLy PYreci pila 
tion Amount. 
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n = Contemporary year 
n-l = Previous year 
n-2 = Second previous year 
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Correlation between 6D Values of Cellulose 
Nitrate and Seasonal Precipitation 


oe 


Season 
Winter 
Summer 
Annual 


Annual 


Winter 
Summer 
Annual 


Annual 


Winter 
Summer 
Annual 


Annual 


ar 5) 
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(Dec-Mar) 


(Apr-Aug) 


2 
(previous) 


(contempo- 
rary) 
(Dec-Mar) 
(Apr-May) 
(previous) 
(contempo- 
rary.) 
(Dec-Mar) 
(Apr-Aug) 


(previous) 


(contempo-- 


rary) 


Slope 
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-0.470.4 
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Intercept 
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=12026 
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tannual (previous) represents average amount of 


precipitation measured from September of n-2 year 


to Auguste ofpn=levyear. 


2annual (contemporary) represents average amount of 
precipitation measured from September o£ snaleyeax 
to August of contemporary year (ye 
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better correlation with SD values of tree rings than summer 
precipitation does. This indicates the importance of 
winter precipitation in determining 6D values of tree rings’. 

Unlike temperature, winter precipitation shows a 
Slightly higher correlation with 6D values of tree rings 
than does mean annual precipitation. This may suggest 
that factors other than the 6D value of precipitation also 
play some role in determining 6D values of tree rings. In 
contrast to correlation with temperature, WS1 and WS3 
show slightly higher correlation with the mean annual 
temperature of the contemporary year than with that of 
the previous year. 

In all cases, SD values of tree rings are negatively 
correlated with the amount of precipitation. The negative 
correlation is due to the amount effect of meteoric water 
on its isotopic composition. . Dansgaard (16) has reported 
several possible contributing factors for the amount effect: 
ibe When a given mass of condensing vapor is considered, 

the isotopic composition of newly formed condensates 

becomes progressively depleted in heavy isotopes 

(D and 180) as the cooling proceeds. Since the 

amount of condensate increases with the degree of 

cooling, the amount effect with a negative correla- 
tion between isotopic compositions of water and the 
amount of precipitation may be expected. 

Aus Fractionation may occur due to isotopic exchange 


between the falling drops of condensate and the 
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environmental vapor. This effect accounts for the 
meelatively highs values of Light rain. sine scucn 
case, vapor below the cloud has not yet been exposed 
to cooling processes and thereby retains high 
values. Thus the isotopic exchange with such vapor 
of high 6 values results in an enrichment of light 
rain in heavy isotopes (D and '°0). In the case of 
heavy rain, vapor composition is determined mainly 
by the liguid phase. 

SA Pract ionatiom bywevaporation.ot fallingsduops.: eLow 
humidity of the low altitude air causes a loss of 
liquid phase by evaporation and thereby considerable 


enrichment of the rain. 


The amount effect is found all the year round at most 
tropical stations and in the summertime at mid latitude, 
but «never at polar stations:(16). The negative correlation 
encountered for three spruce samples from FaMGneOn appears 
to be the result of the amount effects particularly that 
of Praetionat ton by high evaporation in the Edmonton region. 
The high evaporation in turn is related to low humidity in 
this region. 

Table 23 shows the correlation between 6D values of 
juvenile rings and the amount of precipitation. Mean 
annual precipitation of the contemporary year shows a 
better correlation with 6D values of juvenile rings than 


the mean annual precipitation of the previous year, which 
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Table 23. Correlation between §D Values of Juvenile 
Rings and Precipitation 


Sample Season Slope Intercept r 

WS1 Annual (previous) 1 423 ar -17.8118 0.65 
Annual (contemporary) 1.7+0.9 — a3 4 0.81 
Winter Vek yao e ip att 0.73 
Summer do) Sa -185+16 0.78 

WS2 Annual (previous) 6 ay i eth t sedl SPP ETs 9 yee ia 
Annual (contemporary) 1.3+0.1 =16622 0.996 
Winter Ls650.4 wel BRS pe 0.96 
Summer 05-9202 -167+4 0.99 


Annual (previous) represents amount of precipitation 
measured from September of the n-2 year to August 

of the n-1l year. 

Annual (contemporary) represents amount of precipitation 


measured from September of the n-l year to August of 
Eneeconcemporary yearn ae 


is in agreement with the result from temperature data. In 
contrast to temperature, summer precipitation shows a 
slightly higher correlation than winter precipitation. 

The most striking feature is the positive linear relation- 
ship between 6D values of juvenile rings and the amount of 
precipitation whereas the mature rings have a negative 
correlation which is in turn attributed to the amount 
effect of precipitation. The causes of the positive slope 
and higher correlation with summer precipitation which is 


not substantiated by temperature data are not known. They 
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could be an artifact of correlation with a small number 
of vdata. 

Due to poor correlation between 6D values of tree 
rings and precipitation, it is unlikely that one can 
obtain any reliable information about the amount of 
Precipitation from’ D/H ratios of tree rings. “However, 
correlation between 6D values of tree rings and precipi- 
tation supports the results obtained from correlation 
with temperature, i.e. D/H ratios of tree rings mainly 
pettlect D/H “ratio Of Soll water taken by the plant. 


2.4 Hydrogen Isotopic Analysis of Materials other 
than Cellulose in Tree Rings 


Woods are extremely complex chemical systems in 
which each chemical constituent has a different isotopic 
composition. Thus it is necessary to find a material that 
best preserves climatic information. Alternatives to the 
use of cellulose may be extractive-free wood, hemicellu- 
lose’ and lignin. Wood extractives are excluded from) the 
isotopic analysis for a number of reasons. Wood extrac- 
tives include a wide range of chemical types and a very 
large number of individual compounds. Some of the major 
chemical types of extractives are 1) terpenes and related 
compounds, 2) fatty acids, 3) aromatic compounts, and 
4) volatile oils (90). The relative amount of wood 
extractives is so small that isotopic determination of 
them is almost impossible. In addition, some extractives 


found in a given growth ring are deposited at different 
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times from the rest of the constituents and thereby des- 
troy the time scale. 

To find a functional relationship between extractive 
free wood and temperature, D/H ratios of extractive-free 
wood of WS1 were determined. To eliminate the problem 
associated with exchangeable OH hydrogen, extractive-free 
wood mills were nitrated. Table 24 compares the SD values 
of nitrated extractive-free wood with those of nitrated 
cellulose (see Fig. 21). Extractive-free wood (denoted 
by wood in the subsequent discussion) has lower §D values 
than cellulose nitrate by approximately 19°/,,. Total 
variation in 6D values of nitrated wood is 7°/,, whereas 
PocteOr wel lulose nitrate is 10°/,'>. (Lt iseof anmterestato 
note that the nitrated wood of juvenile ring (#17) also 
Nacsa low 6D value. From Fig. 21, it is apparent ‘that 
there is no functional relationship between 6D of wood and 
6D values of cellulose. 

Table 25 shows correlations between various 6D values 
of chemical constituents in wood and mean annual tempera- 
ture. Whereas 6D values of cellulose nitrates are well 
correlated with mean annual temperature (r=0.914), SD 
values of nitrated extractive-free woods are not corre- 
lated with the temperature. The extremely low value of 
the correlation coefficient (r=0.045) suggests that the 
hydrogen isotopic composition of extractive-free woods 
does not preserve temperature information. Furthermore, 


it is to be expected that 6D values of whole wood 
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Figure 21. The relationship between the 6D values of 
extractive-free wood from WS1 and the corres- 
ponding cellulose. Both wood and cellulose 
are nitrated to eliminate exchangeable OH 


hydrogen. 
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Table 24. Comparison between D/H ratios of Nitrated 
Cellulose and Nitrated Extractive-Free 
Wood of WSL 


Ring oD CN* 6D NW** 6D Lignin ASD 
OB ee By Ba 
oe is ee °/ 56 
4 = 4 Se =—1 74 33 
5 seiko =158 =202 81 
6 Ae = tb abies) ay 
e = 139 a S38 54 
8 =o =154 = 2.0.0 69 
9 -134 -155 -197 63 
1A) oe Suse y/| =203 69 
sa ak es, =O -194 Se), 
a2 ese -154 —-196 63 
a7? =P O.G =176 =o 24 
* CN = cellulose nitrate 


** NW = nitrated wood 


- Noy. li is a. juvenile ring group 


: ss awe 
A Cy i 
= 1} i ie t% : 
; | ’ leu hae 
1 ‘ 
‘a 4 


hesetzll to qoesee ae . ie 
ag? svi sOSsRe rs 


abi . fkephe oe er eb 


1 


ag + : aoe | 


Table 25. Relationships between 6D Values of Various 
Wood Components and Temperature 


Variable Regression Intercept Correlation 
Par Coefficient Coefficient 
i kay ha © oY Fs Gx) 

eee aioe ee ee Seri) the matorts!) haleneae Al ashe Fe 8 
6Day vs t (9.441.6 -159+4 059.2 
ODay VSit =0..322<.6 -154+7 -0.045 
6 Dany vs Don -0.042+0.25 -161:+34 -0.064 
Lees vs t -19.7+8.5 -144+22 -0.66 
AGD vs t 29+9 -156+25 0.76 

CN = cellulose nitrate 

NW = nitrated wood 

big = slignin 


including extractives would not have sensitive temperature 
information. Thus, the geriaiey of early isotope work 
done on whole wood is in doubt. 

An attempt has been made to investigate the presence 
of a possible thermodynamic thermometer between cellulose 
and lignin. Lignins are complex, crosslinked, three- 
dimensional polymers formed from phenolic units. The 
quantitative extraction of lignin from wood for isotopic 
determination appears very difficult if not impossible. 
However, taking into consideration that extractive-free 


wood is composed of cellulose, hemicellulose and lignin, 
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6D values of lignin may be calculated from 6D values of 
extractive-free wood and cellulose. 

If hemicellulose is assumed to have similar D/H 
ratios to cellulose and that lignin comprises 30% of 
extractive-free wood, the material balance of wood is 
expressed as extractive wood = 0.7 cellulose + 0.3 lignin. 
If the molecular formula of lignin is assumed to be the 
same as phenol (C.HoO), the total hydrogen in one gram of 
wood will be given by 

No. of hydrogen 
0.7 No. of hydrogen in cellulose + 0.3 ims bionia 


molecular we. ‘or cellulose + aemo lect Lanrewet 
Of; lagnin 


a 
Il 


Or 653". [48] 


From the material balance of deuterium, 6D values of 


lignin are calculated 


= + 0. 4 
0.0653 6Dyy = 0-04326D.y + 0.02216D [49] 


ViG 
= = 50 
SD 36 Dag 25Day [50] 


where NW, CN and lig represent nitrated wood, cellulose 
nitrate and Ifgnin respectively. The difference in 6D 


value between cellulose and lignin is given by 


AéD = SD ae OOS Oy tes 3 (SDay7 SD) - oe) 


Both §D values of lignin and ASD values are listed in 


Table 24. Table 25 shows various linear tiunceronat 
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relationships among 6D values of cellulose; Jagnin.s and 
temperature. As is expected from the poor correlation 
between 6D values of aeractive rec wood and mean annual 
temperature, 6D values of lignin, compared to cellulose 
nitrate, are poorly correlated with mean annual temperature. 
The fractionation between cellulose and lignin thus calcu- 
lated shows a positive linear correlation with the mean 
annual temperature, which is contrary to thermodynamics. 
Thermodynamics predicts a decrease in fractionation as 
temperature increases. Furthermore, the correlation co- 
efficient is lower than that of the cellulose nitrate- 
temperature relationship. Even the seemingly high corre- 
lation coefficient (0.8) of the AdD-temperature is mainly 
due to the effect of the high correlation between cellulose 
nitrate and Me tran A better knowledge of the 
chemical formula of lignin would not likely improve the 
correlation between ASD and temperature. Thus, the 
presence of a thermodynamic thermometer in cellulose-lignin 
pair is not guaranteed. 


Aree) Isotopic Comparison between Nitration of Cellulose 
(WCCN) and Direct Nitration of Wood (WCN) 


There are two ways of obtaining cellulose nitrate; 
one is the direct nitration of wood followed by acetone 
dissolution (WCN) and the other is the nitration of cellu- 
lose which has been extracted from wood by the sodium 
chlorite method (WCCN). The preliminary result of nitra- 


tion of WS1 has shown that cellulose nitrates prepared by 
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the direct nitration of wood are systematically more 
depleted in deuterium than those prepared by the nitration 
of cellulose extracted by the sodium chlorite method. 6D 
values of cellulose nitrates prepared by the direct nitra- 
tion of wood mills of WS1 are attached in Appendix 3. The 


precision of these data is +3°/,, whereas that of cellulose 


nitration is +2°75,. The linear regression estimate between 


the cellulose nitrated by WCN method and mean annual tem- 


perature is given by 


= aoe ° 
SD enwEh ates Leskdee: [52] 


The correlation coefficient is 0.69 which is less than 


ena t vor SDuccy7 temperature relationship, (0.0)-. 


The low 6D values of cellulose nitrate by the WCN 
method and the poor correlation with temperature may be 
attributed to the presence of impurities such as coloring 
matter and some type of hemicellulose. Cellulose nitrates 
Prepared by direct nitration of wood are iigh: brownyin 
color indicating the presence of coloring matter and they 
are not as fluffy as nitrated cellulose. Since lipid 
fractions and lignins are drastically depleted in deuterium 
even very small amounts of these incorporated ingene.cel.u= 
lose phase will reduce the overall 6D values of cellulose 
nitrate. 

Nevertheless, the similarity in temperature coeffi- 
cients of cellulose nitrate-mean annual temperature rela- 


tionships between two different nitration methods suggests 
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that two methods measure essentially the same isotopic 
composition of non-exchangeable CH hydrogen in cellulose. 
However, reproducibility of nitration of cellulose appears 
better than Chat of direct nitration of wood. “For “this 
reason, all other cellulose nitrates were prepared by 
nitrating cellulose extracted by the sodium chlorite 


method, 


26 Hydrogen Analysis of Trees from Various Localities 


Table 26 contains cellulose 6D values of trees from 
various growth sites in North America. Thin radial sections 
of tree ring series ranging from 15 to 60 years were used 
for the isotopic determinations. Average values of 6D were 
obtained by taking long series of tree rings so that 
seasonal and yearly fluctuations are filtered out. Care 
was taken to exclude the juvenile rings in any samples to 
eliminate the possible problems associated with them. The 
functional relationship between 6D values and mean annual 
temperature at the growth sites is found to be linear. 


The linear regression equation is found to be 
SDY= 6.0b = 130°) se [53] 


The correlation coefficient is 0.99 indicating there is 

a strong linear correlation between mean annual temperature 
and §D values in cellulose of trees. These data were 
combined with data selected from the literature (28,31) 


and plotted in Fig. 22. Selection was made on the basis 
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Figure 22. The relationship between 6D values of tree 
rings from various sites in North America 
and mean annual temperatures. Each 6D value 
represents the average over several years. 
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Table 26. 6D Values of Cellulose Nitrate in Trees from 
Various Localities in North America and Mean 
Annual Temperatures at the Sites 


Location SD(°/,. SMOW) Temperature °C 
a ee el Ta BE i ee 
Pewvairane, Grove, Callrornia -60 dial 
2. Edmonton, Alberta nae Phpste' 
Jape Olu wermils 11 On, Alberta -141 male Oe, 
4, Yukon I (Decoeli) 5G =AlaS 
Spee Gull ops iEee = eg oe 
Gre sVancouver Island; B.C. 8 elias?) 
iw OCONCO, Wisconsin =8 7 Teen 
Sue Miami, “lorida -2 23 
9. Spring Green, Wisconsin =o Tteal 
10. Houston, Texas -40 20 
bie Kabamalka, Bt. =9'3 8 
DeeweeeuertOrR co 23 26 


of climatic records being available so that mean annual 
temperatures could be calculated. Data from Epstein et al 
(31) were corrected by adding 5°/,,. to compensate for the 
difference in 6D values obtained by the two nitration 
methods. The linear regression equation obtained from 


this reinforced hydrogen isotope data is given by 
Dee oe =) Mo 0. Vaao. s [54] 


Equations [53] and [54] may be compared with equation [55] 


which describes the relationship between mean annual air 
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temperature and mean annual 6D values of meteoric water 
G25 ex 


SDMSES hot aa LOO es ; [55] 


This equation was derived from the equation between §!%0 
values of meteoric water and mean annual surface tempera- 


ture assuming that meteoric water satisfies 
SDNESC OL Onto. [56] 


The §/%80 data used in this derivation were collected from 
coastal and island stations in tropical and’middle, lati-— 
tudes. Thus, the validity of equation [56] may be 
restricted to those coastal stations. Nevertheless, the 
temperature coefficients (5.5 and 6) thus obtained are 
close to that of the meteoric water equation [56]. This 
is in keeping with the finding of Epstein et al (31) who 
reported a linear relationship in 6D values between cellu- 
lose nitrate and environmental water. The linear regres- 
Siow equation was. a slope (of 0.97, anvinterceptgor =27 47-7 
and a correlation coefficient of 0.98. 

The fact that the temperature coefficient is almost 
the same as that of the meteoric water equation indicates 
that the isotopic composition of the soil water is the 
primary factor which controls the D/H ratio of cellulose 
non-exchangeable hydrogen in plants and that soil water 
utilized by plants is essentially the same as meteoric 
water. The intercept of -130°/,, in the equation, 


compared: to —100°7,. Cf the meteoric equation indicates 
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that the cellulose carbon-bound hydrogens preferentially 
concentrate the lighter isotope with respect to the environ- 
mental water by an average 30°/,, which is close to the 


229/o50 Teperted by Epstein et al, (31). 
6 ie Oxygen Isotope Results 


3.1 Relationship between Temperature and 6!°%0O Values 
of Tree Rings 


*80/'®O ratios of cellulose from samples WS1, WS2 and 
WS3 were measured. Each tree was divided into five-year 
groups of'growth rings. (Thus, the 6’°0O value of each group 
represents the mean §!8O value of the five consecutive tree 
fFings. 6°°O values of ‘celf@ulose from WSl, WS2 and WS3 are 
plotted in Fig. 23 together with mean annual temperature. 
In contrast to hydrogen, there seems to be no significant 
shift in 6!8O values between the juvenile and mature 
sections of growth rings for\any of the three spruce samples. 
However, functional analysis of the oxygen isotope data 
with various climatic parameters were carried out only for 
Mature rings to eliminate any such effects. 

To investigate possible functional relationships 
between oxygen isotopic data and temperature, linear 
regression analysis was carried out between 6'°0O values in 
cellulose of tree rings and various groups of temperature. 
The results are shown in Tables 27 and 28. 

In general, §189 values of cellulose do not corre- 


late as well as 6D values with all temperature groups 
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YEAR 


Comparison between 6'°9 values of cellulose and mean 
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TEMPERATURE (°c) 
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Sand 


1963 1983 


annual temperature (September n-] to August of 


contemporary year (n)). A) 
D) Mean annual temperature. 
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Table 27: Correlation between 6'°o Values of Cellulose 
and Various Annual Temperatures (t.) 


5-year Temperature 


Group Relationship CCt 
ee RO LAELONSOI DS oe Fae Com 


Sample WS1 


Previous Year ae 
(Sept-Aug) 


Contemporary Year §!%0 = (0.99940.35)t_+(21.540.97) 0.66 
(Sept-Aug) 


Contemporary Year 6/80 = (1.22+0.46)t_+(20.941.3) 0.62 
(Jan-Dec) 


(1.3640.41)t_+(20.541.15) Oe 7 


Sample WS2 


Previous Year 6) 
(Sept-Aug) 


Contemporary Year 6189 = (0.6940.27)t_+(23.140.78) 0.62 
(Sept-Aug) 


Montemporary Year (6° °O = (0. 7040.38) t_+(23.141.06) gest 
(Jan-Dec) 


(1.0340.28)t_+(22.240.77) 0.76 


Sample WS3 


Previous Year 6*°o = (1.4740.39)to+(20.441.0) Ona 
(Sept - Aug) 

Contemporary Year Ore (0.8640.40)t_+(22.2+1.1) 0275.6 
(Sept-Aug) 

Contemporary Year S220 = (0.7940.58)t_+(22.341.6) 0.40 
(Jan-Dec) 


CG*e= Correlation) cocre1cLent 


Previous year (Sept-Aug) represents mean annual temperature 
measured from September of the n-2 year to August of the 


n=l year. 

~ ts mean annual tem- 
Contemporary year (Sept-Aug) represen 
perature measured from September of the n-l year to August 
of the contemporary year (n). 


Contemporary year (Jan-Dec) represents mean annual tempera- 
ture measured from January to December of the contemporary 


year (n). 
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Table, 28... CoOnrelation, between §'*o Values of Cellulose 


and Various Seasonal Temperatures (6!°O = at+b) 
Sample 5-year Temperature Slope Intercept (fers 
Group od saa core 


a a SR WE SOW Ee Beeeaon 


WS1 Winter (Dec-Mar) OF 5982 0nero 2 ee eee O52 
Summer (Apr-Aug) dst hal 80 Fra Vi) re ish J 0.64 
WS2 Winter (Dec-Mar) Oe 2720S 2702 lid OFS 
Summer (Apr-Aug) O-8920ns4 14,124 5.2 O63 
WS 3 Winter (Dec-Mar) OSL e0reno OSE eons) 0.48 
Summer (Apr-Aug) Oa'3 22 ee 20 eat OZ Orca a 

CC* = Correlation coefficient 


Winter (Dec-Mar) represents winter temperature measured 
from December of the n-l year to March of the contemporary 
year (n). 


Summer (Apr-Aug) represents summer temperature measured 
from April to August of the\contemporary year (n). 


examined. Nevertheless, §!°O values of cellulose provide 
an insight into the role of various seasonal temperatures 
iInfjdetermining the isotopic composition Of creesrings: 
6189 yalues of WSl have a slightly better correlation with 
growing season temperature than with winter temperature. 
This is contrary to the hydrogen isotope results where 6D 
values are better correlated with winter temperature. 
Similarly, in the case of WS2, 6'°O values are slightly 


better correlated with growing season temperature than 
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with winter temperature. 

On the other hand, while the correlation between 6!°o 
values and winter temperature for WS3 is comparable with 
those of WS1 and WS2 the correlation with growing season 
temperature is much lower (0.21). The best correlation 
between 6'°O values and temperature was obtained when the 
5-year mean annual temperature group was chosen to begin 
One yedr prior to the tree ring group. This was also found 
to be true for the hydrogen results. This suggests that the 
oxygen and hydrogen isotopic compositions of tree rings is 
determined mainly by the isotopic composition of stored 
soil water which has a large contribution from the previous 
year's precipitation at any given time. 

However, the effect of seasonal temperatures on 6D 
and §!°8O values of cellulose is somewhat different. Whereas 
6D values correlate better with winter temperature than with 
growing season temperature, 6'°O values have generally com- 
parable correlations with winter and growing season tempe- 
rature. This indicates a possible temperature-dependent 
fractionation of oxygen occurring in the growing season 
during photosynthesis. This possibility will be discussed 


an the mext’ section. 


3.2 Relationship between Precipitation and 6'*°O Values 


of Tree Rings 
Table 29 shows the correlation between 6180 values 
of cellulose and the amount of precipitation. Even though 


all the correlation coefficients are low, it is apparent 


Ret i i or 


i a = 


7? ua on 
oh 


G'S feaN2Sc ittnievaed ena oilirw , base xertdo. . 
jiw eldwwegnon #2 can 103, sxodeteques 19 
aozgeésa thiwornm Adie wotteleates was. ew Ons 1a oF 
‘poise ft eartes 23d. aT) UE anwol sown at ere 
ads oetw besisedo Sh W WEpRweGee | bas seutey, ote ne 7 
pad oF agen es ee SWI wwgme t taunas “a 


i 7 an 
ps 5 


caw tele UCR PHL _ of ct zobxg 
tog <— 7 


2teappue aid? - eatesngts cea hatiial aria “3OR ores 0 
agit 35 wiked ectim orqatoas ‘Hopowbyd bas 


oe 


nee to pordlequeds Sitgodoe. eds yd tiaiemS inte: 
vq ett now sokwediined egsel 2 sed Hobdv yey 
mtd navte Yas os beiapcar:. 

No’ $2 21. tsxecyod Inuoideon Re +oa%%e ads ‘ytevet 


os width  ¢iteee tee serteoumee ei. pealulisn Re seutsv ont ) 


> inl < DS 


daiw eety sh hala ihe Foahien este qavsed eseledzes « 
wo Pileteier sven eopley of's .Siwostequest gonsak er 
. Ls 
Gms? (olen padwosg tte sects Aa ie agotseloxi09 ids 


bi squh-stunemegigs Riche, & sed gaibad eid? orth : 
ih oo 


A 


noseee 2adwore  adg ei re mepyxo 20 nots 10. kde ee 3 


planes ae 


hbeeawoelt Se ifiw goal 


patsy O° S tne vanadate neswied q 


ie 
ae hey né tZ aes ais ee seg ae - 
Higeds cave no. Set hGboeTR Ro snvoms add bas pees, 


ttezsyes ag #1 oe et taste mises — ine fs. 


" fi 


; ~ J 


——— 


— 


; a yo oe a ae Pr a ie 


Table 29. 


Sample 


WS1 


WS2 


WS3 


Note: 


Correlation between 6'°0 of Cellulose and 
Seasonal and Annual Precipitation 


Season 


Annual 


Annual 
Winter 
Summer 
Annual 
Annual 
Winter 
Summer 
Annual 


Annual 


Winter 


Summer 


Annual 


Annual (contemporary) represents mean annual 


(previous) 


(contempo- 
rary) 


(Dec-Mar ) 


(Apr-Aug) 


(previous) 


(contempo- 
rary) 


(Dec-Mar) 


(Apr-Aug) 


(previous) 


(contempo- 
wary) 


(Dec-Mar) 


(Apr-Aug) 


(previous) 


Slope 


-0.44+0.19 
=0- 342052 


ms). Ot bere 


=0 03820. 0%5 


=) S22 OR cc 
oe OSA S faa 8 hye 


ec ce ein 
=0. 0620-06 
=O. 4 oeOrake 

D052 005 


003.7 £009 
(OPED an wy 


Intercept 


308 + 2a. 9 
Phe a ofan SuSAl 


26.42 To 
OF PAN AN: 
PR SA oles Panes | 
28.4+2.4 
OD ieee al Besa 
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25 S205 
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presents mean annual temperature 
calculated from September of the n-2 year to 
August of the previous year (n-l). 


temperature calculated from September of previous 


year (n-1) 


to August of contemporary year 


(Gay fs 


that 61°80 values of cellulose are best correlated with mean 


annual precipitation of the previous year. 


temperature, 


shows higher correlation with 6*°O values than does the 


amount of growing season precipitation. 


however, 


LMICcOnterast to 


the amount of winter precipitation 


This difference 


in the correlation of §!%O with temperature and precipita- 


tion may suggest that in addition to 6°°O of stored soil 
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water, summer temperature may play an important role in 
determining 6,’ *Ouvalueof iceltulosens Thixsuwasenouizoundweto 
beche rcase with D/H ratio values: 

Based on the result that both 6D and 6!%0 of tree 
rings are a reflection of the isotopic ears Ofmsiod. 2 
water, it is difficult to believe that the large difference 
in 6D values between juvenile and mature sections of growth 
rings could be the result of climatic changes which occurred 
at the site since such changes would have affected both SD 


and §'°0 values of cellulose. 


4, Relationship between Ring Width and Isotope Data 


The use of ring width variations to construct past 
climate has received much attention in recent years (44,63). 
Furthermore, based on the fact that in many tree species, 
either early or late wood width is relatively constant and 
that the early and late wood generally have different iso- 
Hopuc scomposi tion. G327 50), Wigley et al (97) have suggested 
that variations in isotopic composition in cellulose of 
tuee Aingsiimay tbe flan party la tfunction, OfMvartauonssin 
rangiwidth. 

To find a possible relationship between ring width 
and isotopic compositions of cellulose, widths of single 
rings of the three spruce samples from Edmonton were mea- 
sured. The tree ring widths can vary not only with fluc- 


tuations in environmental conditions, but also with 


systematic changes in tree age (44). Thus the normal 
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standardization procedure in dendrochronology is to 
eliminate the ring width fluctuations associated with the 
changing age and geometry of the tree. 

Standardization was carried out by transforming ring 
widths into ring-width indices. A growth curve was con- 
structed by taking an 11 year moving average of the ring 
width data. The measured ring widths (Wt) were then 
converted tO ring-width indices (It) by dividing “each 


width for year (t) by the expected growth (Yt), 


Reka We / Vt. EA 


Division by the expected growth removes the trend in growth 
and scales the variance so that it is approximately the 

same throughout the entire length of the time series. Table 
30 shows the equivalent rings. It is apparent that there is 
no significant correlation between the isotopic composition 
(either oxygen or hydrogen) and ring width for the three 
tree samples analyzed. 

In the case of correlation between total ring width 
and 6D values of the equivalent rings, total ring sequences 
including juvenile rings have a negative correlation whereas 
mature ring series have a positive but lower correlation. 
However, the high correlation coefficient of total ring 
sequences has no physical meaning. It is the direct result 
of the fact that the variation in 6D of total ring sequence 
is largely due to the juvenile isotope effect whereas 


variations in $§D of juvenile and mature rings are small. 
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Table 30. Correlation between Ring Width and Isotopic 
Composition in Cellulose 


W vs 6D Weve 76 
6D = aW+b 6°80 = a'Wtb! 


TT ee ES SE ee ee eee eee eee eee 


Sample Sloper (2/7 o5,/.W,) i's Slope); (°/ ae AN) rs 
WS1 (eee 2 Oe25 iggy aS On32 
C2 209 0.07 Pl ee, Oe 
WS2 Ci 1129220 0.39 PARR SUES, O:256 
C2) eG aes O19 ed ei ne OF 0.6 
WS3 (el) 9 soe S45 Ove25 =) 12020 49 0.09 
C2) e328 Or80.5 =— O04 0g. 0.09 


(1) Total ring series including juvenile rings 
(2) Mature rings 


W = Ring width in arbitrary units 


When the juvenile rings are eliminated from analysis, the 
correlation coefficients decrease drastically. This effect 
is not shown in the correlation with oxygen because there 
is no systematic difference in §*°O values between the 
juvenile and mature rings. 

This result is contradictory to the findings of 
Mook et al (77) who reported a negative correlation between 
6D and ring width in an oak tree from Drente, Netherlands. 


It is however difficult to compare their work with the 
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present study because they did not provide the detailed 
analysis. For example, it is not clear whether or not 
their ring width data were standardized and their growth 
rings included juvenile rings. The poor correlation between 
ring width and 6D or 6'°o values of equivalent rings was 
also observed by Gray and Thompson (50). They measured 
early wood width, late wood width and total ring width for 
the ,entire section of a tree. Further, oxygen isotopic 
analyses were carried out on early and late wood from a 
number of selected rings from the tree. In contrast to 
Wigley et al (104), the maximum contribution of ring effect 
to the variation is estimated to be less than 10% of the 
total variation (50). When 5-year groups of rings were 
used, this ring effect tended to minimize even further. 
Table 31 shows the correlation between ring width indices 
and mean annual temperature. Whereas isotopic compositions 
of tree rings are well correlated with temperature, there 
appears to be no significant correlation between tree ring 


width and mean annual temperature. 


Table 31. Correlation between Ring Width and Mean Annual 


Temperature 
(W=aT+t+b) 
Sample Slope (W°C) Correlation, Coetitecrent a) 
WS1 Ol 2+ O07 0.42 
WS2 Aa ees all Om3 
WS 3 -0.04+0.17 O'S 0H 


W = Ring width in an arbitrary Tne ts T = Temperature 
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5. Relationship between §D and §!°O Values of Cellulose 


It is well known that the §!°0O and the 6D in meteoric 


Waters are related by the simple equation (12), 
SY Se. 6 OL 109/52: 


In general, the slope of 8 indicates the simple Rayleigh 
condensation processes under equilibrium conditions. The 
validity of this equation was further confirmed by Dansgaard 
(16) who obtained essentially the same equation by analyzing 
waters from northern hemisphere continental stations. How- 
ever, it is also well known that both the slope and the 
intercept of this equation vary in different regions (54,16). 
For example, meteoric waters at stations such as Nord, 

Alice Springs, Pretoria and Whitehorse, have a slope of 
approximately 5. The reduced slopes are attributed to non- 
equilibrium evaporation from falling raindrops. Thus the 
actual size of the slope may provide important information 
about the degree of kinetic effect at the site of precipi- 
Catron. 

The. same principle can be applied to the leaf water 
system. The leaf water of terrestrial plants is enriched 
oO compared to the soil water and precipitatrons G3 7420)" 
The enrichment is caused by evapotranspiration of leaf 
water. It is a direct function of the relative humidity. 

In order to demonstrate such a relationship between 6D and 
$180 of leaf water, Epstein et al (28) measured the 6D and 


S!80 of the water vacuum distilled from various plant 
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leaves. The slope of 6D vs 6'°O of leaf waters turned out 
to be 2.5. The low slope was attributed to the effect of 
non-equilibrium evaporation on the 6D and §!%0 of water in 
the plant accompanying evapotranspiration through the 
leaves of the plant. 

The question is whether or not the cellulose, the 
final product of photosynthesis utilizing leaf water, 
preserves the same relationship between 6D and 680 
observed for meteoric water. For example, Libby et al (68) 
in their isotopic analysis of an 1800 year ring sequence 
obtained a linear relationship between 6D and §!°0O of whole 
wood. The slope of the equation turned out to be 8, which 
was used as evidence that both 6D and §!°0O of whole wood 
reflect those of meteoric water. 

To resolve this problem, an attempt has been made 
to find a functional relationship between 6D and 6&'°0O of 
cellulose extracted from the three spruce samples (see 
Table 32). The slopes of the linear equations thus 
©btained"tortherthree™Mecamples® are? 2.921372 72266 ,fand 


Bes H2ne respectively.’ * Despite the fact thatmthemthree 
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samples are from the same climatic region, the slopes of 
these three trees are not the same even taking into account 
the large uncertainties in the slopes. 

Taking into consideration that the standard deviations 
of the slopes are large, the slopes of WS2 and WS3 may be 
regarded as essentially the same having a value of approxi- 
mately 6. This may be compared with the slope of 7.7 
observed for meteoric water in the Edmonton region (54). 

The difference in slope between cellulose and meteoric water 
may indicate the degree of kinetic effect due to evapotrans- 
piration occurring on the leaf structures. 

The slope of 2.3 observed for WS1 is more difficult 
torexplain<s In general pjsaclowaslopesofnéDevs 620 rela= 
tionship is associated with differential enrichment in heavy 


180) due to the kinetic isotope ‘effects in 


isotopes (D and 
evapotranspiration in conjunction with low relative humidity. 
Considering that WSl1 and WS2 are from essentially the same 
Site part cis cumbikely cthat relative humidity at the growth 
site of WS1 was much lower than that of the site of WS2. 

The fact that the 6D values of WS1 are systematically 

lower than those of WS2 and WS3 indicates that the incon- 
sistency in the slope of $p-6!°0 relationship between 
individual trees may be elucidated in the context of the 
source of water taken by the individual trees and the 

source of hydrogen and oxygen in tree rings. 


First, the composition of soil water (mixing Hattie 


of summer-winter precipitation) taken by individual trees 
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may be different. Thus soil water taken by WS1 is mainly 
winter precipitation and its isotopic composition is rela- 
tively constant. Since relative humidity in the Edmonton 


region is relatively low and the 6'°0o value of water is 


more sensitive to the kinetic effect due to fast evaporation, 


the slope of 6D-6'*°O equation for WS1 would be small. Thus 
variation in 6D and 6'°0O values of WS1 is thought to be 
caused mainly by relative humidity variation in this 
region. To explain this, a hypothetical growth condition 
may be considered. When a tree stands on a slope, the 
water taken by the tree will be composed of mainly winter 
precipitation because most of summer precipitation will 

run off. On the other hand, water taken by a tree on flat 
or hollow ground will be a mixture of the winter and summer 
precipitation. Then even if the relative humidity at both 
sites is the same, the water taken by the two trees will 
show different slopes provided that summer and winter pre- 
cipitation preserve different slopes. The effect of 
humidity on the slope will also be different for the two 
trees. For example, it may be expected that the slope for 
the former tree would be more sensitive to relative humi- 
dity than that of the latter which is refurbished with 
fresh summer water. Thus the slope of §D-6'°O relationship 
exhibited by individual trees may be different and it may 
not be a simple function of relative humidity. This 
hypothetical consideration may not be far from the actual 


case reported here. In support of this, it is interesting 
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to note that the 6D values of WSl are better correlated 
with winter temperature than either with summer or mean 
annual temperatures whereas WS2 and WS3 show an equal 
contribution from summer and winter temperatures. 

Secondly, the poor correlation between §D and 6'°o 
of cellulose may indicate that if hydrogen and oxygen 
originate from different sources, D/H ratios and }%0/1!o 
ratios of cellulose may be independent of each other and 
thereby represent two different sources of climatic informa- 
tion. For example, unlike hydrogen, oxygen isotopes may 
experience temperature-dependent fractionation during 
photosynthesis. Then the 6D-é'°0O relationship of meteoric 
water would no longer exist in the cellulose system. This 
problem will be discussed further in the following chapter. 
Finally, it is not completely ruled out that the inconsis- 
tent slopes between individuals of the same species are an 
artifact of the poorly defined relationship between 6D and 


Oe Orotmcermilose. 


G Model Study of the Isotope Composition of Cellulose 


As a result of this experiment, it is obvious that 
isotopic compositions in cellulose of tree rings are 
determined not by a single climatic parameter but by 
several parameters such as temperature and: ‘humidity peek 
is worth developing a model which describes the association 
of isotopic variations with variations in climatace Gactors:. 


In isotope modelling, one of the most important 
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procedures is to investigate the source of isotopes in the 
material in question. A clear understanding of the source 
of isotopes may resolve the problem as to what the isotopic 
variation in cellulose reflects, e.g. isotopic composition 
of water, temperature variation at the time of photosynthe- 
sis, and isotopic composition of air etc. As for the 
source of hydrogen in cellulose, it is quite obvious that 
hydrogen comes from water because plants do not absorb 
organic material from the ground and the amount of hydrogen 
gas in air is negligible. Therefore, it is to be expected 
that the hydrogen isotopic composition of water will play 

a primary role in determining the D/H ratios of cellulose 
in plants. However, the source of oxygen in cellulose is 
not an obvious matter. These exist at least two possible 
pathways whereby oxygen can be fixed into cellulose via 


photosynthesis. 


6a Oxygen Isotope Model 


Cellulose oxygen could be derived from atmospheric 


COl-with. or without equilibration with, water.~ Another 


2 
possible pathway involves derivation of oxygen from water. 
Several models based on these assumptions have been 
proposed to account for the oxygen isotopic composition 
of cellulose in plants. For example, Epstein et al (28) 
postulated that one-third of cellulose oxygen in terres- 


trial plants comes from the environmental water and the 


other two-thirds from the atmospheric co. without equili- 
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bration with water: 


N 


$18o Sey ee lewis 
dell 32.0 “co, See 6 °H,0 : 


[58] 


By combining this with the well known meteoric water, 
equation [56], Epstein et al (28) obtained a» relationship 
between 6D and 6'°0 of cellulose. The slope of 6!80 vs. 
6D of cellulose was calculated to be 24 whereas the 
measured slope of the relationship is#l4” (see Fig. 24)). 
The difference in slope between the model and the actual 
data is almost in the order of magnitude. The model 
developed by Epstein et al (28) is not in full agreement 
with the actual data collected by the same authors. On 
the other hand, Ferhi et al (34) proposed a model of oxygen 
isotopes in cellulose based on laboratory experimental 
results. The model developed for the isotopic composition 
of the cellulose of terrestrial plants is composed of 3883 
of leaf water and 62% of atmospheric CO. Which 7S noe Gin 


isotopic equilibrium with either soil water or leaf water: 


oO = ge sere On ett 0-62 6700 [59] 


cell lw CO. 
where cell, lw and co. represent cellulose, leaf water and 
co. respectively. The result is contrary to the result 
obtained by DeNiro et al (19) who have reported complete 
exchange of the CO. with leaf water in the plants of their 


experiment. They undertook laboratory experiments in 


which they grew two sets of wheat plants under conditions 


bes 


+ 7 af i\ a 2iesaw. {: 


Ista ot sont rewcal Stew ota Het nidg “pata 
Adacciisisr = Scipio AAR de, se ubetnge tee) = 
‘4 30 sqote ail aiotulien +0 Ofte ons a4 
eds aces Le wel od hata lgales bid saclut 


.prt esa) Ao ee sitio foe, aie rt oqote ben 
cetee ett ons Lebom galt may ses MORE pid eonate2t6 ‘b 
(abon ox’? shag posi ibe: gabwo silt int deomts 
snsmesips Lio ab sae et (6) Ia as misgega yd & a 
a 

 ,enetans sfte oie ver Dissioe Loe Bete tad on: 


iepyae Jo \Tebdi « baaaaoag re da te Jisat Snad 
S539 =OQsze crOse tates oO beead evolulios at\a 


h, $e. 3] 


Maitigegmtic slad ii oa! wot haigo ove’ ishom ofT fi : 
ae Zo bSeoaneo ae eauehe iaiacesszed %o seolutien ad? 3 

i ton a4’ Hodtiw gBD of uedituouts to Q59 Baw. sotew soll 
:1scne Wesel To ete Uist tet te tokwoeutsdifjupe ods 


= oat 
hd a eee i ty = *0in eget 
on” * ae + ec & 6£.0 Lfeo° — 


r 


hig 2wetow taeol onoliytigs, siaserng: Pe Bas wi ,ifso x 
- ; , > 

2iUBSs a2 OF (ravague ak Sluser ‘edt ‘Ylevisosqeer 
eselahan hat nee evad ont (Gi) ia de. ‘exten vs 


tilad? toa asasigniwid ee rede tant ‘ia bo. 
Li S4aflams tages wpptsiates poussabau You? -30sni 


7 


Pat 


eats »- 7 a 
7 Tr tala hoy 


2a) 
: 


7 
nei 
4 


PSs 
(=) 


oO 


‘ 
> 
oO 


=Be 


' 
QS 
iS) 


-166 


cavi 
nvironmental water, Dewy (Per mi 1) 


Figure 24. 


151 


I} Plamts growing submerged in water 
morine oO Curve ] ane 


non-marine e K 

2) Ail other plants x marine 

3) Gray and Thompson © (30° C) 
Curve 1 
non-marine 


(30°C) 


(per mil) 


Cellulose, 518, 


SMOW 


Relationships between the 6 “Or oteche 
cellulose (from both aquatic and terrestrial 
plants) and the 6D of the meteoric water 
used iby, .the ,plants,..(28)... .Note thatthe 
curve 4 is the simple linear regression 
curve for terrestrial plants whereas curve 

3 was calculated by Epstein based on his 
oxygen -isotope model: 
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Similar in all respects except for a large difference in 
the ~46/}' 0: ratios rot the CO, Supplied to them. Therefore, 
the difference in 6'°O values of the two celluloses can 

be attributed to the difference in 6!°O values of CO, gas 


supplied. This would indicate that the oxygen isotopic 


composition of atmospheric CO. influences that of cellulose 


”, 
on plants. In spite of the large difference in *°0o/*°o 
ratios of the CO, supplied to the two different wheat 
plants; the ?°0/*"oO¥ratios of Gellulose purified from the 
two sets differ by only a small amount. DeNiro and Epstein 
(19) thus concluded that the oxygen derived from CO, under- 
goes complete exchange with the oxygen of the water in the 
plant during photosynthesis and that the '°0/'°O ratio of 
the leaf water is the primary influence on the BIG / ES 
ratio of the cellulose in terrestrial plants. The models 
proposed “by Epstein et al *(2'8) “and Ferhi etal (34) "fart 
to account for this phenomenon. 

Here an attempt has Been made to develop a model 
which accounts for the hydrogen and oxygen isotopic compo- 


sition of cellulose in plants. Photosynthetic processes 


may be described schematically as 


: 60 
co. + 2H.,0 > CH.O + HO + 0. [60] 


Plant photosynthesis is an oxidation-reduction reaction 
between HO and C04, and uphill transfer of four hydrogen 


atoms from H.O to CO, and can be divided into three parts 


(81): 
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1S Enzymatic transformation of HO bo 0. 
+ : 
2H.0 > 2H, 0. [61] 
2 Hydrogen transfer from an intermediate in enzymatic 


sequence 1. to an intermediate in enzymatic sequence 3 
with the help of light activated chlorophyll. 


oe Enzymatic transformation of CO. to CH,0. 


The first step, the enzymatic dehydrogenation of water was 
proven experimentally by Ruben et al (86) who demonstrated 
that the source of 0. produced in photosynthesis is the 
water molecule and not cO.- It is the third stage that 


envolves incorporation of° CO, into organic materials. 


2 
Therefore, it is apparent that the primary source of oxygen 
in cellulose is cO..- It is, however, unlikely that co, 
could undergo direct reduction. Rather co, LSA PLoS 

bound during photosynthesis by a non-photochemical reaction 
to an organic acceptor. Ribulose Diphosphate (RDP) is 
found to be the primary acceptor of CO. and phosphoglyceric 
gaeid (PGA) “is formed as the produce of this: carboxylation. 
PGA thus produced reduces, through a photochemical process, 
Lae panteto. sugar ‘and in part tothe CO. ACCEDLOL Wi RDe. The 
ee atoms of CO. incorporated into PGA may undergo 
equilibration with H,O by means of hydration as PGA reduces 
to carbohydrate through the Calvin cycle (81). The subse- 


quent polymerization of glucose to form cellulose does not 


appear to provide any further opportunity for exchange of 


the oxygen derived from CO.. 
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In this model, it is therefore assumed that the 


oxygen in cellulose originates from atmospheric CO. ehac 


is in complete isotopic equilibrium with water in the leaf 


and that the fractionation between the oxygen in CO, and 


that in the cellulose precursors in the Calvin cycle leading 


to cellulose formation is not temperature dependent. The 


equation for the oxygen isotope model is then defined as 


On = BPA o ce. + E. [62] 
The equation has the form of a linear regression with slope 
of A and intercept of E. The slope A is the fractionation 
factor between the oxygen in cellulose and the oxygen in 
CO, undergoing isotopic equilibration with leaf water. 
Thus the actual value of A may depend upon the type of 
photosynthetic metabolism. This equation can be rewritten 
in terms of a fractionation, factor ate) 


Z 2 


Seon = A(1000 ina + oe + E [63] 


mhnecelcens fractionation factor between CO. and HO, and 
lw denotes leaf.water. ,To.find the,actual value of A and 
E, the result of DeNiro et al (19) was used because this 
is the only experiment which measured 6° "Ou valueswoL both 
leaf water and cellulose. The slope, A, turned out to be 
0.62 and intercept 2.5%. For this calculation, the frac- 


tionation factor 2 was calculated uSing the equation 
CO,-H,O 
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2 3 


oe Ul tke Oe ie Oo. COUT Ey <a loner injhiola ge [64] 
Thus the equation has the form 
a0 oe = 0.62(1000 Liteo.cuset on Paws ar ed aa A 
= 620 acer ee ROMA Rac, oe ty Oo [65] 


Table 33 compares the 6'°0O values calculated by different 
models with the experimental data obtained by DeNiro et al 
(19). As was mentioned earlier, all other models fail to 
account for 6'*O values of wheat plants in DeNiro's 


experiment. 


Table 33. Comparison of Various Oxygen Isotope Models 
with Measured 61°0O Values of Cellulose 


6240 mea=. 6450 620 ne) MELE) 


SM EM FM1 FM2 
Sie ad Cee wets fee is. sae 
oo 
ete ee a eee 
Experiment l IESG PPM ie OE PASE TY 21a 20.23 pra eeal 
Experiment 2 554: 3.0e Les ee DMB es 104.4 30.4 


SM present model 
EM = Model of Epstein et al (28) 


Model of Fehri et al (34) 


FM1 


FM2 = Model of Fehri et al (34) revised 
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61,2 Isotopic Enrichment of the Leaf Water 


In general, water in leaves of land plants are 
enriched in heavy isotopes (*%0 and D) compared to water 
in branch or soil water (45,20,37). This enrichment is due 
to evapotranspiration. The variation in isotopic composi- 
tion associated with evapotranspiration is given by a 


formula derived by Dongmann et al (20): 


es ie — * ous 
qt . (6 {e*x + Ses OS €,)h}) [66] 
where «* = a-l, 6 = the enrichment of leaf water relative 


to branch water, A = the enrichment of water vapor relative 
to branch water, Soe non-equilibrium separation factor 
and. t = Tne) ;,h)- 


For steady state (dd/dt=0), the equation reduces to 


s* He )h [67] 


O 
II 


t (o,7€ 


k k 


= * ras 
See : os + ¢* + Ey i Gor= o's [68] 


where Sow is the isotopic composition of branch water. 
The uptake of soil water by the roots and its transporta- 
tion to the branches causes no fractionation (37). Thus 


the isotopic composition of branch water in equation (68) 


can be replaced by those of soil water: 


a = eee + ¢* + Ey + (S,7€,)h. [69] 


This equation describes the enrichment of leaf water in 
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the heavy isotope due to evapotranspiration. DeNira et al 
(18), however, have shown that the 6!°0O value of leaf 


water is influenced by that of CO,- This effect could be 


a result of either oxygen isotope exchange between CO. and 


H,0 in the leaf or production of H,0 from photosynthetic 


processes in which CHO, dissociates into CH.0 and H,0 


whose oxygen comes from CO.- in natural’ systems; this co 


effect has not been noticed simply because 6!°0O value of 


2 


atmospheric CO, is relatively constant throughout the 
world’ (6) and the magnitude of the effect seems” to *be small 


CUimar /yo) "FL the Cow effect “i's ‘tCaken™intolconsiderationy 


2 
the equation will have the form 


Jh + pa Opa i [70] 
Z 


= * 
Siw cee + ¢* + Ey + (6 


ee 
The numerical value of x calculated from the result of 
Deniro et al. (19) turns: outsto be 0.0143. “In this -cal- 
Culation, “5 "YO Valite ‘of inlet CO, was used fOr ea in 

2 
the equation. Combining equations [65] and [70],. the 


equation for the oxygen model becomes 


oO. = 620 lna + 0.62 1S One ee et 


ell CO.-H,0 


18 
- ON dace mere gin 
Ex + (5, ey, )h + x6 co, El 
To test the validity of this equation, a number of 
parameters such as €*, Eye Sn and h are required. The 


values of e* at different temperatures are well defined. 
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For example, at normal temperature e«* is 9°/,, (16). The 
kinetic separation factor Ex uSimore diffieult+to estimate. 
The value of ey depends upon the type of aerodynamic boun- 
dary layer at which the kinetic separation occurs (20). The 
value of €, appears also to be species dependent. For 
example, deciduous trees such as oak and birch have ey value 


of about 15°/,, whereas coniferous trees like fir and spruce 


Ravetewnvalue of 24°/2.0(3%)).9 As for the valtiefof A, it jis 


k 
important to note that there are two sources of water vapor 
in the environment of the leaf (20): well mixed water vapor 
from elsewhere and the water vapor from evapotranspiration 
of the surrounding plant and soil. Thus two limiting cases 
can be considered. The first is the case where the leaf is 
isolated from other vegetation. Furthermore, assuming that 


atmospheric water is in equilibrium with soil water, the 


A is determined by the equilibrium separation e* 


On See O10 ai 204 C)N. 
In the second case, the water vapor near the leaf originates 
exclusively from evapotranspiration. Under steady state 
conditions, the water vapor has the same isotopic composi- 
tion as the soil water. Thus on = 0. In natural systems, 
however, On will lie between these two extreme values. For 
example, the On value measured in the Edmonton region is 
eee 22), POLSLe Ls (3G) mnas demonstrated rather constant 
On values of about -8°/,, regardless of the species of 
trees. Using these parameters, the oxygen model is cali- 
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brated with the actual data obtained from laboratory experi- 


ments by DeNiro et al (19). Table 34 compares the SIO Ore ae 


measured with Seay See calculated from the model. 


Table 34. Comparison between the Actual and Model Values 
Of 6) "OP OfecCellulose 


18 18 
6) Orell measured 6 Onell calculated 
Experiment 1 ZOO, 2.9 26.8 


Experiment 2 Boe a es ONL 35.4 


The equation thus obtained is given as 


Vi = + 0.62(6'°0_ + 25-23H + 
On oy ee a eo sano OOo ae 5-23 
18 
0.0143 6 Oa Nees Ress) ye 
2 
where € = 9°/,.,, oD ig NG S76. (and SAR =a oe oy 


Untortunately, DeNiro et al (19) *didynor: keep the 
relative humidity constant in their experiménc.”) [he Gela— 
tive humidity varied from 30% at the beginning to 100% at 
the end. Thus the relative humidity of 77% is taken as the 
average during growth experiment, which appears to be a 
reasonable choice. The result obtained from the model is 


in full agreement with experimental results of DeNiro et 
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OFs Hydrogen Isotope Model 


In the case of the hydrogen model, it is assumed that 
there is no temperature dependent fractionation during 
photosynthesis. This assumption is based on the fact that 
the relationship between 6D values of cellulose nitrate 
and mean annual temperature is very similar to that between 


6D values of meteoric water and mean annual temperature: 


6 aye) ee el sOF tae ea 


Peet 


z = OOK 74 
SDE 5.6 ta TOO oe [74] 


In particular, the almost identical slopes indicate that 
the isotopic composition of the source water is the primary 
factor which controls the D/H ratio of cellulose. MThus the 


equation for the hydrogen model is given as 


SD +e* + €, + (6,7-€,)h + E fete] 
sw 


where Fk isthe biological factor. To.calibrate ehis 
hydrogen model, the data of DeNiro et al (18) was used. 


The parameters for this calibration are as follows: 


Sow 4S08/ 56 Hele) 
ie TA) ox Me) 
he= Pe o 

EP SISIONS: 


6D thus calculated are in excellent agreement with the 
cell 
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measured values (see Table 35). 


Table 35. Comparison between the Measured and Calculated 
D/H Ratios of Cellulose 


eee SD measured. 4s oDrcalemlated im 
Experiment 1 = 3 o/ oe Le / ote 
Experiment 2 15/5 mek ayes 2 


The equation thus obtained is given by 


= oD = 30h 92. (76) 


6.4 Calibration of the Isotope Models with Natural Systems 


Finally, the validity -of aes models must be checked 
with the data from natural systems. For this purpose, a 
number of trees from various sites in North America were 
examined. Table 36 contains all the climatic information 
necessary for model calculations. Temperature is the mean 
daytime temperature during the growing season (usually 
Apral to August).« Most of sthe climaticeinformation) was 
obtained from meteoric stations near the growth sites 
(1,80,87,100). However, in some cases climatic information 
was not available for the exact sample site. In this case, 
the data from the nearest meteorological station were used 


with corrections where applicable. 


In applying this model to natural systems, some of 
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the values of the parameters involved in the model equa- 
tion must be chosen. For example, the kinetic separation 
factor Ey for the oxygen model was chosen as 24°/,, (38) 
because most of the trees analyzed were conifers. The 
value of ee is assumed constant at -8°/,,. The term 
et Pc, becomes 0.6°/,, because the oxygen isotopic 
composition of atmospheric CO. is: cOnstentyat 4287 2. 


throughout the world. Thus the working equation for oxygen 


model reduces to 


18 a gi8 if 
6 Ocell Oot Gos shoe FOE Giz. ( Ost 33-6 B22 rey wa beer 
In the case of hydrogen, the parameters used in this calcu- 
lation are as follows: 

ete = 80°/.. 
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Hence the equation for the hydrogen model becomes 


=06D Ut lO4 = Oana e 7 eee [78] 


pO ell SW 


The equation for both oxygen and hydrogen are in excellent 
agreement with the actual data (Table 37) PEGs 25 4Snows 
the strong linear correlation between Onna measured 
and 6/80 calculated from the model. The correlation 

e 


ell 
coefficient is 0.98. Fig. 26 shows the linear relationship 
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Figure 25. 


19 21 23 26 27 29 31 33 35 37 


OVP HO) (%- ) MEASURED 


SMOW 


The relationship between §'°O values of 
cellulose measured and calculated with the 
oxygen isotope model for trees from various 
growth sites in North America. 
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Figure 26. The relationship between 6D values in cellulose 


nitrate measured and calculated using the 
hydrogen model developed for trees from various 
growth sites in North America. 
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Table 37. Correlation between various pairs of isotopic 
and climatic parameters 


parameter pair slope intercept correlation 
coefficient 

eee eee gee See oe 

6D measured vs 6D TO. OZ -0.3+1.6 0.998 
model 

6D cellulose vs 6D 0.88+0.06 -13+6 0.974 
water 

6D cellulose vs 
Summer temperature 10:1 ~306219 O96 0 

6D cellulose vs 13.4 tH: -174+79 0.329 
relative humidity 

epcellultose vs" oO 740.5 -3.6+6.02 0.974 
water 

bpecellulose vss’ °O 14+2 -459+55 0.894 
cell 

6D water vs summer DETER OG eae} —-326+19 0.968 
temperature 

6D water vs relative 2384111 0 -241+78 0445 
humidity 

pO model vs. 6°70 O98 207.05 ip [oie Ayia 
model 

6Orce lulose vs 0.38+0.08 Bale L 0.788 
Savon water 

SO cellulose Vs 0.54+0.12 15832207 0.792 

summer temperature 

6° Orcellulose vs -1.7+0.8 2842545 -0.066 
relative humidity 

i Oowater vs summer I 39208) -41.642.3 02967 
temperature 


61°80 water vs 
relative humidity 29.3214 —31,129..7 Os 527 
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between 6D cell measured and 6D cell calculated by the 
model. The correlation coefficient is 0.998. The large 
correlation coefficients for both oxygen and hydrogen 
models indicate that the assumptions made in this modelling 
are generally true. 

To further illustrate the role of individual para- 
meters in the model equations, a number of additional 
relationships are given in Table 37. 6D values of cellulose 
have a strong linear relationship with that of soil water 
whereas 6!°O values of cellulose have poor correlation 
with those of soil water (see Figs. 27,28). In fact, 6D 
values of cellulose even have an excellent correlation with 
61®Q values of water (see Fig. 29). These indicate that 
6D values of cellulose mainly reflect §D values of soil 
water but 6'°O values of cellulose are influenced by other 
taccors#in. addition “og Ov values of ‘soil water. In 
support of this, 6D values of cellulose show good corre- 
lation with the growing season temperature (see Fig. 30), 
which can be attributed to the good correlation between 
6D. and temperature (see Fig. 31). Despite the fact that 
6° °O values of soil water@aére relatively well correlated 
with summer temperature (see Fig. 32), the S2n80 tia luestic£ 
cellulose are less sensitive to the growing season tempera- 
ture (see Fig. 33). There are a number of reasons for 
this poor correlation. Hirst, the contribution or the 
6189 value of water to the 6'°0O value of cellulose is only 


62% whereas that of hydrogen is 100%. Second, the CO.-H.0O 
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Figure 27. The relationship between $D values of cellulose 


nitrate in trees from various sitessin North 
America and 6D values of meteoric water used by 


the plants. 
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Figure 29. The relationship between 6D values of 
cellulose nitrate in trees from various 
sites in North America and 6!°O values of 
meteoric water used by the plants. 
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Note that the slope is close to that of the 
meteoric water line. 
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The relationship between growing season tem- 
perature and 6D values of cellulose nitrate 
in trees from various localities in North 
America. 
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Figure 31. The relationship between growing season tem- 
perature and 6D values of meteoric water. 
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Figure 32. 


SUMMER TEMPERATURE (°c) 


The relationship between growing season tem- 
perature and 6'°O values of meteoric water. 
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Figure 33. 
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The relationship between growing season tem- 
perature and 6'°0O values of cellulose in tree 
rings from various localities in North America. 
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Fractionation factor has a negative correlation with the 
temperature whereas the regression coefficient between 6!°0 
value of water and temperature is positive. Finally, oxygen 
is sensitive to the kinetic effects of evapotranspiration. 
Thus 6*°O values of cellulose become less sensitive to the 
growing season temperature due to the conflicting effects 


of the growing season temperature on a 6 Orvalue 


“a ' 
CO, HO 


of water and kinetic effect of evapotranspiration. 
Fig. 34 shows the linear relationship between rela- 
tive humidity and ASD which is the enrichment of cellulose 


relative to soil water defined by 


ASD = oO m 6D. [78] 


The striking feature of Fig. 34 is the linearity between 
ASD (enrichment of cellulose with respect to soil water) 
and relative humidity. In sharp contrast to Epstein's 
finding that the deuterium enrichment of cellulose relative 
to environmental water is constant (2291/7 )o) te S appa 
rent that the enrichment of cellulose is a linear function 
of relative humidity. Thus if the differencein 6D values 
between soil water and cellulose in tree rings is known, 
the relative humidity of the growth site can be obtained 
and vice versa. 

However, it is interesting to note that both the 6D 
and the 6!80 of cellulose in tree rings are poorly corre- 
lated with relative humidity (see Table 37). This may be 


due to the fact that unlike meteoric water, relative 
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Figure 34. 


RELATIVE HUMIDITY 


The relationship between relative humidity and 
enrichment of cellulose nitrate in deuterium 
with respect to soil water. 
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humidity has no significant functional relationship with 
geographic” factors™such’astatitude~and-alLerecude,;-«which 

is partly shown in poor correlations of relative humidity 
with either the isotopic composition of soil water or with 
air temperature. For example, air temperature, in general, 
decreases as the latitude’ or altitude *increases'.in «the 
northern hemisphere. However, no such generalization is 
found with relative humidity. 

Despite the fact that oxygen isotope fractionation 
is sensitive to relative humidity, Gray et al (52) obtained 
linear relationship between 6'°0O values of cellulose and 
mean annual temperature for various terrestrial and peat 
deposits. The virtually identical slopes shown by two 
different types of plants representing different humidity 
environments were attributed to the constant relative 
humidity effect by Gray and Thompson (52). However, it 
also may be interpreted as a lack of a functional relation- 
ship between relative humidity and 60 values of cellulose 
due. to random distribution of relative humidity throughout 
North America. Thus the scatter of isotope data in Fig. 12 
may be attributed mainly to random distribution of relative 
humidity throughout the regions. 

However, when the isotopic variation of cellulose 
from a single tree is considered, the relative humidity 
plays an important role because variation in relative 
humidity at a given site appears to have a functional 


relationship with the change of air temperature at the site. 


os 


: a —- 5 


Stley giisadisales conta dean gnenttineta: on 8 sd Y9ibh 


inane ,shittisls Ss sbganeeb ‘es Rove sxosa88 staan 4 
yiibiqud aviceie? to eabpeeteeaee wom ai" “pworte YE 
itiw t6 sajg6ew ices Te Hod aegis slavparaiaadl sas cate 
fstensn at . so otetoemed zis vo LGMERs _ . Stouts rc 


_ & 


of ni essssioni shudeele zo shusiael aft ze 


nolijsexslissensp nd5ve om ,rerewon = 
ate | oetsbiauel al ch iain asiw & 
[opat2? Sqosteb names tard paek wd edsenat 
(cE) [5s 99. YESS) VSR RM, Suita oF ewk sium 
- seoiuitos 23 eduthy OMB) meahied gbilenonsaten 
: -s p- 
$e9q bos Isiatesasrset ebGzwey 205 ovegacipgeed a 
{wt yd owede peddle tRokdilobs Kt banerake sit vets 
vwschcmod aasystt ib pas ive sied einsly 3 esqy? inex 
gyiseie: anatenos Sat ae ed uclitetas —_— ataswe nontty 
(aevewoH (2) sogemont Bs verte Yd roote (abi a 
neizsiga lsmoitoact (6530 soe 5 BS betergresat od vee 
a ifep, 20 Zewiley ang Bris ke Brace meinen neswied ¢ / 
JUORoOVoGsS Veo er pitta ley, to: aebsad tel h mone os li ab 
-ph4 mi Sye5 agnsee1 + seat8oe ane audT “sukneuk sxom 
aviseias 20 AGitids fate niebriat ad ore besudiasts ed-yem 
sau: ites snd ‘svoseuetd3 x9 tbtaudt 
sSaoYiles to ott eai saqusea sit nedw \ 28VOWOH ) 
VWiEGifen: svigeres ont" ee BE esis sipnte 5 wt 
sViisis. ini nol teisayogenkoad ef6% sass soal ts aves 4 
ténofisnrt..6 SVaTt (OS, SIesyas asia mevig & 38 ves Eka 


waeih ens 46 ees the te ee anz. thee iano shin 


5 


L283 


bid Relationship between 6D Values and 6!°O Values of 
Cellulose inferred from the Model Study 


From the models developed, it is obvious that the 
relationship between 6D and 6'°O of meteoric water does not 
hold for cellulose in plants. Thus, a new relationship 
between 6D and 6'°O of cellulose is sought. The equations 


used in this calculation are as follows: 


18 Zs 18 fen 
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[80] 
SD = g6'®o = l0o4h + 69 [81] 
cell sw 
Opera les Shere 2a of 3 [82] 
SW SW 


First, the case where trees from various localities are 
involved is considered. 


Eliminating 6D, and ENEO. yields 
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Since 6D el has almost no functional relationship with 
e 

relative humidity or Chesfractionation factor "CO,-H,0' 

the humidity term and. fractionation factor can be 

regarded as constant with large random errors. If rela- 


tive humidity and daytime temperature are arbitrarily 


chosen to be 70% and 25°C respectively, the equation 


becomes 
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SD all SIO OF 432. [84] 


The linear regression equation describing §D and §!°0 of 
cellulose in trees from various sites in North America was 


calculated (see Fig. 35): 


5) Bee ee ast - 459, [35] 
The slopes of the two equations are identical within the 
error limits. Furthermore, the close agreement of the 
intercepts indicates that the models developed in this 
study do indeed account for the isotopic composition of 
cellulose in plants. 

Even if there is no systematic geographical variation 
in relative humidity, the variation of relative humidity 
for a given site appears to have a functional relationship 
with the air temperature. Thus the relationship between 
6D and §!80 values of cellulose in a series of tree rings 
from a single tree is likely different from the one obtained 
above. @Figures 36 and 37 show the diurnal and daily varia- 
tions of air temperature and relative humidity at Julich 
(37,20). From these.figures, it 1S apparent that, both the 
diurnal and daily variations: of relative humidity. are: not 
constant and must be eliminated from the equation. Elimi- 


nation of relative humidity reduces equation [81] to 
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6D abs On 29 724.756 Os 242. [86] 
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Figure 35. 
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The relationship between the §!°0O values of 
cellulose and the §D values of cellulose 
nitrate in trees from various localities in 
North America. 
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Since the variation of §!°0O of soil water for a given site 
is relatively small and approximately one-half of 6/°%o 
variation of cellulose, the equation reduces to 
— rus 18 
SD ea Gan ES) O25 shia) “aaa [87] 
In support of this result, the spruce trees from Edmonton 


have slopes. of 2.3) 7 and 5.3. which-lie within. the, range 


of the calculated slope. 


6. 6 The Best Climatic Information from Tree Rings 


It is apparent that relative humidity influences 
both the hydrogen and the oxygen isotopic composition of 
cellulose in tree rings. However, the derivaton of 
humidity TneoRme elon from isotopic data of cellulose in 
Plants appears to be unreliable and impractical. Despite 
the fact that the fractionation ASD between 6Doai1 and 
6D. is well correlated with relative humidity, 6D of cellu- 
lose has a poor correlation with relative humidity. This 
indicates that unless both 6D of cellulose and 6D of soil 
water are known, 6D of cellulose cannot provide any 
reliable information about relative humidity. In the case 
of paleoclimate reconstruction, the 6D of soil water is 
another unknown variable. Therefore, the claim that the 
isotopic composition of cellulose can be used as paleo- 


hygrometer (35) appears not to have strong support either 


from theory or from frelagdaca. 
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On the other hand, despite the fact that the hydrogen 
model does not contain a temperature-dependent fractiona- 
tion term and the temperature-dependent fractionation of 
CO. involved in the oxygen model is small, temperature 
(both the growing season and the annual) has a relatively 
good correlation with isotopic compositions of cellulose. 
In particular, temperature shows excellent linear correla- 
tion with 6D values of cellulose as well as 6D values of 
Soil water.,.In spite.of the relative humidity effects on 
sbavalues of cellulose, 6D values of cellulose ‘still show 
excellent correlation with temperature. This suggests that 
not only 6D of soil water but also the relative humidity 
has a functional relationship with temperature at the site. 

A number of researchers have shown that diurnal 
changes in relative humidity are caused by variation in 
temperature at the site and that there exists a negative 
correlation between them (37,20). The average regression 
coefficient calculated from their data turns out to be 
approximately »-0.06°97.,./2C. .Utilizing,this value, an 
attempt has been made to examine this temperature effect 
Viasrelative humiditysonndD of cellulose jin wthe model 
developed. here: Firstyorvsall, oD values wot cellulose jin 
trees from various localities are considered. Since the 
geographical relative humidity distribution is not a 
function of temperature, the overall 6D values of cellulose 
from various sites become a function of temperature only. 


This is indeed confirmed by the empirical equation 
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5 = 5.5. SFOS KS 2 [88] 


a cart 


which has the same temperature coefficient (5.5°/,,/°C) 


as that of the meteoric water equation 


— = rs a fe} 
6D éDe, 5.6t. LOU at aces [89] 
Secondly, when 6D values of cellulose from a single tree 
are considered, the negative temperature coefficient 
associated with relative humidity must be taken into con- 


Sideration. Thus, 


SD = 5) Gt aut Ge bok. [90] 
a SS} 


cell 


In the case of the oxygen model, the equation reduces to 


gS up eve gal sat sk ee (Ong ah tl [91] 
The resultant temperature coefficients would vary from 5.6 
tor il 27for the. hydrogen model and from 0.43 to 1.5 for the 
oxygen model depending upon the contribution from the 
growing season temperature terms. The dual temperature 
scales in equations [90,91] may introduce complications 

fo -aunterpretation of the 1sotope-Lemperacure relationship. 
In general, en and t. vary in the same direction reinforc- 
ing temperature coefficients. This phenomenon accounts 
for the increased temperature coefficients of both the 
hydrogen and oxygen isotopic compositions of the three 


spruce samples from Edmonton (see Tables 19 and 27). The 
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Measured temperature coefficient of hydrogen equations 

for WSL)®@WS2"'and? WS3s-are® 723713" and@ 415 respectively. The 
oxygen temperature coefficients for the equivalent samples 
are” 1.36, 1.03 and° 1.47 respectively. 

It is apparent that §D values of cellulose in tree 
rings reflect mainly 6D values of soil water. Superimposed 
upon that, temperature effects via relative humidity also 
play an important role in determining 6D values of cellu- 
lose. The dual temperature scales involved in equations 
[90] and [91] may sometimes result in scattering of data 
if the relative humidity has no functional relationship 
with temperature. But the high temperature coefficients 
Shown by individual trees from Edmonton indicate that 
temperature effect via relative humidity is generally in 
phase with temperature effect via hydrogen isotope frac- 
tionation of soil water and thereby increases the overall 
temperature coefficient. 

In spite of the dual temperature scales, §D values 
of cellulose can be generally interpreted in terms of 
temperature. Therefore temperature is the best information 
one can derive from isotopic analysis of cellulose in 
Plant materials. ‘The signiticance of equations. )90)] sand 
[917] is that, unlike relative humidity, temperature infor- 
mation can be obtained without any prior knowledge of the 
isotopic composition of the soil water. This is particu- 
larly important in the case of paleoclimatic reconstruction 


because the isotopic composition of the soil water for a 
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given site is not known. 

The results sof) this study suggest that cellulose in 
tree rings contain isotope thermometers (D and !%0) whereby 
paleotemperature variation of the environment can be 
reconstructed. However, care should be taken in interpre- 
tation of the isotope data. The isotope data obtained 
from a tree may not be the representative of other trees. 
For example, WS1l and WS2 exhibited systematic difference 
in 6D values of cellulose for the time equivalent rings 
despite the fact that both trees came from practically the 
same site. This suggests that there may exist biological 
fractionation in isotopic composition between individuals 
as well as species. When paleotemperature reconstruction 
is sought, it is desirable to analyze sequences of tree 
rings from a single tree. When isotope data from more 
than one tree are combined, it must be ascertained that 
the two trees are overlapped for some growth period and 
that the time equivalent tree rings have similar isotopic 


compositions or corrections must be made. 


ray PE 


(o** ‘Bre oes siieiba 
ed aso jnaoncmrma Sf Geogoeeeeaey, 
ad bideth Stee -, revewal : 

Jarek saomags wid 20 


és wearw 


cere 


etgisint mh .ae 
hoagiiesde astat aaannas any 
@eots usigo fo sv istsinaeesmaes SHY sd Som, gam set 5 a 
sone tet ih oinbee: saver B pStidbdxe | — atl nial » alge xe 
‘ppre> treiay sues: Site ris: yes eatery So seutsy a 
and \YCieoijos td moss smss aay? tod send goat arid ediqn 
laptpoicid Jains -yeo eee SWNy Sstegpee. aceP. adiiee-s: 
eleubivyires dewaen 19is aoa 2itorees ai nodssaotae. 
fi Misasiccs suse iagmgae tide. cite ei nega aS i 
[e0rpSss sav iene. oF sities ei +i <dipcom: at 
9909 sharihe ® moms ot 


saad iol ess) 
 sgodtoes pany 

6 40 Jab ta Reais HES ous Seo | 

‘ serie t dipwre: ets eoo33. ow: 


$eit4d | Renal 
bOrte¢ rswott) Gane -to7 
7a sash Reds sosibeyiupe emis ata 3 
chan ad abit ecebebegngo to onois Leeda 


i< 
= 


CHARTER: *¥ 


SUMMARY AND CONCLUSIONS 


A new method of measuring both D/H ratios of !°0/!*o 
ratios in organic compounds has been developed. 

This technique involves pyrolysis of organic compounds 
in a nickel reaction vessel. The pyrolysis method 
gives reproducible and accurate results with a variety 


of organic compounds and standard waters. The aver- 


(oe) 


age precision of hydrogen isotopic analysis is +23. 
At present the method is limited to the organic 
compounds whose carbon-oxygen ratio is less than l, 
due to formation of tar during early stage of pyro- 
lysis: However, «introduction Of an voxidizindyagens 

or free oxygen into the pyrolysis may remove this 
limitation. 

Use of the nitration technique eliminates the 
exchangeable OH hydrogen without affecting the 
original isotopic composition of non-exchangeable CH 
hydrogen in cellulose. Pyrolysis of cellulose nitrate 
in the nickel reaction vessel shows no indication of 
tar formation indicating completion of reaction. 
Nitration of cellulose extracted from wood by the 
sodium chlorite method is superior to the direct 
nitration of wood followed by acetone dissolution in 
terms of reproducibility and correlation with clima- 


tic parameters. The direct nitration of wood 
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produces cellulose nitrate whose 6D pants rs 
Systematically lower than that of nitrated cellulose 
extracted by sodium “chlorite method: 

6D values of cellulose nitrate in a long sequence of 
tree rings from single trees contain two different 
isotope scales. All three spruce samples WS1, WS2 
and WS3 from Edmonton clearly show a sudden shift in 
6D values between the juvenile and mature sections of 
tree rings. In all cases, juvenile sections are 
systematically depleted in deuterium compared to 
mature sections. The reason for this is not com- 
pletely understood. It may be attributed to unknown 
hormonal or biological metabolism associated with the 
tree growth, with results in changes in relative 
amount and structure of chemical constituents in 

tree rings. Such changes in turn may give rise to 
large variation of D value. When paleoclimate recon- 
struction is sought, the juvenile sections of tree 
ring sequence should be eliminated from the analysis. 
The differentiation of juvenile section from mature 
one is possible using ring width data. 

The isotopic composition of cellulose in tree rings 
preserves a functional relationship with climatic 
parameters. ‘6D ‘values’ in”cellulose *nitrate™ftrom'a 
long sequence of tree rings show a strong linear 
correlation with temperature at the growth site. On 


the other hand, 6'°O values of cellulose are found to 
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to be less well correlated with temperature. The 
choice of mean annual temperature is rather subject- 
ive. The mean annual temperature calculated from 
September of 2 years prior to the growth year to 
August of the previous year shows the best correla- 
tion with isotope data of tree rings. In other words, 
it is the mean annual temperature of the previous 
year which shows the best correlation with both 6D 
and 6 | °O.values.o+ tree rings. 

For the three spruce samples analyzed, winter 
temperature plays a more important role in determin- 
ing, 6D-and 6!20 of cellulose. than the temperature, of 
the growing season. This conclusion together with 
the previous one indicates that the isotopic compo- 
sition of soil water is the primary factor control- 
ling the isotopic composition of tree rings and that 
soil water taken by the trees consists mainly of 
winters precipitation. This conclusion. 1s, further 
supported by the fact that the amount of precipitation 
shows a negative correlation with the isotope data. 
The negative correlation in turn is due to the amount 
effect of meteoric water on the isotopic compositions 
of* precipitation? P However," the® preferentialyassoc= 
iation of the isotopic composition with a particular 
seasonal temperature may be a growth site effect that 
is controlled by factors such as lithology and 


permeability of soil strate, and the topography of 
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the site. 

Among the various chemical constituents in wood, 
cellulose shows the best correlation with temperature. 
oD values of extractive free wood show no significant 
correlation with temperature. 6D values of lignin 
calculated from the isotope data of cellulose and 
extractive free wood also show poor correlation with 
temperature. Hydrogen isotope fractionation between 
lignin and cellulose shows poor correlation with 
temperature indicating the absence of a thermodyna- 
mic thermometer in the lignin-cellulose pair. 

Ring width data show no significant correlation with 
isotope data or with mean annual temperature. 
Hydrogen and oxygen isotope models have been 
developed to examine the role of individual climatic 
parameters in determining the isotopic composition 

of tree rings. The model equations describing the 


isotopic compositions of tree rings are given by 


= * + (6 -e )h+E 91 

SD ool] SD te +e, (O. ey) [91] 
18 - 3 18Q t+e*+te, +(a_-e, /h}+E 

65° 0.611 = 620 maseoperieie) 6215 °°O_ te*te, + (a, €,) } 921 


The models developed here are in excellent agreement 
with the isotope data. The working model equations 
which describe most of the natural isotope signature 


in trees is given by 
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6D = 6D. - 104h + 69 [93] 


6°°0 61) = 620 In Co.-H.o + 9-62(6'%0_ +33.6-32h)+2.5 [94] 
2ir2 


It is apparent that 6D values of non-exchangeable 
hydrogen in cellulose of tree rings reflect mainly 

6D values of soil water. However, in addition to 

this growing season temperature via relative humidity 
also plays an important role. There is no indica- 
tion of temperature-dependent hydrogen isotopic 
fractionation occurring during photosynthesis. 

The oxygen in cellulose appears to originate from 

CO, which is in complete equilibrium with leaf water. 
Thus the oxygen model contains temperature-dependent 
fractionation factors. Due to the conflicting effects 
of temperature via CO,-H,0 fractevonatirom, isoil 
water-vapor fractionation and relative humidity, 
oxygen shows a poorer correlation with temperature 
than does hydrogen. 

The hydrogen isotope enrichment of cellulose with 
respect to soil water is a linear function of relative 
humidity at the site. Thus, relative humidity of a 
given site can be estimated by measuring the 6D of 


tree rings and the 6D of water taken by the tree. 


They are related by the equation 


a = = - + 69 95 
AéD 6D 6Do 104h [95] 
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The relationship between 6D and 6§!80O values in 

cellulose varies depending upon the association of 
relative humidity with temperature. When relative 
humidity is linearly related to temperature at the 


site the relationship is given by 


6 (5.242)6!80 + RB. [96] 


een peu 


This is generally the case where isotope data are 
obtained from a single tree. When relative humidity 
has no such functional relationship with temperature 


this becomes 


= 18 u 
BE aie lye ee 2 [97] 


The isotope data collected from a variety of growth 
Sites belong to this category. 

The best information one can obtain from isotope data 
of tree rings is mean annual temperature. However, 
interpretation may be complicated by the effect of 
relative humidity. The relationship between 6D cell 


and temperature is given by 


-4.6t +6.2t +E [938] 


ee tl a s 


where the term 6.2t is the contribution from relative 
humidity. In general, the changes in the growing 
season and mean annual temperature are in the same 
direction. Thus trees from single site will have a 


temperature coefficient bigger than 5.6. But mean 
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6D values of tree rings from various sites will have 
a value! close M645 76. 

The fact that the two trees from the same locality 
showed systematic difference in 6D values of cellulose 
indicates that absolute temperature scale is not 
guaranteed from isotope data. Nevertheless, the 
temperature variation can be obtained from D/H 

ratios of non-exchangeable hydrogen in cellulose in 

a sequence of tree rings. When paleoclimate 
reconstruction is desired, hydrogen isotopic composi- 
tion of cellulose nitrate ina long sequence of tree 
rings from a single tree should be analyzed. When 
isotope data from more than one tree are combined, 
the correction has to be made in such a way that the 
time equivalent ring section have the same isotopic 


composition. 
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RECOMMENDATION 

The ultimate purpose of iberene studies in climato- 
logy is to obtain paleo-temperature variation over the last 
few thousand years. The information will form a basis 
which enables a better understanding of the natural 
mechanisms of climatic variation to be obtained and thereby 
future climate predictions made possible. The present 
study clearly shows that isotope records in tree rings can 
provide such information. Therefore, detailed isotopic 
CELE EO in a long sequence of tree rings should 
be measured. Spectral analysis of isotope data thus 
obtained will show the cyclic nature of climatic variations. 
The sampling interval of tree rings, however, should be 
chosen in such a way that the desired periodicity can be 
resolved. 

Present isotope technique provides only temperature 
variation. As for the absolute temperature scale, it is 
apparent that further work must be done. Absolute tempera- 
tures can be determined only from two co-existing substances 
formed in equilibrium with each other. Such systems may 
be found in cellulose. There are two types of hydrogen in 
cellulose. One is carbon-bound hydrogen and the other is 
oxygen-bound hydrogen. Carbon bound hydrogen is stable 
whereas oxygen bound hydrogen exchanges isotopically with 
water. However, about 25% of hydroxyl hydrogen which is in 


crystalline regions is stable. Thus if a significant 
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Bee eicnation exists between carbon bound hydrogen and 
non-exchangeable hydrogen in crystalline regions, an 
absolute temperature scale may be obtained by measuring 
the isotopic composition of CH hydrogen by the nitration 


method and OH hydrogen by the equilibration method. 
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Hydrogen Data (SMOW) 
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APPENDIX. 2 


Oxygen Data (SMOW) 


Sample year group Temp°C 61°8O(SMOW) 
1 WSL 1883-1888 Le 36 23-66 
2 1888-1893 2.60 24.13 
3 1893-1898 Leak 24.20 
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APPENDIX 3 
Direct Nitration of Wood (W-C-N) 


WS1 (#0 #17) 


6D (SMOW) Temp°C 

1} -146 Pe eiall 
Zz -139 3.41 
3 -142 Sells 
4 -150 2.04 
5 -137 2.89 
6 =<137 CT ae | 
"3 -148 24 
8 -139 314 
9 -143 VL ete) 
10 -142 Deel 
Epps -148 ee Sie) 
1 -141 PLAT Rs: 
eS -142 Sy Aad) 
14 -157 PARASITE 
LS -162 Jo 4 
6 -160 25.100 
Lid -171 E06 


ASD = +2.4 (Average) 


Ze 


ae 


; : ; f hoe 7 
E 4EGMRSTA tae 7 7 i" , 
(Die i hone . 706 de ; ; ‘ ; ph : _ rf y 


a" , 

oe. 

Srgnat . — (MOMEGa 
=3) 


3 
>> Di wy 

bow 3 

b 


; 7  fAl~= 


=. 
~ ’ 
io OT oe re or > ee : 
a 
- 
La) 
s 


Nie t2l- 
be Sat- 
4 ~ Rik 4 
ae ifi- 


= 


(apesevay &. my = aang 


- 4 
~ 
SO 
"fy 
i! , 
a 
’ =! y 


raat A : 
aia 
oe 


roe 


